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TECHNICAL  REPORT  SUMMARY 


Sections  A-C  describe  the  research  forming  the  Ph.D  dissertation  of  William  R.  Walter 
which  was  supported  in  part  by  contract  F19628-89-K-0022.  Section  A  (Walter  and  Priest¬ 
ley,  1991)  reviews  two  of  the  most  widely  used  P  wave  spectral  models  of  earthquakes  and 
two  leading  models  of  explosions,  to  assess  their  implications  for  high  frequency  (>  1  Hz) 
seismic  discrimination  of  underground  nuclear  explosions  from  earthquakes.  The  models 
are  then  compared  with  North  American  events  and  a  limited  data  set  of  Central 

Asian  events.  The  earthquakes  are  consistent  with  a  constant  stress  drop  model  which  de¬ 
cays  as  at  high  frequencies.  The  deeper  (>  700  m)  Nevada  Test  Site  (NTS)  explosions 
also  show  a  fall  off  between  10  and  30  Hz,  in  contrast  with  the  shallower  explosions 
which  have  a  steeper  decay.  Near  regional  recordings  of  the  September  14,  1988  Soviet 
JVE  explosion  show  a  higher  corner  frequency  and  lower  1  to  4  Hz  P  wave  spectral  ratios 
than  predicted  by  either  explosion  model,  and  than  predicted  for  similar  size  earthquakes. 
This  high  corner  frequency  may  be  influenced  by  source  complications  such  as  spall  and 
tectonic  release 

Section  B  (Walter  and  Patton,  1990)  describes  some  simple  forward  modeling  to  quan¬ 
tify  both  the  moment  and  the  mechanism  of  the  tectonic  release  of  the  Soviet  JVE.  Using 
10-20  second  Love  and  Rayleigh  wave  data  the  tectonic  release  moment  is  found  to  be 
from  one  fifth  to  one  tenth  of  the  explosion  moment.  To  resolve  the  ambiguity  in  tectonic 
mechanism,  synthetics  seismograms  were  computed  to  compare  the  relative  amplitudes  of 
the  SH  and  SV  parts  of  the  SmS  phase.  Although  the  character  of  the  SmS  phases  are  not 
completely  matched,  the  strike  slip  model  proposed  for  the  tectonic  release  comes  closer  to 
matching  the  observed  relative  amplitudes  than  the  possible  reverse  faulting  mechanism. 

Section  C  (Walter  and  Brune,  1992)  describes  a  spectral  model  for  the  seismic  radiation 
of  a  circular  crack  failing  in  tension.  The  model  may  be  a  useful  first  approximation  for 
modeling  seismic  sources  with  some  volumetric  component,  such  as  magma  injection,  spall, 
or  slip  on  a  non-planar  fault. 

Section  D  (Priestley  et  al,  1992)  describes  a  study  of  the  tectonics  of  the  region  sur¬ 
rounding  the  South  Caspian  basin,  a  relatively  aseismic  block  within  the  Alpine-Himalayan 
Belt,  which  is  surrounded  by  zones  of  high  seismicity.  Earthqucike  focal  mechanisms  and 
depths  were  determined  for  16  earthquakes  in  the  region  and  these  used  to  investigate  the 
active  tectonics.  Earthquakes  SW  of  the  Caspian  Sea  have  shallow  thrust  mechanisms 
showing  that  the  continental  crust  of  NW  Iran  is  overthrusting  the  "oceanic-like”  crust 
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of  the  south  Caspian  Basin.  Both  high  angle  reverse  faulting  and  left-lateral  strike  slip 
faulting  are  observed  south  of  the  Caspian  Sea  in  northern  Iran.  These  mechanisms  show 
oblique,  left-lateral  shortening  and  overthrusting  of  the  south  Caspian  Basin  by  north- 
central  Iran,  but  that  this  oblique  regional  convergence  is  partitioned  into  pure  strike-slip 
and  reverse  faulting.  Earthquakes  e«ist  of  the  Caspian  Sea  also  show  a  mixture  of  high 
angle  reverse  and  strike-slip  faulting  mechanisms  and  thus  may  be  another  excimple  of  the 
partitioning  of  oblique  slip  into  strike-slip  and  thrust  motion.  Normal  faulting  mechanisms 
dominate  in  the  belt  of  seismicity  which  extends  across  the  central  Caspian  Sea.  It  is  im¬ 
probable  that  these  events  represent  the  motion  between  the  southern  Caspian  Basin  and 
Eurasia  for  they  imply  a  sense  of  motion  which  is  incompatible  with  the  observed  topogra¬ 
phy  and  folding  in  the  sediments.  Two  shallow  earthquakes  in  this  belt,  one  a  small  event 
and  the  other  a  large  second  event  of  a  multiple  earthquake,  have  thrusting  mechanisms 
suggesting  that  shortening  occurs  as  the  continental  crust  of  the  northern  Caspian  is  thrust 
over  the  ” oceanic-like”  crust  of  the  southern  Caspicin  Basin.  Shortening  is  also  suggested 
by  the  orientations  of  folds  in  the  sedimentary  cover  south  of  the  central  Caspian  Sea  seis¬ 
mic  belt.  The  normal  faulting  events  appear  to  occur  in  the  deep  basement  of  the  material 
being  overthrust.  We  suggest  that  this  shortening  does  indeed  represent  a  NNE  motion 
of  the  Caspian  Sea  relative  to  Eurasia,  but  that  the  motion  is  slow  and  has  not  produced 
many  earthquakes.  The  '•ombined  effect  of  the  overthrusting  of  the  south  Caspian  Basin 
in  the  southwest  and  south,  and  north  will  lead  to  the  eventual  destruction  of  the  south 
Caspian  Basin  and  the  possible  formation  of  an  intermediate  depth,  dipping  seismic  zone 
within  the  continental  interior,  similar  to  that  presently  observed  in  the  Hindu  Kush. 
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SECTION  A 


High-Frequency  P  Wave  Spectra  from 
Explosions  and  Earthquakes 

William  R.  Walter  and  Keith  F.  Priestley 
Seismolcgical  Laboratory,  Mackay  School  of  Mines,  University  of  Nevada,  Renu 


Abstract 

Two  explosion  P  wave  spectral  models  (Sharpe,  1942;  Mueller-Murphy,  1971) 
and  two  earthquake  P  wave  spectral  models  (Archambeau,  1968,1972;  modified  Brune 
1970,1971)  are  reviewed  to  assess  their  implications  for  high  frequency  (>  1  Hz) 
seismic  discrimination  between  earthquakes  and  explosions.  The  importance  of  the 
comer  frequency  scaling,  particularly  for  models  with  the  same  high  frequency  spectral 
decay  rate,  is  demonstrated  by  calculating  source  spectral  ratios  (a  potentially  impor¬ 
tant  regional  discriminant)  for  these  models.  We  compare  North  American  events  and 
a  limited  data  set  of  Central  .Xsian  events  with  these  spectral  models.  We  find  North 
American  earthquakes  are  consistent  with  a  constant  stress  drop  modified  Brune  model 
between  10  and  30  Hz.  Shallow  (<  7(K)  m  depth)  Pahute  Mesa  explosions  at  the 
Nevada  Test  Site  have  a  high  frequency  spectral  decay  between  10  and  30  Hz  greater 
than  the  predicted  by  the  explosion  models.  Near  regional  recordings  of  the  Soviet 
Joint  Verification  Experiment  (JVE)  explosion  show  a  higher  corner  frequency  and 
lower  1  to  4  Hz  spectral  ratios  than  predicted  by  either  explosion  model.  The  higher 
corner  frequency  of  the  Soviet  JVE  appears  not  to  be  due  to  attenuation,  or  receiver 
effects,  and  may  represent  a  need  for  different  corner  frequency  scaling,  or  result  from 
source  complications  such  as  spall  and  tectonic  release.  A  regional  recording  of  the 
Soviet  JVE  (NEIC  mj,  =  6. 1)  is  shown  to  have  a  lower  1  to  4  Hz  spectral  ratio  than  a 
smaller  earthquake  (NEIC  m(,  =  4.6)  recorded  on  a  nearly  recipnxal  path. 


Introduction 

The  United  States  and  the  Soviet  Union  have  been  observing  the  terms  of  the 
Threshold  Test  Ban  Treaty  (TTBT)  and  the  Peaceful  Nuclear  Explosion  Treaty  (PNET) 
since  1976.  These  treaties  limit  the  maximum  yield  of  underground  nuclear  testing  to 
150  kilotons  (kt).  Large  explosions  (>20kt)  can  be  readily  detected,  and  discriminated 
from  earthquakes  with  seismic  waves  recorded  at  teleseismic  (>3(XX)  km)  distances  (cf. 
Dahlman  and  Israelson,  1978).  Smaller  explosions  detonated  under  these  treaties  or 
under  a  future  Comprehensive  Test  Ban  Treaty  (CTBT)  or  lower  yield  TTBT  require 
detection  and  discrimination  using  regional  (<3(XX)  km)  seismic  data.  Regional 
waveforms  offer  the  potential  of  increased  signal  bandwidth  for  detection  and 
identification,  particularly  at  high  (>  1  Hz)  frequencies.  Various  studies  (i.e.  Pomeroy 
et  al.  1982;  Bennett  and  Murphy,  1986;  and  Taylor  et  al.,  1989)  suggest  that  spectral 
ratios  of  some  regional  phases  hold  promise  for  discrimination.  However,  as  Taylor  et 
al.  (1989)  note,  "...  a  number  of  uncertainties  regarding  the  lack  of  a  physical 
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understanding  of  how  the  spectral  discriminants  work  makes  their  utilization  question¬ 
able  at  this  point".  Under  a  CTBT  or  lower  yield  TTBT  there  may  be  a  need  to  detect 
and  identify  events  from  a  geographic  region  in  which  we  have  little  or  no  previous 
experience.  In  addition  it  will  be  necessary  to  understand  under  what  circumstances  a 
particular  detection  and  identification  method  might  fail.  For  these  reasons  it  is  impor¬ 
tant  to  develop  a  physical  understanding  of  both  the  source  and  the  propagation 
characteristics  of  regional  waves. 

Separating  source  and  propagation  effects  can  be  difficult,  particularly  for  high 
frequency  waves.  The  observed  seismic  displacement  amplitude  spectra  Q(to)  can  be 
considered  to  be  the  product  of  a  series  of  linear  filters: 

Q  (CO)  =  S(o))  A(co)  G(co)  1(0))  (1.1) 


where  S(0))  is  the  source  spectrum,  A(co)  is  the  attenuation  operator,  G(0))  is  the 
Earth’s  response,  l(co)  is  the  instrument  response,  and  co  =  27tf  is  the  angular  fre¬ 
quency.  While  1(0))  is  easily  determined,  any  attempt  to  uniquely  isolate  S(0))  is  sub¬ 
ject  to  limitations  in  our  knowledge  of  A(o))  and  G(o)).  The  attenuation  operator  is 
given  as  A(0))  =  e-“‘'2Q(w)  j  travel  time  and  Q(o))  is  a  dimensionless 

measure  of  the  energy  lost  in  each  sei.smic  wavelength.  In  most  regions  Q(o))  has  been 
found  to  be  a  constant  or  a  slowly  increasing  function  of  frequency.  At  higher  fre¬ 
quencies  the  effect  of  the  attenuation  operator  becomes  increasingly  important,  strongly 
trading  off  with  estimates  of  the  source  spectral  amplitude  (cf.  Anderson,  1986).  Many 
high  frequency  (  >  1  Hz)  studies  of  P  wave  attenuation  assume  a  particular  form  for 
the  source  spectrum  and  determine  the  attenuation  by  adjusting  Q(o))  until  the 
observed  spectrum  matches  the  assumed  source  spectrum  (e.g.  Der  et  al.,  1985,  Hough 
et  al.,  1988,  Sereno  et  al.,  1988).  The  Earth  response  G(a))  is  also  subject  to  more 
uncertainty  at  high  frequencies  since  shorter  wavelength  signals  are  affected  by  small 
scale  features  in  the  Earth.  G((o)  includes  source  region  effects  such  as  free  surface 
reflections,  geometrical  spreading  which  may  be  frequency  dependent  for  some 
regional  pha.ses  (e.g.  Sereno,  1990),  scattering,  and  receiver  site  effects.  Recent 
borehole  studies  at  several  hardrock  .site;s  have  found  that  attenuation  and  re.sonance 
within  a  shallow  weathered  zone  (typically  tens  of  meters  thick)  can  have  a  pro¬ 
nounced  effect  on  observed  high  frequency  P  wave  spectra  (Malin  et  al.,  1988;  L.  Car- 
roll,  J.  Fletcher,  H.  Liu,  and  R.  Porcella,  unpublished  data,  1990) 

As  a  result  of  the  difficulties  of  uniquely  determining  the  source  spectral  ampli¬ 
tude,  various  models  for  both  earthquakes  and  explo.sions  have  been  proposed,  each 
with  different  implications  for  the  problems  of  seismic  discrimination.  For  example, 
the  Sharpe  (1942)  and  Mueller-Murphy  (1971)  explosion  source  spectral  models  are 
richer  in  high  frequency  P  waves  than  the  Archambeau  (1968,  1972)  earthquake  source 
spectral  model,  when  normalized  to  have  the  same  low  frequency  level.  This  was  the 
basis  for  a  proposal  by  Evernden  et  al.,  (1986)  that  the  solution  to  the  problem  of 
detection  and  identification  of  low-yield  underground  nuclear  explosions  was  available 
using  high  frequency  P  waves.  However,  if  the  Brune  (1970)  earthquake  source  spec¬ 
tral  model  modified  for  P  waves  (Molnar  et  al.,  1973)  is  compared  with  these 


explosion  models  the  source  spectral  differences  are  much  smaller. 

In  this  paper  we  will  first  review  these  four  competing  spectral  models  and  their 
implications  for  the  seismic  discrimination  issue,  then  discuss  efforts  to  distinguish 
between  various  earthquake  and  explosion  P  wave  source  spectral  models  using  data 
from  western  North  America.  Finally  we  will  compare  those  results  with  results  from 
near  regional  (<  800  km)  seismic  data  collected  in  the  vicinity  of  the  Soviet  Union’s 
East  Kazakhstan  Test  Site  (KTS).  Although  data  for  Central  Asia  is  limited,  we  find 
significant  differences  in  high  frequency  (  >lHz  )  P  wave  spectra  for  Soviet  explo¬ 
sions  compared  with  both  shallow  explosion  spectra  observed  at  the  Nevada  Test  Site 
(NTS)  and  from  explosion  spectral  models  based  on  North  American  experience. 


Theoretical  Models 

Simple  source  spectra  are  typically  characterized  by  three  values:  the  low  fre¬ 
quency  asymptote  who.se  level  is  proportional  to  the  seismic  moment  (or  source 
strength),  the  high  frequency  asymptote  whose  slope  depends  on  the  source  displace¬ 
ment  time  history,  and  the  intersection  of  these  two  asymptotes,  the  comer  frequency, 
f(.,  which  is  related  to  the  source  dimension.  On  a  plot  of  log  displacement  amplitude 
versus  log  frequency  the  low  frequency  asymptote  has  a  constant  value.  The  high  fre¬ 
quency  asymptote  is  proportional  to  where  must  be  greater  than  1.5  to  have 
finite  energy,  A  larger  r]  indicates  a  greater  high  frequency  decay.  Our  discussion  of 
spectra  for  various  source  models  will  be  in  temis  of  these  three  parameters. 

The  generation  of  seismic  waves  by  an  underground  nuclear  explosion  is  a  com¬ 
plex  phenomena.  The  explosion  region  can  be  conceptually  modeled  with  increasing 
distance  from  the  shot  point  as  a  vaporization  cavity,  a  strongly  non-linear  region 
characterized  by  a  shock  wave,  and  a  linear  region  (Rodean,  1971).  Sharpe  (1942) 
nuxieled  the  seismic  waves  from  explosions  as  an  arbitrary  pressure  function  applied  to 
the  interior  of  a  spherical  surface  at  an  "elastic  radius",  R(,],  defined  as  the  distance 
from  the  detonation  point  to  the  point  at  which  the  medium  begins  to  respond  linearly 
felastically).  In  the  Sharpe  model  all  of  the  effects  of  the  non-linear  region:  rock 
vaporization,  cavity  expansion,  shock  wave  propagation  and  decay  to  an  elastic  wave, 
are  contained  in  the  fomt  of  the  pressure  function.  Sharpe  (1942)  discussed  several 
pressure  histories,  but  we  will  follow  Evemden  et  al.  (1986),  and  use  a  step  or  Heavi¬ 
side  function  in  pressure  acting  at  the  elastic  radius  as  the  "Sharpe"  model.  The 
Sharpe  far-fiek!  displacement  amplitude  spectrum  is  given  by; 


^  (u"-().75R,,’oj2)2  +  (aR,,to)'J 

where  is  the  amplitude  of  the  pressure  step  function,  p.  is  the  modulus  of  rigidity,  a 
the  7’  wave  velocity,  and  R  is  the  hypocentral  distance.  Evernden  et  al.  (1986)  give  a 
relationship  between  the  elastic  radius  and  the  yield  as: 


(CO) 
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R..,  =  1.61  X  10^ 


meters, 


where  is  in  pascals,  k  is  a  coupling  constant  and  W  is  the  explosion  yield  in  kilo- 
tons.  Evemden  et  al.,  (1986)  give  values  of  P„  =  20  MPa  (200  bars),  k  =  9  for  KTS 
and  Po  =  13  MPa  (130  bars),  k  =9.4  for  NTS. 

Mueller  and  Murphy  (1971)  modeled  explosions  with  a  different  pressure  function 
at  the  elastic  radius.  They  used  a  pressure  function  based  on  free-field  obsersations  of 
several  underground  nuclear  explosions.  This  pressure  function  has  the  form; 

P(t)  =  ( P,  e-Y‘ -  P2 )  H(t)  (1.4) 

where  Pj,  P2,  and  y  are  functions  of  depth  and  the  shot  point  geology  which  have  been 
empirically  determined  from  U.S.  explosions.  The  far-field  amplitude  spectnam  of  this 
model  is  given  by: 


2  I  T  D  D-  D  2 1 


„  R,,"  P,"  +  2Pip2  P2 


4  p  R  [  +(£)2 


_ _ 

(a2-().75R\ia)2)2  +  (aR^,w) 


If  we  assume  a  depth  -  yield  scaling  relationship  appropriate  for  the  NTS  of  h  =  120 
meters  and  a  granitic  medium,  then  the  empirically  determined  parameters  for 
this  model  are  (Stevens  and  Day,  1985): 

R„i  =  2.38  X  10^  „  .r  meters,  (1.6) 


Y  =  1.09x  lO"*  — ^  hertz. 


P,  =  3.77:;  lO^h  -  P2  pascals. 


Y^O.87 

1^2  3.94  X  10*'’ — z—rpr  pascals.  (1.9) 

R.i^h*^ 

Other  models  for  explosions  have  been  ba.sed  on  fitting  polynomials  to  near-field 
reduced  displacement  potentials  (Haskell.  1967.  von  Seggern  and  Blandford  1972, 
Helmberger  and  Hadley  1981).  The  revision  of  Haskell’s  (1967)  mcxlel  proposed  by 
von  Seggern  and  Blandford  (1972)  gives  speeir.i  th.it  are  nearly  identical  to  the 
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Mueller-Murphy  (1971)  mcxlcl. 

The  generation  of  seismic  waves  by  earthquakes  has  been  the  subject  of  many 
kinematic  and  dynamic  modeling  efforts.  One  the  most  widely  used  models  of  the 
seismic  source  spectra  is  characterized  by  a  single  c''>rner  frequency  and  a  high  fre¬ 
quency  decay  proportional  to  U)”‘  (Aki,  1967,  Brune  1970).  In  this  paper  we  examine  a 
specific  model,  proposed  by  Brune  (1970),  in  which  the  earthquake  is  modeled  as 
a  circular  shear  crack  that  slips  instantaneously.  Here  we  assume  that  the  P/S  comer 
frequency  ratio  is  equal  to  the  P/S  wave  velocity  ratio  as  some  studies  show  (e.g. 
Molnar,  1973).  The  far-tield  displacement  amplitude  spectrum  for  this  model  is  given 
by  (Brune,  1970,  1971): 


raa 


0^a»  =  Re4.— 


P  0)^  + 


(1.10) 


In  (1.10)  a  is  the  stress  drop,  r  is  the  radius  of  the  shear  crack,  R0^  is  the  radiation 
pattern  factor,  and  (0^,  is  the  angular  corner  frequency  which  is  given  by  0)(.  =  2.34  a /r. 
We  will  refer  to  this  model  as  W2P  (co“‘  decay  for  P  waves). 

Archambeau  (1968,  1972)  developed  a  model  for  earthquakes  in  which  the 
sei.smic  waves  are  generated  by  the  relaxation  of  pre-existing  stress  over  a  spheroidal 
volume.  The  high  frequency  decay  for  the  far-field  P  wave  displacement  spectrum  of 
this  model  (with  rupture  velocity  equal  or  less  than  the  shear  wave  velocity)  is  propor¬ 
tional  to  co”^.  We  will  refer  to  this  model  as  W3P.  We  follow  Evernden  et  al.  (1986), 
in  using  only  the  asymptotic  values  for  the  far-held  spectra  of  the  W3P  model.  Al  low 
frequencies  the  spectra  has  a  constant  low  freciuency  value  given  by: 


5  o  Peo 
'Tlpr/R 


(1.11) 


from  zero  frequency  up  to  a  comer  •recpiency: 


1 

27iL 


( 3  a’  Vj) 


1/3 


(1.12) 


after  which  the  spectrum  decays  as  In  equations  (1.11)  and  (1.12)  L  is  the  max¬ 
imum  rupture  dimension  and  v,  is  the  mpiure  vekx'ity.  In  deriving  equation  (1.11)  we 
have  assur  d  a  Ik'isson  ratio  <4  1/4  and  set  the  W3P  moment  equal  to  the  expression 
for  the  zc!  :rec|uency  level  of  a  double-couple  point  source  (Walter  et  al,,  1988). 

Making  use  of  cqu;itions  (1.2  1.12)  we  can  cotnpare  the  high  frequency 

beh;tvior  of  these  models,  l  igure  la  shows  the  theoretical  far-held  amplitude  spectra  of 
the  .Sharpe  and  Mueller-Murphy  (MM)  models  calculated  for  explosions  with  a  yield 
of  !(K)  kt  and  1  kt.  For  comparison  the  two  earthquake  models  have  al.so  been  calcu¬ 
lated  for  the  same  low  fre(.|uency  values  as  the  Sharpe  model.  For  the  W3P  mrxlel  we 
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Fig.  1.  A  compari,son  of  the  Muellcr-Murphy  (M-M)  and  Sharpe  explosion  models 
with  the  W2P  and  W3P  earthquake  models,  (a)  Displacement  amplitude  spectra 
calculated  for  two  different  low  frequency  values.  The  uppermost  explosion  curves 
were  calculated  for  an  explosion  yield  of  100  kt.  The  lower  explosion  curves  were 
calculated  for  a  yield  of  1  kt.  The  earthquake  models  have  the  same  low  frequency 
levels  as  the  Sharpe  curves,  (b)  1 .0/4.0  spectral  ratio  versus  moment  for  the 
models.  In  both  figures  the  explo.sion  models  were  calculated  u.sing  equations  (1.1) 
and  (1.5- 1.9)  with  a  P  wave  velocity  of  4.6  km/s.  For  the  Sharpe  model  Po  =  200  bars 
and  k=9  and  for  the  M-M  model  equations  (1.6-1 .9)  were  used.  The  earthquake 
models  were  calculated  u.sing  equations  (1.10-1.12)  with  a  P  wave  velocity  of  6  km/s 
and  a  stress  drop  of  100  bars.  For  (b)  a  .single  relation  appropriate  for  the  earthquake 
.source  parameters  was  u.sed  to  calculate  the  moment  from  the  low  frequency  level. 
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have  assumed  that  the  rupture  velocity  is  equal  to  the  shear  wave  velocity.  This  gives 
the  maximum  high  frequency  values  for  the  W3P  model  since  at  higher  rupture  veloci¬ 
ties  the  spectrum  decays  as  and  at  lower  velocities  both  the  corner  frequency  and 
spectral  amplitudes  above  the  comer  frequency  decrease  (Evemden  et  al.,  1986,  Figure 
16).  (When  the  W3P  model  has  rupture  velocities  slower  than  the  shear  wave  velocity 
the  asymptotic  form  of  the  spectral  shape  is  no  longer  as  good  an  approximation  since 
the  transition  between  the  flat  spectrum  at  low  frequencies  and  the  to“^  falloff  at  high 
frequencies  occurs  over  a  wider  frequency  band).  Figure  la  shows  the  large  difference 
in  high  frequency  (>  10  Hz)  energy  content  between  the  W3P  and  explosion  models. 
Differences  between  the  W2P  model  and  the  explosion  models  are  smaller  ana  pri¬ 
marily  due  to  slightly  higher  comer  frequencies  of  the  explosion  models. 

One  of  the  more  promising  regional  discriminants  based  on  studies  of  Western 
U.S  earthquakes  and  explosions  are  spectral  ratios  of  regional  phases  (Murphy  and 
Bennett,  1982;  Bennett  and  Murphy,  1986,  Taylor  et  al.  1988).  The  differences  in  the 
high  frequency  falloff  between  the  four  models  discussed  above  predict  very  different 
source  P  wave  spectral  ratios.  Figure  lb  shows  the  log  of  the  spectral  ratio  of  1.0  Hz 
to  4.0  Hz  energy  as  a  function  of  moment,  using  the  same  parameters  as  Figure  la. 
For  small  moments  the  ratios  are  "nity  since  both  frequencies  are  on  the  constant  part 
of  the  spectrum.  For  large  moments  the  ratio  is  a  constant  value  greater  than  unity, 
since  both  frequencies  are  on  the  decaying  pan  of  the  source  spectra.  For  these  higher 
moments  the  W3P  1.0  to  4.0  Hz  ratio  is  larger  than  for  the  other  models  due  to  the 
co"^  decay  rate,  and  suggest  the  source  spectral  ratio  would  be  a  reliable  discriminant 
between  eanhquakes  and  explosions,  assuming  the  W3P  and  either  explosion  model  is 
valid.  For  moments  from  about  10^^  to  10^^  dyne-cm,  the  frequencies  fall  on  either 
side  of  the  source  comer  frequency,  and  the  behavior  of  the  spectral  ratio  curves  in 
Figure  lb  reflect  the  different  comer  frequency  scaling  of  each  model.  If  the  W2P 
model  for  eanhquakes  is  valid,  the  corner  frequency  scaling  of  both  the  W2P  and 
either  explosion  model  will  play  a  large  role  m  the  success  of  a  spectral  ratio  discrim¬ 
inant  based  purely  on  differences  in  source  spectra. 

From  equations  (1.10  -  1.12)  we  can  find  a  relation  for  both  the  W2P  and  W3P 
models  between  low  frequency  level,  stress  drop  and  comer  frequency  if  we  hold  the 
medium  dependent  parameters  and  ob.serving  distance  constant; 

D'^fw^O)  oc  (1.13) 

(f/)- 


Given  two  earthquakes  with  the  same  moment  (i.e.  same  low  frequency  level)  but 
different  stress  drops,  the  one  with  the  high  stress  drop  has  both  a  higher  comer  fre¬ 
quency  and  more  high  frequency  energy  than  the  event  with  the  lower  stress  drop.  A 
similar  relationship  can  be  found  for  explosions  from  equations  (1.2)  and  (1.5)  with  the 
pressure  replacing  the  stress  drop; 


f2'’(co— >0)  « 


P. 
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(1.14) 
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where  i  ==  0  for  the  Sharpe  model  and  i  =  2  for  the  Mueller  Murphy  model.  For  explo¬ 
sions  the  comer  frequency  is  inversely  proportional  so  the  eh.stic  radius  (Murphy. 
1977J- 


a 

2nR^, 


(1.15) 


Equations  (1.13)  and  (1.14)  show  that  (for  earthquakes  and  explosions  v.ith  the  same 
low  frequency  spectral  level)  the  ratio  of  explosion  corner  frequency  io  earthquake 
comer  frequency  depends  on  the  ratio  of  explosion  pressure  to  earthquake  stress  drop. 
For  example  the  Sharpe  model  with  =  200  bars  will  have  about  the  same  comer  fre¬ 
quency  as  W2P  model  high  stress  drop  earthquakes  (100  -  2(X)  bars).  In  this  case  a 
spectral  ratio  discriminant  will  not  work  very  well  (see  Figure  lb).  For  extremely  high 
stress  drop  events  (>  8(K)  bars),  the  corner  frequencies  of  the  W3P  and  Sharpe  model 
would  be  about  equal,  but  the  difference  in  falloff  slope  should  still  allow  a  spectral 
ratio  discriminant  to  work  if  one  of  the  frequencies  chosen  is  sufficiently  higher  than 
the  corner  frequency.  The  Mueller-Murphy  pressure  P2  given  in  equation  (1.9) 
depends  on  both  yield  and  depth  complicating  any  simple  comparisons.  Nevertheless,  it 
is  apparent  from  Figure  1  that  in  the  range  from  1  to  100  kt  (for  the  parameters 
assumed  in  Figure  1).  high  stress  drop  (100  -  200  bars)  W2P  earthquakes  will  not 
discriminate  very  w'ell  from  these  Mueller-Murphy  explosions  or,  the  basis  of  source 
spectral  ratios. 


North  American  Data 

In  a  previous  paper  (Walter  et  al.,  1988.  included  as  Appendix  A)  w'e  tested  some 
of  the  specific  predictions  of  the  V/2P,  W3P,  and  Sharpe  models  with  high  sample  rate 
digital  data  recorded  at  small  hypocentral  distances  (2  -  40  km).  The  earthquake  data 
consist  of  events  from  the  San  Jacinto  stfike-slip  fault  region  near  Anza,  California,  the 
subduction  zone  near  Oaxaca,  Mexico  and  the  Long  Valley  caldera  and  surrounding 
area  near  Mammoth  Lakes,  California.  The  explosions  are  NTS  explosions  from 
Pahute  Mesa.  To  minimize  the  effects  of  attenuation  we  compared  the  variation  of  the 
earthquake  high  frequency  amplitude  with  moment  for  each  region  using  a  constant  set 
of  stations  and  events  with  similar  depths  and  ray  paths.  The  determination  of  moment 
based  on  low  frequency  spectral  amplitudes  is  relatively  insensitive  to  attenuation  and 
thus  any  increase  in  earthquake  high  frequency  spectral  amplitude  with  moment,  within 
each  region,  can  be  attributed  to  a  source  effect.  When  comparing  the  absolute  spectral 
amplitude  of  explosions  with  earthquakes,  however,  we  need  to  take  into  account 
attenuation,  site,  and  path  effects. 

Figure  2a  displays  10  Hz  P  wave  spectral  amplitude  as  'i  function  of  moment  for 
a  variety  of  earthquakes  and  NTS  Pahute  Mesa  undcrgiound  nuclear  explosions. 
Theoretical  curves  have  been  drawn  for  the  Sharpe  explosion  mcxiel  and  lor  the  W2P 
and  W3P  earthquake  models.  The  sigmas  denote  the  strcr.s  drop  in  bars  tor  each 
earthquake  curve.  In  Figure  2a  the  earthquake  10  Hz  amplitude'  arc  consistent  with 


Anza,  Oaxaca  and  NTS  Events 


Log  Moment  (dyne-cm) 


Figure  2a.  P  wave  displacement  spectral  observations  of  earthquakes  (solid 
symbols)  and  NTS  explosions  plotted  as  a  function  of  moment,  and  compared 
with  theoretical  curves.  All  symbols  and  curves  are  normalized  to  a  hypocentral 
distance  of  10  km.  No  attenuation  correction  was  applied  to  the  data.  Events 
where  more  than  one  station  recording  was  used  have  reletively  larger  symbol 
size.  Theoretical  curves  are  for  constant  stress  drop  denoted  by  a  in  bars. 

The  W3P  model  is  given  for  constant  rupture  velocity  equal  to  the  shear  wave 
vel(x:ity.  The  Sharpe  curve  was  calculated  assuming  Po  =  130  bars  and  k=9.4. 
values  appropriate  for  NTS.  (a)  The  10  Hz  P  wave  displacement  spectral 
amplitude  as  a  function  of  moment,  (from  Walter  et  al.,  1988) 
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both  earthquake  models  and  the  explosions  show  about  an  order  of  magnitude  scatter 
near  the  theoretical  Sharpe  amplitude  prediction.  Figure  2b  compares  30  Hz  P  wave 
amplitudes  versus  moment  w'ith  the  same  theoretical  curves.  The  earthquakes  show  a 
continued  increase  in  30  Hz  amplitude  with  moment  consistent  with  an  fo  *'  falloff 
model  and  in  disagreement  with  an  0)“^  model  for  larger  moments. 

The  explosion  spectral  amplitudes  are  from  both  hardrock  and  borehole  sites 
(hypocentral  distances  15-25  km)  and  sediment  sites  (2-12  km).  One  explosion  was 
recorded  on  both  hardrock  and  sediment  sites;  the  10  and  30  Hz  spectral  amplitudes  at 
the  hardrock  sites  were  approximately  a  factor  of  two  less  than  at  the  closer  sediment 
sites.  On  the  basis  of  limited  velocity  and  attenuation  data,  Walter  et  al.  (1988) 
estimated  that  the  sediment  recordings  should  approximate  the  source  spectral  ampli¬ 
tude  at  10  Hz  and  underestimate  it  by  about  a  factor  of  3  at  30  Hz.  The  earthquake 
amplitudes  are  all  from  seismic  stations  located  at  hardrock  sites  (hypocentral  distances 
5  -  40  km).  Based  on  estimates  of  the  attenuation  at  Anza  by  Hough  et  al  (1988),  we 
expect  the  10  and  30  Hz  spectral  amplitudes  to  underestimate  the  source  spectral 
amplitudes  by  about  a  factor  of  2  to  3.  Site  effect  studies  at  Anza  show  an  increase  in 
10  and  30  Hz  spectral  amplitude  by  a  factor  of  about  5  above  the  free  surface  effect. 
Therefore  we  expect  the  Anza  spectral  amplitudes  shown  in  Figures  2a  and  2b  to 
slightly  overestimate  the  source  spectral  amplitudes. 

The  explosions  in  Figure  2b  separate  into  two  groups:  one  showing  two  orders  of 
magnitude  less  30  Hz  P  wave  energy  than  predicted  by  the  Sharpe  model,  the  other 
closer  to  the  Sharpe  prediction.  The  group  with  less  30  Hz  energy  was  detonated  at 
depths  between  600  and  700  m  while  the  group  with  more  30  Hz  energy  was 
detonated  at  depths  of  800  to  I  UK)  m.  Thus  the  shallower  group  has  a  high  frequency 
spectral  decay  between  10  and  30  Hz  that  is  much  greater  than  the  predicted  by 
the  Sharpe  model.  The  shallower  explosions  were  detonated  near  the  level  of  the  water 
table  (about  650  m)  at  Pahute  Mesa  and  we  suggested  that  the  location  of  the  explo¬ 
sion  relative  to  the  water  table  may  cau.se  the  difference  in  falloff  slope  between  the 
deep  and  shallow  events  (Walter  et  al.,  1988).  These  results  are  similar  to  the  change 
in  high  frequency  slope  observed  by  Denny  (1990)  and  used  by  Taylor  and  Denny 
(1990)  to  model  spectral  differences  observ'ed  between  NTS  and  KTS  explosions.  Tay¬ 
lor  and  Denny  (1990)  show  that  for  weak  porous  rock  the  radiated  shock  wave  divides 
into  a  two  wave  .system.  In  terms  of  the  Shaqie  model  this  is  equivalent  to  introducing 
a  rise  time  in  the  pressure  function  at  the  elastic  radius,  so  the  far-field  P  spectra  then 
decays  as  0)“^  at  high  frequencies. 

Figures  2a  and  2b  imply  that  Anza  earthquakes  have  more  30  Hz  P  wave  energy 
than  NTS  Pahute  Mesa  explosions  detonated  near  the  water  table,  in  contrast  to  the 
Evernden  et  al.  (1986)  proposal.  Chael  (1988)  has  noted  that  the  slope  of  the  Pg  spec¬ 
tra  discriminate  between  NTS  explosions  and  western  U.S.  earthquakes  in  the  opposite 
sense  of  the  Evernden  et  al.  (1986)  proposal  since  the  explosions  had  steeper  high  fre¬ 
quency  spectral  slopes  than  the  earthquakes  .Murphy  and  Bennett  (1982)  and  Taylor  et 
al.,  (1988)  used  low  to  high  frequency  spectral  ratios  for  P^,,  Pg,  and  Lg  and  noted  that 
the  NTS  explosions  had  less  high  frequency  energy  than  the  western  U.  S.  earth¬ 
quakes.  In  contrast,  Taylor  and  Marshall  (1990)  compared  KTS  events  to  shallow 
Asian  earthquakes  and  found  that  the  explosions  had  more  high  frequency  P  wave 
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Fig,  2b.  The  30  Hz  P  wave  displacement  spectral  amplitude  as  a  function 
of  moment  (from  Walter  et  al.,  1988) 
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energy  than  the  earthquakes.  The  materia!  pri)penies  and  depth  to  tiie  water  table  at 
KTS  are  very  different  from  NTS.  This  raises  the  questions  of  vvluit  are  the  comer 
frequency  scaling  and  high  frequency  decay  slofie  for  Soviet  KTS  espf.wious,  and, 
whether  Soviet  explosions  can  be  discriminated  from  earthquakes  in  cither  the  manner 
of  Evemden  et  al.  (1986),  or  the  manner  of  Chael  (19S8).  We  seek  a  preliminary 
answer  to  these  questions  by  examining  the  regional  seismic  data  recorded  m  the  vicin¬ 
ity  of  the  Soviet  East  Kazakh  Test  Site. 


Central  Asian  Data 

Soviet  JVE.  Figure  3  shows  the  location  of  the  September  14,  1988  Soviet  JVE 
explosion  (NEIC  01^=6.!,  about  650  m  depth)  and  the  four  near  regional  seismograph 
stations  which  recorded  the  event  and  w'hose  digital  seismograms  are  available  to  us. 
The  three  clo.sest  stations,  Karkaralinsk  (KKL,  A  =  255  km),  Bayanau!  (BAY,  A  =  255 
km)  and  Karasu  (KSU,  A=  160  km)  were  in.stallcd  as  part  of  a  ccroperative  data  collec¬ 
tion  effort  between  the  Natural  Resources  Defense  Council  (N'RDC)  in  the  L'.  S.  and 
the  Soviet  Academy  of  Sciences  (SAS)  in  the  C.  S.  S  R  (Priestley  et  al  ,  1990).  Both 
short  period  (I  s)  and  broadband  (15  s  free  period)  high  sample  rate  (.2(X)  sample/s) 
seismographs  recorded  the  Soviet  JVE  at  the  NRDC-SAS  sites.  The  seismometers 
were  located  in  the  same  vaults  as  the  previous  NRDC-SAS  cooperative  experiment 
from  1986  -  1987  (Berger  et  al.,  1987).  All  three  of  those  stations  are  located  on  simi¬ 
lar  granitic  intrusions  of  Paleozoic  to  early  Mesozoic  age  (Leith,  1987).  The  seismo¬ 
graph  located  at  Talgar  (TLG,  A  =  74{)  km)  is  a  digital  seismograph  operated  by  the 
SAS  and  has  a  free  period  of  1.6  seconds  (Priestley  et  al..  1990)  and  a  flat  velocity 
response  to  about  15  Hz.  TLG  is  situated  on  the  Northern  Hanks  of  the  Zaili-.Matau 
mountains  and  is  sited  within  Precambrian  and  lower  Paleozoic  crystalline  nvks. 

Seismograms  of  the  Soviet  JVE  .‘‘■om  the  vertical  short  perifKi  components  at  the 
three  NRDC-SAS  sites  are  shown  in  Figure  4.  Tlie  KKL  seismogram  shows  an  emer¬ 
gent  low  frequency  first  arrival  followed  about  three-cjuaners  of  a  second  later  by  a 
low  amplitude  higher  frequency  arrival.  The  BAY  .seismogram  shows  a  single  arrival 
containing  high  frequencies.  The  stations  KKL  and  BAY  are  past  the  P^  -  Pg  cross¬ 
over  distance  predicted  from  the  Eastern  Kazakhstan  vekK’ity  stiaicturc  (Antonova  et 
al.,  1978;  Leith  1987;  Priestley  et  al.,  1988),  and  we  assume  the  first  low  amplitude 
arrival  is  P^.  The  relatively  low  frequency  content  of  the  first  KKL  P,.  arrival  is  con¬ 
sistent  with  this  arrival  being  a  head  wave;  however,  the  relatively  high  frequency  con¬ 
tent  of  the  second  low  amplitude  KKL  arrival  and  the  BAY  P^,  arrival  are  more  con¬ 
sistent  with  their  being  turning  rays.  The  Pj.  wavetrain  is  dominated  by  the  large 
amplitude  dilatational  arrival  on  the  vertical  and  radial  ct)mponents  following  the  P„ 
arrival  by  about  1.5  .seconds.  The  dilatational  motion  of  this  arrival  suggests  a  super¬ 
critical  reflection  from  the  Moho.  .Synthetic  reflectivity  seismograms  computed  for  the 
DSS  velocity  structure  (Ixith,  1987)  agree  in  both  the  P„  -  P..|P  time  interval  and  with 
the  large  amplitude  and  relatively  simple  shape  of  P.,,P  nb>eiv<'ii  at  BA>  au  i.  particu¬ 
larly,  at  KKL  (Walter  and  Patton,  1990,  f  igure  I  ).  The  closest  recording  of  the  Soviet 
JVE  was  made  at  KSL  within  the  P„-P..  crossover  ihstance.  Both  the  short  and  long 
period  records  show  a  3  Hz  rcson;ince  associated  with  tlic  sue  (Bercer  et  al..  1988; 


Fig  3.  Map  showing  hxraiions  of  Soviet  events  and  stations  used  in  this  study. 
The  explosions  are  denoted  by  stars  and  the  earthquakes  are  denoted  by  squares. 
The  .September  14,  I98S  Soviet  JVE  explosion  is  the  easternmost  star.  The 
.sei.smic  recording  ,s(ati<ins  are  .shown  as  triangles  with  their  three  letter 
abreviations. 
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Priestley  et  al.,  1990).  The  filtered  KSU  trace  in  Figure  4  shows  a  large  amplitude 
arrival  following  the  first  arrival  by  about  two  seconds,  which  is  consistent  with  the 
time  for  the  arrival  at  this  distance  range. 

P  wave  spectra  of  the  Soviet  JVE,  both  corrected  and  uncorrected  for  attenuation, 
from  data  recorded  at  BAY  and  KKL  are  shown  in  Figure  5.  The  spectra  from  KSU 
are  contaminated  by  the  site  resonance  and  tire  not  shown.  In  order  to  estimate  the 
effect  of  any  site  resonance  at  KKL  we  computed  the  ratio  of  the  surface  to  borehole 
(66  m  deep)  spectral  amplitude  from  the  first  six  seconds  of  the  P  wave  from  the  m^,  = 
4.6  earthquake  of  May  26,  1987.  The  surface  to  borehole  spectral  ratio,  shown  in  Fig¬ 
ure  6  displays  little  evidence  of  any  near  surface  site  effect  between  1  and  10  Hz,  but 
there  is  an  amplification  of  the  surface  spectra  of  about  2  to  4  in  the  U)  -  30  Hz  range. 
Thus  we  expect  site  effects  at  KKL  to  have  little  or  no  effect  on  the  shape  of  the  spec¬ 
tra  less  then  10  Hz  and  a  small  effect  on  decay  rates  above  10  Hz.  The  first  arrival  of 
this  earthquake  was  not  recorded  at  BAY  but  a  comparison  of  the  surface  and  borehole 
(99  m  deep)  P  coda  in  the  7.2  to  6.0  km/s  group  velocity  window  shows  no  evidence 
of  surface  resonance  in  the  1  to  6  Hz  frequency  band  where  the  surface  recording  has 
a  reasonable  signal-to-noise  ratio.  Thus  we  do  not  expect  the  BAY  spectral  shape  to  be 
strongly  effected  by  site  effects  between  1  and  6  Hz.  The  borehole  instruments  at  the 
NRDC-SAS  sites  had  been  removed  by  the  time  of  the  Soviet  JVE. 

Spectra  for  an  equal  length  data  window  prior  to  the  P^  or  Pg  arrival  were  com¬ 
puted  to  indicate  the  noise  level  and  are  also  shown  in  Figure  5.  The  spectra 
uncorrected  for  attenuation,  should  give  a  minimum  estimate  for  the  source  corner  fre¬ 
quency  and  a  maximum  estimate  for  the  source  spectral  decay  slope.  The  uncorrected 
Pn  and  Pg  spectra  at  BAY  give  similar  comer  frequencies  of  2.6  -  4.0  Hz  and  decay  as 
0)'^  at  higher  frequencies.  The  uncorrected  P^  spectra  at  KKL  does  not  show  a  clear 
corner  frequency,  perhaps  because  of  the  complicated  multiple  arrival  as  discussed 
above.  The  uncorrected  P„  spectrum  at  KKL  has  an  apparent  comer  frequency  of  about 
1.5  Hz  and  decays  as  from  2  to  about  7  Hz  steepening  to  ^  from  7  to  25  Hz. 
The  corner  frequency  detemiined  by  the  intersection  of  the  a.symptote  with  the 
low  frequency  level  is  about  3  Hz.  To  correct  the  Pj,  spectra  for  attenuation  we  used 
the  results  of  Sereno  [1990]  who  determined  attenuation  to  be  given  by 
Q(0  =  3(X)f®-^  by  assuming  a  source  spectrum  similar  to  the  W2P  model  for  both 
earthquakes  and  mine  blasts,  and  fitting  spectra  from  data  recorded  at  the  three 
NRDC-SAS  sites.  Sereno  (1990)  notes  a  frequency  independent  value  of  Q  =  1175  fit 
the  data  equally  well.  To  correct  the  Pg  spectra  for  attenuation  we  have  used  the  Pg  Q 
of  20(X),  estimated  by  Given  et  al.  (1990)  to  be  a  minimum  value  for  frequencies 
above  10  Hz  on  the  basis  of  chemical  blasts  recorded  at  the  NRDC  stations  in  Kaza¬ 
khstan.  Given  et  al.  (1990)  note  that  the  differing  ray  paths  strongly  affect  the  Pg 
amplitude  so  that  Q>20(X)  should  be  regarded  as  a  tenuous  estimate.  The  attenuation 


corrected  BAY  P^,  Pg  and  KKL  Pg  spectra  for  the  Soviet  JVE  have  about  the  same 
comer  frequencies  as  the  uncorrected  spectra.  The  high  frequency  decay  slopes  of  the 
BAY  Pn  and  Pg  spectra  decay  approximately  as  in  the  range  of  5  -  25  Hz.  The 
KKL  Pg  spectra  decay  as  0)“'  in  the  range  from  1.5-7  Hz  and  as  ”  in  the  range 
from  7  -  25  Hz. 
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KKL 
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Fig.  4.  Short  period  vertical  seismograms  of  the  Soviet  JVE  recorded 
at  the  four  regional  seismograph  sites  shown  in  figure  3.  Two  expand¬ 
ed  plots  of  the  KSU  P  wave  are  shown,  the  upper  plot  is  the  original  data 
and  the  lower  trace  is  the  same  data  after  low  pass  filtering  with  a  four 
pole  butterworth  filter  with  a  comer  at  1 .5  Hz  in  order  to  remove  the  3 
Hz  site  resonance. 
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JVE  SPECTRA  AT  KKL 
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Fig.  5.  Spectra  of  Soviet  JVE  recorded  at  .stations  BAY  (top)  and  KKL  (bottom). 
Windows  u.sed  in  computing  spectra  arc  also  shown.  Dotted  lines  indicate  spectra 
of  an  equal  length  window  preceding  the  signal  to  indicate  the  noi.se  level.  Solid 
lines  are  the  uncorrected  signal  spectra  and  dashed  lines  have  been  corrected  for 
attenuation  as  described  in  the  text.  Corner  frequency  is  indicated  for  both  the  Pn 
and  Pg  pha.se.s  at  BAY  and  for  the  Pg  pha.se  at  KKL. 
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Fig.  6.  Ratio  of  surtace-lo-boreholc  P  wave  spectra  of  the  first  six 
seconds  of  the  May  26,  1987  event  (NEIC  mb  =  4.6,  740  km  distant) 
recorded  at  KKL.  The  borehole  seismometers  were  66  m  below  the 
surface.  The  specU'a  were  smoothed  once  with  a  five-point  running 
mean  prior  to  computing  the  ratio. 
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One  of  the  interesting  features  of  the  near  regional  recordings  of  the  Soviet  JVE 
is  the  large  high  frequency  SH  phases  apparent  on  the  transverse  components.  1  he  SH 
pulse  on  the  BAY  transverse  component  has  a  frequency  of  about  1  Hz.  The^c  phases 
have  been  modeled  as  tectonic  release  by  Walter  and  Patton  (1990).  It  is  often 
assumed  that  tectonic  release  has  a  negligible  contribution  at  frequencies  of  1  Hz  or 
greater,  but  these  seismograms  rai.se  questions  about  whether  such  assumptions  are 
valid  for  KTS  explosions. 

Comparison  of  Soviet  Earthquakes  and  Explosion.  In  order  to  isolate  differences 
in  source  spectra  from  propagation  effects,  we  attempt  to  minimize  the  effects  of 
attenuation  as  much  as  possible  by  comparing  the  spectra  of  two  Soviet  explosions 
with  the  spectra  of  two  earthquakes  recorded  over  similar,  but  nearly  reversed  paths. 
The  two  earthquakes  (whose  location  is  shown  in  Figure  3)  occurred  on  July  21,  1986 
and  May  26,  1987.  Both  events  were  =  4.6  (NEIC),  and  were  recorded  at  KKL. 
The  1986  event  was  about  610  km  distant  (NEIC  depth  33  km)  and  the  1987  event 
was  about  740  km  distant  (NEIC  depth  20  km).  The  two  explosions  are  denoted  by 
solid  stars  in  Figure  3.  The  easternmost  explosion  in  Figure  3  is  the  Soviet  JVE;  the 
other  explosion  is  an  m^,  =  4.9  (NEIC)  event  w'hich  occurred  at  Degelen  on  October 
18,  1988.  Both  were  recorded  at  the  Soviet  Academy  of  Sciences  station  Talgar 
(TLG)  near  Alma-Ata,  about  740  km  away. 

Figure  7  compares  instrument-corrected  spectra  from  each  event.  For  these 
near  reciprocal  propagation  paths  the  explosions  have  a  relatively  greater  high  fre¬ 
quency  decay  compared  with  the  earthquakes.  Fitting  the  log-log  spectrum  from  5  to 
20  Hz  with  a  least  squares  line  we  obtain  fall  off  slopes  of  about  and  °  from 
the  1986  and  1987  earthquake,  respectively.  In  contrast  the  JVE  explosion  has  a  larger 
decay,  about  between  5  and  15  Hz.  The  smaller  Degelen  explosion  shows  a  simi¬ 
lar  high  frequency  decay  where  the  signal  is  above  the  noise.  The  explosions  have 
high  apparent  comer  frequencies,  about  5  Hz.  The  earthquakes  do  not  show  a  well 
defined  apparent  corner  frequency  but  show  gradually  larger  decay  rates  with  increas¬ 
ing  frequency. 


Di.scussion 

Since  we  have  only  a  limited  Soviet  data  set,  any  conclusions  are  of  necessity, 
preliminary.  The  steeper  high  frequency  P  wave  spectral  decay  observed  for  the 
Soviet  explosions  compared  with  these  ob.served  for  Central  Asian  earthquakes  (Figure 
7)  with  nearly  reciprocal  paths,  is  similar  to  that  observed  for  Nonh  American  earth¬ 
quakes  and  explosions  (Chael,  1988;  Walter  et  al.,  1988).  However,  although  compar¬ 
ing  data  along  nearly  reciprocal  paths  minimizes  attenuation  effects  related  to  the  pro¬ 
pagation  path,  we  have  not  eliminated  effects  due  to  the  differences  in  source  depth 
nor  the  differences  in  .seismograph  site  effects.  'I'he  attenuation  corrected  near  regional 
JVE  spectra  show  falloff  slopes  between  and  0)“^.  If  the  attenuation  correction  is 
valid  and  Soviet  earthquakes  falloff  as  to^^,  then  the  significant  difference  in  slope 
observed  in  the  near  reciprocal  path  data  between  the  earthquakes  and  explosions  may 
be  partly  due  to  depth  dependent  attenuation  effects,  seismograph  site  effects  at  TLG, 
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Fig.  7.  Pn  spectra  ot  Soviet  explosions  and  earthquakes  recorded  on  nearly 
reciprocal  paths.  The  top  two  traces  are  the  September  14,  1988  Soviet  JVE 
explosion  and  the  October  18.  1988  Degelen  explosion  recorded  atTLG. 
The  bottom  two  traces  are  the  May  26.  1987  and  July  21,  1986  earthquakes 
recorded  at  KKL.  Earthquakes  depths  are  from  NEIC,  explosion  depths  are 
estimates  assuming  NTS  depth  scaling. 


or  some  a  combination  of  both. 

Recordings  of  the  Soviet  JVE  show  more  high  frequency  energy  in  the  2  to  5  Hz 
band  than  predicted  by  the  Sharpe  or  Mueller-Murphy  models.  Priestley  et  al,  (1990) 
used  the  Lg  amplitude  at  the  near  regional  NRDC-SAS  stations  to  estimate  a  yield  of 
about  120  kt  for  this  event.  Sykes  and  Ekstrom  (1989)  used  a  combined  mt,  -  Mj  mag¬ 
nitude  yield  relation  to  obtain  a  similar  value  113  Kt.  The  comer  frequency  for  the 
Soviet  JVE  (assuming  a  120  kt  yield)  predicted  by  the  Sharpe  model  from  equation 
(1.3)  and  (1.15)  is  about  1.4  Hz.  The  apparent  corner  frequency  predicted  by  the 
Mueller-Murphy  model  from  equations  (1.9)  and  (1.15)  for  a  depth  of  650  m  is  about 
0.95  Hz.  Because  of  the  peaking  near  the  corner  frequency  in  the  Mueller-Murphy 
model,  the  intersection  of  the  high  and  low  frequency  asymptotes  gives  a  slightly 
higher  value  for  the  corner  frequency,  about  1.7  Hz.  The  average  comer  frequency  of 
the  Soviet  JVE  (using  the  asymptotic  comer  frequency  for  the  KKL  Pg  spectrum,  and 
the  BAY  Pn  and  Pg  spectra)  is  about  3  Hz.  If  we  assume  that  this  comer  is  not  due  to 
spall,  site  effects  or  tectonic  release,  then  both  the  Sharpe  model  and  the  Mueller- 
Murphy  model  would  appear  to  under  predict  the  comer  frequency  from  Soviet  explo¬ 
sions.  However,  the  hu'ge  SH  phases  on  the  transverse  components  suggest  that  com¬ 
plications  to  the  explosion  source,  such  as  spall  and  tectonic  release,  may  be  adding 
significant  high  frequency  energy  to  this  event.  Walter  and  Patton  (1990)  used  surface 
wave  amplitudes  from  the  NRDC-SAS  stations  to  estimate  the  tectonic  release  moment 
at  about  one  fifth  to  one  tenth  that  of  the  explosion.  If  the  tectonic  release  source  spec¬ 
tra  resemble  the  typical  earthquake  spectra  (i.e  ii:.i  oia  lo  a  comer  and  then  decaying) 
it  would  be  difficult  to  explain  the  h‘gh  P  wave  corner  frequencies  on  this  basis. 
Recent  theoretical  calculations  lor  spall  in  a  velocity  model  appropriate  for  KTS 
(Barker  et  al,  1990,  McLaunhlin  et  al.,  1990)  give  a  peaked  spectrum  at  about  1  -  4  Hz 
for  the  JVE.  Synthetic  calculations  tui  P,,  ni.J  Pg  by  McLaughlin  et  al.,  (1990),  indi¬ 
cate  the  spectral  amplitude  at  1  -  4  Hz  may  be  dominated  by  spall.  Without  near 
source  acceleration  data  however,  it  is  difficult  to  detennine  whether  spall  is  a 
significant  contributor  of  1  -  4  Hz  energy  for  the  Soviet  JVE.  The  smaller  explosion 
shown  in  Figure  7  presumably  contains  tectonic  release  and  spall  signals  that  are 
different  from  the  JVE  event,  yet  they  have  similar  apparent  comer  frequencies  at  the 
TLG  station. 

The  relatively  high  frequency  content  of  the  near  regional  Soviet  JVE  spectra  can 
also  be  seen  by  comparing  1  to  4  Hz  spectral  ratios  with  theoretical  calculations.  We 
calculated  the  ratio  of  the  0.75  -  1.25  Hz  to  3.0  -  5.0  Hz  energy  for  the  attenuation 
corrected  spectra  in  Figure  5.  Tiie  BAY  P„,  Pg  and  KKL  Pg  regional  phases  give  spec¬ 
tral  ratios  of  2.2,  2.8,  and  5.7  respectively  The  Mueller-Murphy  and  Sharpe  models 
predict  larger  spectral  ratios  of  6.6  and  8.3  respectively.  For  comparison  the  W2P 
earthquake  model  gives  a  spectral  ratio  of  7.2  (100  bar  stress  drop)  or  12.3  (10  bar 
stress  drop),  and  the  W3P  earthquake  model  gives  a  spectral  ratio  of  about  64  (for 
either  1(X)  or  10  btu"  stress  drop),  when  constrained  to  have  about  the  same  low  fre¬ 
quency  level  as  the  Soviet  JVE.  Thus  the  Soviet  JVE  also  shows  both  smaller  explo¬ 
sion  spectral  ratios  than  predicted  by  theory,  and  than  expected  for  earthquakes. 
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We  also  calculated  the  same  spectral  ratios  of  the  spectra  of  the  m^,  =  6. 1  Soviet 
JVE  recorded  at  TLG  and  the  mi,  =  4.6  earthquake  that  occurred  near  the  TLG  station 
recorded  at  KKL  using  the  previously  assumed  P„  attenuation  correction  (Sereno, 
1990).  Again  the  explosion  has  a  small  spectral  ratio  of  about  0.25,  which  is  much 
less  than  the  ratio  of  about  4.2  calculated  for  the  earthquake,  even  though  the  earth¬ 
quake  is  about  1.5  m(,  units  smaller.  The  earthquake  spectral  ratio  value  of  4.2  is  con¬ 
sistent  with  a  10  bar  stress  drop  W2P  model,  but  the  explosion  value  is  again  much 
lower  than  either  explosion  model  predicts.  Because  the  paths  are  nearly  reciprocal 
the  ratio  of  the  two  different  spectral  ratios  does  not  depend  upon  the  specific  distance 
attenuation  correction  chosen,  but  it  may  be  influenced  by  depth  dependent  attenuation 
or  site  effp^'ts  at  TLG  as  noted  above.  This  may  be  reflected  in  the  less  than  one 
explosion  spectral  ratio  value,  indicating  that  there  is  some  peaking  of  the  TLG  spectra 
near  4  Hz.  Such  peaked  spectra  are  not  observed  at  the  closest  stations.  Overall  these 
observations  give  similar  results  to  Taylor  and  Marshall  (1990)  who  found  that  KTS 
explosions  gave  lower  0.5  -  1.0  Hz  to  2.0  -  3.0  Hz  spectral  ratios  when  compared  with 
shallow  Central  Asian  earthquakes  at  the  United  Kingdom  teleseismic  arrays. 


Conclusions 

The  North  American  earthquake  data  presented  in  this  paper  are  consistent  with  a 
constant  stress  drop  W2P  model,  and  have  more  high  frequency  energy  than  shallow 
NTS  explosions  at  Pahute  Mesa  which  show  a  spectral  decay  greater  than  (o“^  above 
10  Hz.  In  contrast  the  attenuation  corrected  near  regional  recordings  of  the  Soviet  JVE 
explosion  have  a  spectral  decay  between  co~^  and  from  10  to  25  Hz.  These  results 
are  approximately  consistent  with  the  hypothesis  of  Taylor  and  Denny  (1990),  who 
found  that  spectral  ratio  data  from  U.S.  explosions  could  be  fit  by  an  explosion  model 
that  decays  as  at  high  frequencies,  KTS  (Shagan  River)  explosion  spectral  ratios 
could  be  fit  with  a  model  that  decays  as  o)~^,  and  both  western  U.S.  and  Central  Asian 
earthquake  spectral  ratios  were  consistent  with  an  W2P  model.  In  addition,  we  found 
the  near  regional  recordings  of  the  Soviet  JVE  have  a  higher  comer  frequency  than 
predicted  by  both  the  Mueller-Murphy  model  for  granite  and  the  Sharpe  (using 
Evernden  et  al.  parameters  for  KTS)  model.  This  higher  corner  frequency  does  not 
appear  to  be  due  to  attenuation,  path  or  site  effects.  The  higher  corner  frequency 
causes  a  lower  observed  explosion  1  to  4  Hz  spectral  ratio  for  the  Soviet  JVE  than  is 
predicted  by  the  explosion  models.  Comparing  the  Soviet  JVE  explosion  with  an  earth¬ 
quake  recorded  on  nearly  reciprocal  paths  we  find  a  similarly  low  1  to  4  Hz  spectral 
ratio  for  the  explosion  and  a  higher  ratio  for  the  smaller  mj,  earthquake,  which  may  be 
partly  due  to  depth  dependent  attenuation,  and  site  effects  at  TLG.  Whether  the  higher 
corner  frequency  for  the  Soviet  JVE  is  due  to  a  need  for  a  different  explosion  scaling 
relationship  at  KTS  than  contained  in  the  models  reviewed  here,  or  is  due  to  complica¬ 
tions  in  the  simple  source  models,  such  as  spall  and  tectonic  release,  needs  to  be 
resolved  in  order  to  confidently  u.se  a  di.scriminant  based  on  differences  between  earth¬ 
quake  and  explosion  spectra. 
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Abstract 

The  regional  broadband  seismic  recordings  of  the  Soviet  Joint  Verification  Exper¬ 
iment  underground  nuclear  explosion  show  evidence  of  tectonic  release  in  the  form  of 
Love  waves  at  long  periods  and  high-frequency  S^S  phases  on  the  transverse  com¬ 
ponents.  We  present  the  results  of  some  simple  forward  modeling  of  this  data  in  order 
to  quantify  both  the  moment  and  mechanism  of  the  tectonic  release.  Using  10-20 
second  Love  wave  amplitude  and  phase  data,  we  estimate  the  tectonic  release  moments 
to  be  3.5x10'*’  N-m  for  a  north  striking  right-lateral  strike-slip  fault,  and  7.0xl0'^ 
N-m  for  a  northwest  striking  reverse  fault.  Using  the  results  of  these  end  member  tec¬ 
tonic  models,  the  Rayleigh  waves  can  be  satisfied  with  explosion  moments  of 
1.7xl0'*  N-m  (F  =  0.31)  and  2.4xl0'^  N-m  (F  =  0.44),  respectively.  To  resolve  the 
ambiguity  in  tectonic  mechanism  we  used  the  reflectivity  method  to  generate  synthetics 
and  compare  the  relative  amplitudes  of  the  SH  and  SV  pans  of  the  S^S  phase.  While 
the  results  favor  the  strike-slip  mechanism  of  tectonic  release,  application  of  source 
inversion  techniques  to  this  data  and  other  regional  data  will  be  needed  to  evaluate 
more  realistic  tectonic  release  models  involving,  for  example,  oblique  fault  slip  and 
dip. 


IntrtKluction 

On  .September  14,  198H,  the  Soviet  Union  conducted  an  underground  nuclear 
explosion  as  part  of  the  Joint  Verification  Experiment  (JVE).  Under  an  agreement 
between  the  Soviet  Academy  of  Sciences  and  the  Natural  Resources  Defen.se  Council, 
the  University  of  Nevada,  Reno  recorded  the  explosion  at  three  seismic  stations  sur¬ 
rounding  the  Soviet  test  site  at  near  regional  distances  (160-260  km).  A  description  of 
the  data,  instruments  and  recording  sites  is  given  in  Priestley  et  al.  (1990). 


Theoretically,  spherically  symmetric  explosions  should  generate  only  P  and  SV 
waves.  Observations  of  significant  SH  waves  from  explosions  have  been  explained  by 
two  mechanisms;  earthquake  faulting  triggered  by  the  explosion  (Brune  and  Pomeroy, 
1963),  and  release  of  tectonic  prestress  in  a  shattered  zone  surrounding  the  explosion 
(Press  and  Archambeau,  1962).  In  either  case  an  equivalent  double  couple  mechanism 
can  be  found  for  this  "tectonic  release".  Evidence  of  the  tectonic  release  from  the  JVE 
at  the  near  regional  stations  is  supported  by  both  the  long  period  Love  waves,  and  the 
short  period  phases  arriving  with  a  group  velocity  of  about  3.5  km/s  on  the  transverse 
components,  particularly  at  BAY.  We  have  identified  the  most  prominent  of  these 
phases  as  S^S,  the  Moho-reflected  shear  wave,  as  shown  in  Figure  9.  Murphy  and 
Archambeau  (1986)  observed  a  short-period  coherent  SH  phase  from  the  Rulison 
underground  nuclear  explosion  on  the  transverse  component  of  a  very  close  station,  23 
km  away,  but  such  short  period  evidence  of  tectonic  release  is  not  commonly 
observed. 

Many  investigators  have  noted  anomalous  surface  waves  from  Soviet  underground 
explosions  in  Eastern  Kazakhstan  (e.g.  Rygg,  1979;  Helle  and  Rygg,  1984;  Given  and 
Mellman,  1985).  Tectonic  release  from  a  reverse  fault  striking  northwest  is  the  favored 
explanation,  since  some  large  F-value  (Toksoz  et  al.,  1971)  events  have  reversed  Ray¬ 
leigh  wave  polarities  at  all  stations.  The  45°  dipping  reverse  fault  is  the  only  mechan¬ 
ism  that  can  account  for  the.se  polarity  anomalies  at  all  azimuths. 

In  this  study  we  try  some  simple  forward  modeling  to  determine  the  size  and 
orientation  of  the  tectonic  release.  We  assume  the  mechanism  is  either  pure  strike-slip 
or  pure  reverse  on  a  45°  dipping  plane.  First  we  use  long  period  Love  and  Rayleigh 
amplitude  and  phase  information  to  estimate  the  tectonic  and  explosive  moments. 
There  is  an  ambiguity  however,  since  a  reverse  fault  on  a  45°  dipping  plane  and  a 
strike-slip  rotated  45°  with  respect  to  the  reverse,  produce  the  same  four-lobed  Love 
wave  radiation  pattern.  At  low  F-values  Rayleigh  Waves  cannot  resolve  the  ambiguity 
due  to  the  shallow  source  depth  of  explosions  (Patton,  1988).  To  resolve  the  ambiguity 
we  use  the  reflectivity  method  to  generate  synthetics  and  compare  the  relative  ampli¬ 
tudes  of  the  SV  and  SH  parts  of  the  S^S  phase. 


Surface  Wave  Modeling 

We  used  only  two  of  the  three  stations,  BAY  and  KKL  (see  Figure  3,  in  Chapter 
1)  in  our  analysis  becau.se  of  long  period  instrument  problems  at  KSU.  The  azimuths 
of  the  two  stations  from  the  JVE  differ  by  about  38  degrees.  Each  of  these  stations 
was  about  260  km  away  from  the  JVE.  The  data  were  corrected  for  the  effects  of  the 
instrument  at  long  periods  using  the  stopwatch  measurements  of  the  free  periods  (about 
22  s  at  BAY,  15  s  at  KKL)  and  as.suming  0.7  of  critical  damping.  The  corrected  data 
were  then  rotated  to  vertical,  radial  and  transverse  components  using  the  apparent  back 
azimuth.  We  used  arrival  times  from  co-located  University  of  Wisconsin  instruments 
that  had  an  omega  receiver  clock  (H.  Meyer,  personal  communication)  and  the  origin 
time  given  by  NEIC. 
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Fig,  8.  Data  and  synihciits  at  260  km.  The  first  traces  are  the  vertical 
( V)  and  transverse  (T)  component  recordings  of  the  Soviet  JVE  at  KKL 
and  BAY.  The  middle  traces  are  the  same  components  calculated 
for  an  explosion  plus  a  strike-slip  model  of  tectonic  release  in  the  DSS 
vehK’ity  model.  The  last  traces  were  calculated  for  an  explosion  plus 
a  reverse  fault  model.  Both  data  and  synthetics  have  been  low  pass 
nitered  with  a  four  pole  Butterworth  filter  at  2.5  H/. 


Phase  velocities  and  eigenfunctions  were  computed  for  the  velocity  mcxicl  BAY- 
Japan  (Priestley  et  al.,  1988)  using  the  method  of  Takeuchi  and  Saito  (1972).  The 
phase  spectra  of  the  Love  and  Rayleigh  waves  were  found  by  taking  the  rouriei 
transform  of  a  sixty  second  window  starting  at  a  group  velocity  of  about  3.8  km/s. 
These  phase  spectra  were  then  corrected  for  the  travel  path  using  the  calculated  phase 
velocities.  The  resulting  source  phase  delay  spectra  are  shown  in  Figure  10.  The 
theoretical  values  for  an  explosive  source  of  Rayleigh  waves  with  step  function  time 
dependency  is  3  71/ 4  radians  (0.375  cycles).  Note  the  rea.sonable  match  of  the  data  to 
the  theoretical  for  the  Rayleigh  waves,  particularly  at  BAY.  This  gives  us  confidence 
in  our  velocity  model  and  our  subsequent  interpretation  of  the  Love  wave  source  phase 
delay.  The  theoretical  values  for  a  double  couple  source  with  step  time  function 
dependency  is  7t/4  or  5  7C/4  radians  (0.125  or  0.625  cycles)  for  Love  waves,  depend¬ 
ing  upon  the  polarity  of  the  lobe  of  the  radiation  pattern.  The  Love  wave  data  in  Fig¬ 
ure  10  show  values  of  0.675  cycles  for  both  KKL  and  BAY  implying  they  are  on  the 
same  lobe  of  the  radiation  pattern.  Using  the  additional  information  that  the  long 
period  Love  wave  amplitudes  are  nearly  equal  constrains  the  orientation  of  a  suake-slip 
mechanism  to  between  about  355°  and  15°  in  strike,  or.  for  a  reverse  mechanism,  to 
about  310°  to  330°.  In  the  following  analysis  we  assume  a  strike-slip  mechanism 
striking  due  north  and  a  reverse  mechanism  striking  due  northwest. 

The  eigenfunction  were  combined  to  compute  synthetic  spectra  of  the  fundamen¬ 
tal  mode  surface  .  es  at  the  appropriate  distances  and  azimuths  of  the  two  stations. 
Assuming  th'=’  ;  onic  relea.se  is  co-located  and  simultaneous  with  the  explosion  we 
adjusted  the  moment  for  an  assumed  mechanism  until  a  reasonable  fit  was  obtained  in 
the  periou  range  of  5-20  seconds.  Figure  11  shows  the  data  and  synthetics  for  Love 
wavev  at  BAY  assuming  a  north  striking  strike-slip  mechanism  and  a  seismic  moment 
of  3.5  X  10 N-m.  A  similar  fit  is  obtained  at  KKL.  Using  an  explosion  moment  of 
l./x  10 N-m  with  the  strike-slip  tectonic  mechanism  fits  the  Rayleigh  wave  data  rea¬ 
sonably,  as  shown  in  Figure  1 1 .  As  noted  previously,  we  can  generate  the  same  syn¬ 
thetic  Love  wave  spectra  using  a  reverse  mechanism  and  a  moment  twice  as  large  as 
for  the  strike-slip  case  (7.0x10*^  N-m).  Assuming  a  reverse  mechanism  leads  to  an 
explosive  moment  of  2.4x  lO’^  N-m.  These  results  correspond  to  F  =  0.31  for  the 
strike-slip  case  and  F  =  ().44  for  the  reverse  case. 


Body  Wave  Modeling 

To  attempt  to  resolve  the  ambiguity  between  the  two  mechanisms  of  tectonic 
release,  we  looked  at  the  SH  and  SV  parts  of  the  phase.  Although  the  SH  radia¬ 
tion  patterns  are  the  same  for  the  two  mechanisms,  the  SV  patterns  differ,  causing 
different  relative  SV  and  SH  S^S  amplitudes  on  the  vertical  and  transverse  com¬ 
ponents.  We  used  the  reflectivity  method  (Muller,  1985)  and  the  same  velocity  struc¬ 
ture  as  for  the  surface  'vave  modeling  along  with  the  deep  sei.smic  sounding  (DSS) 
model  used  as  a  starting  model  by  Priestley  et  al.  (1988).  The  DSS  model  has  no  low 
velocity  zone  and  a  stronger  Moho  contrast.  The  resulting  differences  in  the  synthetics 
between  the  two  velocity  mcxlels  are  shown  in  Figure  12.  Both  models  give  very 
similar  long  period  (>  10  s)  pha.se  velocities,  but  differ  at  shorter  pericxi.s.  The 
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Frequency  (Hz) 


Fig.  10.  Observed  and  .synthetic  (a)  Love,  and  (b)  Rayleigh  wave  spectra  at 
BAY.  The  solid  line  is  the  amplitude  spectra  of  60  s  of  the  surface  waves, 
and  the  dotted  line  shows  the  preceding  60s  for  comparison.  The  dashed 
line  is  the  synthetic  spectra  of  the  fundemcntal  mode  surface  wave. 
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assumption  of  a  constant  high  Q  vckx'ity  structure  resultcti  in  large  amplitudes  of  the 
higher  mode  surface  waves  which  are  not  observed  in  the  data.  Ob.servations  of  the 
dispersion  of  Rg  phases  in  New  [England  give  crustal  models  with  very  low  Q  values 
in  the  uppemiost  kilometer  (Saikia  et  al.,  1990).  The  eastern  region  of  the  U.S.  is  often 
thought  to  be  analogous  to  Eastern  Kazakhstan.  For  these  rea.sons  we  used  an  attenua¬ 
tion  .structure  of  0=80  in  the  top  1  km,  Q  =  4(K)  from  1  to  2  km  depth  and  Q  =  2000 
from  2  km  down  for  F  waves,  and  3/4  of  these  values  for  S  waves.  In  all  ca.ses  we 
assumed  a  source  depth  of  640  m  and  a  step  source  time  function. 

Figure  9  shows  a  comparison  of  the  data  and  reflectivity  synthetics  at  KKL  and 
BAY.  The  KKL  Pn  and  P,„P  relative  amplitude  and  character  are  matched  quite  well 
by  the  strike-slip  synthetics.  The  reverse  mechanism  predicts  a  slightly  larger  P^^P  than 
is  observed.  The  S^S  phase  for  the  reverse  mechanism  shows  far  too  much  amplitude 
on  the  vertical  component  relative  to  the  transverse.  The  strike-slip  mechanism  comes 
closer  to  the  correct  amplitude  ratio,  but  the  data  show  an  even  higher  frequency  on 
the  transverse  than  the  synthetics.  The  strike-slip  mechanism  is  also  a  better  match  of 
the  relative  vertical  to  transverse  S^,S  amplitude  ratio  at  BAY.  However  the  P^^P  and 
Sfj,S  character  at  BAY  is  not  well  matched,  since  the  data  show  a  less  coherent  Pn,P 
arrival  and  a  stronger  Sll  arrival  than  the  synthetics,  possibly  due  to  the  limita¬ 
tions  of  the  simple  source  assumptions. 


Discussion 

The  orientation  of  the  tectonic  release  is  important  in  both  the  problems  of  deter¬ 
mining  the  yield  of  an  underground  explosion  using  seismic  data,  and  discriminating  it 
from  naturally  occurring  earthquakes.  Both  and  m^  measurements  can  be  affected 
signiHcantly  due  to  tectonic  release.  For  example  the  reverse  mechanism  usually 
assumed  in  Eastern  Kazakhstan  lowers  the  value  of  M,  compared  with  that  of  the 
explosion  by  itself.  Sykes  and  Ekstrom  (1989)  correct  their  estimate  by  adding  a 
lactor  of  0.15  due  to  tectonic  release  on  an  assumed  reverse  fault  for  the  Soviet  JVE 
event.  This  correction  of  changes  their  yield  based  purely  on  from  101  kilotons 
(kt)  to  145  kt.  Our  observations  however,  imply  that  the  strike-slip  model  is  more 
appropriate  for  this  particular  explosion,  so  that  the  correction  used  by  Sykes  and 
likstrom  (1989)  may  overestimate  the  M.,  yield  of  the  Soviet  JVE. 

Although  previous  large  lectonic  release  events  show  a  reverse  mechanism,  there 
is  some  evidence  for  a  strike-slip  stress  regime  in  the  East  Kazakhstan  region.  Pooley 
et  al.  (19X3)  analyzed  a  suspected  earthquake  that  occurred  in  1976  less  than  1(X)  km 
from  the  JV'E  site  and  found  it  fit  an  oblique-slip  on  a  northwest  striking,  steeply  dip¬ 
ping  plane.  They  also  noted  the  nearby  presence  of  the  Chingiz  fault,  a  right-lateral 
strike  slip  fault  striking  \4()'  W. 

The  results  of  this  study  are  based  on  some  simple  forward  mtxleling  of  the 
regional  broadband  seismic  recordings  at  two  stations.  Two  end  member  models  of 
tectonic  release  were  iiUHleled:  vertical  strike-slip  and  reverse  faulting  on  a  45°  dipping 
fd.ine  We  used  the  synthetic  seismognims  to  dctemiine  the  relative  excitation  of  the 
SV  arnl  Sll  parts  of  the  S,|,S  phase  for  each  of  these  mrxJels.  Although  the  chttracter  of 


the  Sn,S  phases  are  not  eonipletely  matched,  the  strike-slip  mtxlel  of  tectonic  release 
comes  much  closer  to  matching  the  relative  amplitudes  observed  in  the  data  than  the 
reverse  model.  ITte  modeling  included  the  simplifying  assumptions  of  step  function 
source  time  dependence  and  the  co-location  and  simultaneity  of  the  explosion  and  tec¬ 
tonic  events.  Additional  complexities  such  as  spall  were  not  considered.  The  use  of 
other  regional  data  and  more  tectonic  release  models,  such  as  oblique  fault  slip  and 
dip,  will  be  needed  to  further  evaluate  the  tectonic  relea.se  from  the  Soviet  JVE.  A 
source  inversion  of  this  near  regional  data  will  be  the  subject  of  future  work 
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Abstract 

A  model  for  the  average  far-field  P  and  S  wave  amplitude  spectra  from  a  circular 
crack  failing  in  tension  is  presented.  The  model  explicitly  avoids  specifying  the  details 
of  the  rupture  process,  instead  the  spectral  amplitude  is  determined  by  applying  physi¬ 
cally  realistic  constraints  to  an  idealized  spectral  shape  for  each  wave  type.  The  ideal¬ 
ized  spectral  form  is  given  by:  J  where  is  a  parameter  represent¬ 

ing  the  corner  frequency  and  the  parameter  determines  the  overall  scaling  of  the 
spectra.  The  P  and  S  spectral  models  are  constrained  by  fixing  the  low  frequency 
level  of  each  to  be  equivalent  to  a  point  source.  They  are  further  constrained  by 
equating  the  total  radiated  energy  with  the  available  elastic  energy.  A  final  constraint 
between  the  comer  frequencies  of  the  P  and  S  spectra  is  needed  to  completely  deter¬ 
mine  the  four  free  parameters.  For  the  type  of  model  considered  here:  an  equidimen- 
sional  fault  with  greatest  displacement  in  the  center  we  expect  the  P  to  S  corner  fre¬ 
quency  ratio  to  range  between  1  and  1.73  depending  on  the  details  of  the  rupture. 
This  model  gives  an  average  S/P  amplitude  spectrum  ratio  that  is  about  2.1  at  very  low 
frequencies  and  between  2.1  and  0.7  at  very  high  frequencies,  depending  on  the  P  to  S 
comer  frequency  ratio.  Applying  the  same  criteria  to  a  circular  shear  crack  we  show 
that  the  average  S/P  amplitude  is  about  7.1  at  very  low  frequencies  and  between  7.1 
and  2.4  at  very  high  frequencies,  again  depending  on  the  S  to  P  comer  frequency  ratio. 
The  tensional  crack  thus  has  a  lower  average  S/P  ratio  than  the  shear  crack,  and  such 
low  S/P  ratios  may  be  an  identifying  characteristic  of  tensile  or  tensile-equivalent  rup¬ 
ture.  Using  the  shear  and  tensional  models  we  construct  composite  spectra  that  have 
significantly  smaller  S/P  ratios  than  the  shear  cracks  alone  although  the  tensional 
events  radiate  less  than  10  per  percent  of  the  total  radiated  energy.  The  tension  crack 
average  spectral  model  presented  here  may  be  useful  as  a  first  approximation  for 
modeling  seismic  .sources  with  a  volumetric  component  in  them,  such  as  magma  injec¬ 
tion,  or  slip  on  a  non-planar  fault. 
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Conceptually  most  of  the  sources  of  naturally  occurring  seismic  radiation  in  the 
earth  have  been  modeled  as  shear  slip  on  a  perfectly  planar  surface.  However  observa¬ 
tions  of  exposed  faults  show  they  are  often  not  continuous,  not  straight  and  not  smooth 
(e.g.  Aviles  and  Scholtz,  1986;  Segall  and  Pollard,  1980).  Shear  slip  occurring  on  such 
non-planar  surfaces  might  be  expected  to  contain  a  significant  amount  of  motion  nor¬ 
mal  to  the  shear  slip  direction.  For  example  a  fault  plane  containing  a  bend  should 
tend  to  form  a  new  fault  at  the  bend  in  order  to  avoid  a  stress  singularity  (Andrews, 
1989).  During  shear  slip  this  triple  junction  must  undergo  a  volume  change,  which 
might  either  be  an  opening  or  some  complex  local  deformation  in  the  vicinity  of  the 
junction.  Other  sources  of  seismic  radiation  might  be  more  explicitly  related  to 
volumetric  changes  such  as  magma  injection  or  the  loss  of  compressional  strength  in  a 
region.  As  a  first  approximation  to  these  volumetric  sources  of  seismic  radiation  we 
develop  a  model  of  the  average  spectra  radiated  by  a  circular  crack  failing  in  tension. 

While  shear  crack  models  of  seismic  radiation  by  earthquakes  have  proven  very 
successful,  some  observations,  particularly  small  S  to  P  ratios  at  high  frequencies 
remain  inadequately  explained  by  such  models.  Recent  observations  of  anomalous  S 
to  P  ratios  near  Anza,  California  (Vernon  et  al.,  1991)  and  Guerrero,  Mexico  (Castro  et 
al..  1991)  seem  to  be  due,  at  least  in  part,  to  a  source  effect.  These  investigators  have 
speculated  that  motions  normal  to  the  fault  plane  may  be  causing  the  anomalously  low 
ratios.  The  idea  of  eanhquakes  containing  some  tensional  motions  goes  back  to 
Haskell  (1964),  who  suggested  that  in  order  to  explain  the  observed  ratio  of  short 
period  S  to  P  ratios  earthquakes  should  be  modeled  as  the  sum  of  subevents,  some 
shear  and  some  tensional.  This  idea  was  later  extended  by  Blandford  (1975)  to  explain 
other  observations  such  as  the  robustness  of  the  mt,:M5  earthquake/explosion  discrim¬ 
inant  and  high  P/S  corner  frequency  ratios. 

For  some  physical  processes  expected  to  generate  seismic  waves,  source  types 
other  than  shear  on  a  plane  are  expected.  For  example  magmatic  processes  have  been 
modeled  as  a  compensated  linear  vector  dipole  which  can  be  decomposed  into  a  tensile 
crack  plus  a  pure  volume  expansion.  Some  ob.servers  have  found  volumetric  com¬ 
ponents  in  moment  tensor  inversions  of  with  long  period  data  in  volcanic  volcanic 
regions  (e.g.  Julian,  1983;  Aki  1984  )  but  the  interpretation  of  these  results  as  non¬ 
shear  sources  is  controversial  (Wallace  1985).  As  pointed  out  by  Wallace  (1985)  it  is 
impossible  to  resolve  volumetric  change  in  a  source  using  wavelengths  longer  than  the 
source  depth.  Recently  Patton  and  Zandt,  (1991)  have  found  a  small  (17  per  cent)  but 
significant  volumetric  component  from  surface  wave  moment  tensor  inversions  of 
earthquakes  within  the  Basin  and  Range. 

A  number  of  investigators  have  shown  both  theoretically  (Odens  and  Martin, 
19X5,  Ixjmnitz-Adler  1991)  and  ob.servationally  (Tolstoi,  1967)  that  slip  on  a  rough 
surface  generates  motion  luxinal  to  the  plane  of  slip  leading  to  reduced  effective  nor¬ 
mal  stress  and  thus  reduced  friction.  This  reduction  of  effective  normal  stress  has  been 
proposed  as  playing  a  iintjor  role  in  stick-slip  sliding  expected  during  earthquakes 
(Brune,  John.son  and  .Slater,  1990;  Melosh,  1979).  These  nornial  motions  might  reason¬ 
ably  be  approximated  by  a  number  of  tensile  sub-events  added  to  the  overall  shear 
event.  These  tensile  sub-events  might  change  the  observed  spectra  from  that  expected 
for  a  pure  shear  source. 
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In  this  paper  we  will  first  detemjine  a  model  for  the  average  amplitude  spectra  of 
seismic  radiation  from  a  circular  tensional  failure.  In  contrast  with  other  tensional 
failure  models  (e.g.  Haskell,  1964,  Blandford,  1975)  we  do  not  attempt  to  explicitly 
define  the  rupture  location  and  velocity,  but  rather  assume  a  physically  realistic  spec¬ 
tral  shape  and  use  basic  energy,  point  .source  and  comer  frequency  relations  to  con¬ 
strain  the  models  free  parameters.  Because  the  rupture  and  slip  time  histories  remain 
unspecified  we  cannot  make  any  conclusions  regarding  the  change  in  spectral  shape 
with  azimuth  due  to  a  directivity  effect  (e.g.  Boatwright,  1980).  The  spectral  model 
used  here  represents  an  average  over  azimuth.  In  fact  at  relatively  high  frequencies  it 
has  been  proposed  that  the  rupture  velocity  undergoes  abrupt  changes  due  to  stress 
(e.g.  Madariaga,  1983)  or  frictional  (e.g.  Segall  and  Pollard.  1980)  heterogeneities. 
Therefore  simple  constant  mpture  models  of  directivity  may  not  be  applicable  at  high 
frequencies,  and  the  average  amplitude  used  here  may  be  more  appropriate. 

The  approach  used  here  is  similar  to,  but  more  general  than  that  used  by  Brune 
(1970)  and  Randall  (1973)  for  a  circular  shear  crack.  Brune  used  an  energy  balance 
criteria  but  detemiined  only  the  S  wave  spectrum.  Randall  (1973)  assumed  1(X)  per 
cent  conversion  of  elastic  energy  to  radiated  energy  and  fixed  the  P  to  S  comer  fre¬ 
quency  at  the  ratio  of  the  wave  velocities.  For  our  tensional  crack  model  we  allow 
both  to  vary.  We  then  apply  the  same  analysis  to  a  circular  shear  crack  and  compare 
the  spectra  from  the  two  types  of  events.  Finally  we  construct  some  composite  models 
containing  both  shear  and  tensile  faulting. 


Model 

We  shall  assume  an  idealized  far-field  displacement  spectral  amplitude  model  for 
the  average  seismic  radiation  from  a  tensile  crack  of  the  form: 
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where  Q^,  co/,  and  (oj*  are  parameters  to  be  determined  by  physical  constraint 
equations.  The  phase  spectrum  remains  unspecified  so  that  in  the  time  domain  a 
variety  of  pulse  shapes  will  have  this  amplitude  spectrum.  While  the  choice  of  mcxiel 
is  somewhat  arbitrary,  this  model  is  about  the  simplest  fomi  of  a  displacement  spec¬ 
trum  that  has  the  correct  asymptotic  behavior.  It  approaches  a  constant  value  (neglect¬ 
ing  near-field  terms)  as  co^O  that  is  proportional  to  the  moment,  and  as  co— >0°  it  falls 
off  as  co^^  which  is  steeper  than  the  co~*^  needed  to  keep  the  radiated  energy  finite. 
This  model  has  been  applied  successfully  to  many  earthquakes  being  mcxleled  as  shear 
cracks.  We  have  no  apriori  reason  to  use  a  different  spectral  mcxiel  in  the  case  of  a 
tensile  crack  failure  and  therefore  will  assume  this  model  in  the  following  analysis.  In 
an  appendix  we  show  that  the  same  analysis  can  be  applied  to  spectral  mcxlels  similar 
to  those  described  by  equation  (3.1)  but  with  a  different  high  frequency  falloff  rate. 
The  adoption  of  a  different  high  frequency  falloff  rate  slightly  changes  the  comer  fre¬ 
quency  scaling  constants  but  does  not  otherwise  affect  our  analysis. 
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The  geometry'  for  the  tensional  failure  considered  here  is  that  of  a  circular  crack 
in  the  interior  of  an  infinite  elastic  medium.  The  crack  occupies  a  circle  of  radius  a  in 
the  X]-X2  plane  as  shown  in  Figure  13.  The  crack  is  formed  by  axially  symmetric 
tension  along  the  Xj-axis  (o,).  which  is  perpendicular  to  the  plane  of  the  crack.  The 
displacements  and  elastic  energy  change  from  a  circular  or  "penny  shaped"  crack  fail¬ 
ing  in  tension  was  first  solved  by  Sack  (1946)  in  oblate  spheroidal  coordinates.  Sned¬ 
don  (1951)  gives  a  solution  in  the  more  tractable  cylindrical  coordinates.  The  net 
opening  in  the  Xvdirection  across  the  crack  is  given  by  (Sneddon,  1951); 
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where  we  have  assumed  the  Lame  constants  X  and  p  are  equal  (Poisson’s  ratio  is  1/4). 
From  equation  (3.2)  it  is  apparent  that  the  maximum  displacement,  which  we  denote  as 
e,  occurs  at  the  origin.  The  average  displacement  can  also  be  obtained  from  equation 
(3.2)  and  is  2/3  c.  We  note  that  the  displacements  given  by  equation  (3.2)  are  indepen¬ 
dent  of  the  stress  o,  and  o^.  W’e  also  note  that  the  di.splacements  of  the  circular  crack 
in  a  uniaxial  tensii'inal  field  are  equivalent  to  those  generated  by  a  circular  crack 
formed  by  a  uniform  internal  pressure. 

The  first  constraint  equation  is  provided  by  the  elastic  energy  available  to  be  radi¬ 
ated.  The  energy  difference  between  the  intact  medium  and  one  containing  a  circular 
tension  crack  is  (Sneddon.  PJ.'Sl): 
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(3.3) 


Equation  (3.3)  represents  the  total  energy  available  for  seismic  radiation,  the  rupture  of 
atomic  bonds  and  any  other  energy  consumptive  processes  that  occur  during  the  tensile 
failure.  The  seismically  radiated  energy  is  given  bv: 
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vshere  (t  and  [)  are  the  P  and  S  wave  velocity,  p  represents  the  medium  density,  and 
R  is  the  distance  lonn  the  crack  to  where  the  seismic  waves  are  recorded.  Using  the 
s[K‘ctral  sha[ies  troni  cqii.inon  I'.li  the  first  constraint  equation  is  obtained  by  setting 
; .U  equal  lo  !.v4i  via  a  con^tani  f  which  represents  the  fraction  of  the  elastic  energy 
that  eoes  into  seismic  railiatiorr 
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Fig.  12.  The  coordinate  system  for  the  tensile  crack  model.  The  circular  crack  of 
radius  a  lies  in  the  x  l-x2  plane,  centered  at  the  origin.  In  side  view  the  deformation 
in  the  x.1  direction  is  elliptical  with  the  maximum  opening  of  E .  occurring  at  the 
center  of  the  crack.  We  shall  assume  a  »  6  .so  that  the  crack  is  approximately  planar. 
The  radiation  is  ob.served  at  a  point  at  distance  R  and  angle  Q  from  the  origin  and  a 
distance  R’  from  a  specfic  point  .specified  by  r  andO  on  the  plane  of  the  crack. 
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It’  the  circular  crack  were  opened  by  a  sudden  internal  pressure  of  magnitude  03 
instead  of  external  tension  the  displacement  in  the  xj  direction  is  still  described  by 
equation  (3.2).  The  work  done  by  this  internally  driven  crack  is  given  by  the  pressure 
times  the  area  times  the  average  defonnation: 

^  ^23 

^  ->  2n‘d-a^£  a^ 

W  =  (03)(7ta^)(^e)  = - ; -  =  2— - .  (3.6) 

3  3  n 

This  is  twice  the  energy  in  the  equation  (3.3).  'I'hus  only  half  of  this  energy  goes  into 
deforming  the  medium  and  the  rest,  the  same  amount  as  in  equation  (3.3)  is  available 
for  seismic  radiation. 

The  second  and  third  constraint  equations  come  from  assuming  that  the  source 
approaches  a  point  source  as  the  frequency  of  seismic  radiation  goes  to  zero.  We  shall 
limit  ourselves  to  considering  only  those  cases  where  a  >>  e  so  that  the  failure  can  be 
approximated  as  occurring  on  a  plane  at  X3  =  ().  We  will  also  assume  that  the  source  to 
receiver  distance  is  large  compared  with  the  radius  of  the  crack  (R  »  a)  and  that  we 
are  in  the  far-held  for  frequencies  of  interest.  The  moment  density  tensor  is  then 
given  by  (Aki  and  Richards,  1980;  equation  3.20): 

'  X  0  o' 

m  =  0  X  0  Au3(r,  4),  t).  (3.7) 

.0  0  X  +  2|iJ 

The  far-held,  whole-space  displacement  is  then  given  by  Aki  and  Richards  (1980) 
equations  3. 19  and  14.4: 
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where  mp^^  is  an  element  of  the  moment  density  ten.sor,  G  represents  the  Green’s  func¬ 
tion,  S  the  area  of  the  crack  and  the  asterisk  represents  a  convolution.  Each  y 
represents  a  directional  cosine  of  the  radiation  from  the  source,  y^  =  0R'/9x|.  Since  we 
assumed  that  R  »  a  the  directional  cosines  are  approximately  constant  and  indepen¬ 
dent  of  r  and  0.  and  can  be  taken  outside  of  the  integral.  Solving  for  the  magnitude  of 
the  separate  P  and  S  displacements  and  taking  the  Fourier  transform  of  the  results  we 
obtain: 
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for  the  far-field  P  wave  spectra,  and 


^  f[Aii,(r.(?,a})e'‘'*^''l^dS  (3.10) 
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for  the  far-field  S  wave  spectra,  where 


ALi3(r,  0,  (o)  =  J  Aii^Cr,  0,  t)  e'“‘dt.  (3.11) 

—  oo 

In  the  above  equations  9  is  defined  in  Figure  13.  As  the  O)— >0  the  far-field  spectra 
given  in  equations  (3.9)  and  (3.10)  approach  a  constant  value,  proportional  to  the 
integral  of  the  final  offset  of  Auj  over  the  area  of  the  tensile  crack  and  independent  of 
the  time  history.  Using  the  final  offset  for  a  circular  fault  given  by  equation  (3.2)  we 
can  equate  the  model  in  equation  (3.1)  with  the  far-field  spectra  given  by  equations 
(3.9)  and  0.10)  at  co  =  0  to  obtain; 
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iaP(co=())i  =  ( i-t-2cWe)-  -  (3.12) 
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and 


1qS(oj  =  0)I  =  (2 sin 9 cos 9)  ■  (3.13) 
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The  factors  of  l-t-2cos^9  and  2. sin 0 cos 6  appearing  in  the  above  equations  are  the 
point  source  radiation  patterns  of  a  tensile  crack.  ITie  radiation  pattern  amplitudes  are 
shown  in  Figure  14,  note  the  S-wave  nodes  in  the  plane  and  perpendicular  to  the  plane 
of  the  crack.  The  P  wave  radiation  pattern  does  not  have  any  radiation  pattern  nodes. 
In  order  to  u.se  equations  (3.12)  and  (3.13)  to  determine  the  energy  content  of  the 
seismic  radiation  a  mean  .square  value  of  the  radiation  pattern  must  be  determined. 
Averaging  over  the  whole  focal  sphere; 

^  ■)  wD  2  47 

<  (  1  -t-  2  cos’9  )"  >  =  — ^ — r  f(  1  +2  cos"9  )^  — —  sin  9  d9  =  (3.14) 
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and 
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<  (  2  cos9  sin9  )^>  =  — ^ — r  f(  2cos9  sin9)^  — —  sin  9d9  =  (3.15) 
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P  Wave 


S  wave 


Fig.  i  P  and  S  point  soiiao  radiation  patients  lor  the  tensile  crack.  The  crack 
plane  is  indicated  hy  the  circle. 
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Equations  (3.5),  (3.12)  and  (3.13)  represent  3  constraints  on  the  four  unknowns  in 
our  model,  given  by  equation  (3.1).  One  further  constraining  equation  is  needed  before 
we  can  solve  for  the  unknowns  in  the  model.  The  linal  constraining  equation  must 
relate  the  P  and  S  corner  frequencies  to  each  other.  The  physical  rea.son  for  corner  fre¬ 
quencies  can  be  seen  by  examining  the  fomi  of  equations  (3.9)  and  (3.10).  These  equa¬ 
tions  are  in  the  fomt  of  one  Fourier  transfomi  in  time  and  two  in  space.  Thus  comer 
frequencies  arise  because  of  the  finite  duration  of  faulting  t.s  well  as  the  finite  area  of 
faulting.  For  an  equidimensional  fault  with  mpture  velocity  on  the  order  of  the  shear 
wave  velocity,  the  corner  frequencies  associated  with  the  two  spatial  transforms  and 
the  one  time  transform  should  be  similar  and  consistent  with  our  spectral  model 
assumption  given  by  (3.1). 

As  discussed  by  Sato  and  Hirasawa  (1973),  a  seismic  source  which  undergoes 
greater  displacement  over  a  longer  time  in  the  center  of  the  crack  than  at  the  edges 
leads  to  a  larger  P  wave  corner  frequency  than  the  S  wave  comer  frequency.  Molnar  et 
al.  (1973)  showed  that  in  the  case  of  instantaneous  rupture  the  P  to  S  comer  frequency 
ratio  is  equal  to  the  P  to  S  wave  velocity  ratio.  We  can  represent  the  P  wave  comer 
frequency  as  some  constant,  g  times  the  S  wave  corner: 

0)^’’ =  g  co|  (3.16) 

where  we  expect  g  to  fall  in  the  range  from  l<g<a/(3  or  between  1  and  V3=  1.73. 
Using  (3.16)  as  our  last  constraint  equation  we  solve  (3.5)  for  the  S  wave  comer  fre¬ 
quency; 
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(3.17) 


The  parameter  f  represents  the  fraction  of  elastic  energy  which  goes  into  seismic  radia¬ 
tion  and  might  be  expected  to  vary  greatly.  In  constant  rupture  velocity  models  such 
as  Sato  and  Hirasawa  (1973)  this  fraction  increases  with  faster  rupture  velocities.  In 
general  the  energy  radiated  will  be  related  to  the  form  of  Au3(r,  0,  t)  which  might  be 
quite  complex  and  not  adequately  modeled  by  constant  rupture  velocity  models.  Errors 
in  our  choice  for  the  parameter  f  will  not  change  the  low  frequency  spectral  level  but 
change  the  comer  frequency  by  f’^  and  the  high  frequency  level  by  f^^.  This  effect  is 
illustrated  in  Figure  15.  In  this  paper  we  will  assume  that  the  energy  partition  between 
radiation  and  energy  consumptive  processes  such  as  rupture  of  atomic  bonds  is  approx¬ 
imately  equal  and  f  =  0.5  .  Let  us  examine  .some  specific  cases  for  the  comer  fre¬ 
quency  ratio  given  by  (3.16). 

Case  g  =  o/p.  If  a  rupture  occurs  rapidly  enough  so  that  the  fault  can  be 
treated  as  radiating  essentially  instantaneously,  then  the  P  and  S  comer  frequencies  are 
proportional  to  the  vcIcK'ity  divided  by  the  fault  dimension.  This  can  be  seen  directly 
by  substituting  a  delta  function  time  dependence  of  Au^  into  equations  (3.9)  and  (3.10) 
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Frequency  (Hz) 


Fig.  14.  Tcn.silc  P  wave  spcelra  for  several  values  of  the  parameter  f,  the  fraction 
of  elastic  energy  radiated.  All  spectra  were  calculated  for  a  1km  radius  crack, 
opening  1  cm. 
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Displacement  Amplitude  (cm-s) 


Frequency  (Hz) 


Fig.  15.  P  and  S  spectral  pairs  for  tensile  crack  model  with  three  different  values 
of  g,  the  P  to  .S  comer  frequency  ratio.  All  spectra  were  calculated  for  a  1  km 
radius  crack  with  I  cm  of  maximum  opening.  Upper  two  pairs  of  spectra  were 
multiplied  by  factors  of  10  and  100  for  display  purposes. 
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(cf.  Molnar  et  al.,  1973).  In  this  case  the  P  to  S  corner  frequency  ratio  will  be  equal  to 
the  P  to  S  velocity  ratio,  or  g  =  V3=I.73.  Using  equation  (3.17)  we  find  an  S  wave 
comer  frequency  of  (o^  =  1.59 p/a  and  a  P  wave  corner  frequency  of  =  1.59a /a 
The  P  and  S  wave  spectra  for  this  case  are  shown  in  Figure  16. 

Case  g  =  /  .  If  the  effective  duration  of  rupture  is  long  compared  to  the  radius 
over  the  seismic  wave  velocity,  then  the  observed  comer  frequency  might  be  deter¬ 
mined  completely  by  the  mpture  time.  This  is  equivalent  to  treating  the  crack  as  a 
point  source  (R'  =  R)  for  the  frequencies  of  interest,  so  that  the  comer  frequencies  are 
detennined  only  by  the  rupture  time.  In  this  case  the  corner  frequencies  are  equal  and 
g  =  1.0.  Using  (3.17)  and  (3.16)  the  S  and  P  comer  frequencies  are 

co^  =  (0^  =  2.05 p/a.  Spectra  for  this  case  are  shown  in  Figure  16 

Case  g  =  1.4.  We  would  now  lil.e  to  consider  an  intermediate  case  when  g  = 
1.4.  While  this  choice  is  somewhat  subjective  we  are  motivated  by  similar  values  for 
the  P  to  S  wave  ratio  obtained  by  a  variety  of  investigators  for  a  circular  shear  crack 
which  involves  solving  integrals  completely  analogous  to  those  appearing  in  equations 
(3.9)  and  (3.10)  (e.g.  Sato  and  Hirasawa,  1973,  Molnar  et  al.,  1973,  Madariaga,  1976). 
The  S  and  P  wave  corner  frequencies  are  (o^  =  l.SOp/a  and  -  2.51  p/a.  Using 
equations  (3.9)  and  (3.10)  the  far-field  P  and  S  wave  amplitude  spectra  is  given  by: 
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The  P  and  S  wave  spectra  are  shown  in  Figure  16. 

The  lithostatic  pressure  increases  with  depth  approximately  as  300  bars/km  in  the 
cru.st.  At  the  5-35  km  depths  where  most  crustal  earthquakes  occur  the  lithostatic  pres¬ 
sure  will  be  of  the  order  of  thousands  of  bars.  During  an  earthquake  one  might  expect 
some  regions  to  be  under  local  tensile  stress  during  the  event,  perhaps  due  to  shear  slip 
on  a  non-planar  fault  or  some  other  mechanism,  which  then  returns  to  compressive 
stress  after  the  event  is  over.  In  such  a  case  we  can  model  a  circular  region  subject  to 
short  duration  of  tensile  stress  as  an  opening  and  closing  crack.  Making  a  simplifying 
assumption  that  the  opening  and  closing  crack  displacement  pulses  are  equal  but  oppo¬ 
site  in  sign,  the  riu"  field  displacement  is  given  by: 


u(t)  =  u„(t)-u„(t  +  x) 


(3.20) 


where  u,,  is  the  far  -field  liisplacement  of  the  opening  crack  and  x  is  the  time  the  crack 
remains  open.  This  simplifying  assumption  can  lead  to  radiated  energy  that  is  about  20 
per  cent  higher  than  twice  the  energy  of  the  opening  crack  for  particular  values  of  x. 
However,  the  values  of  x  we  will  use  later  in  the  paper  give  an  energy  for  the  opening 


and  closing  case  that  is  very  nearly  twice  that  of  the  opening  case  alone,  justifying  this 
simple  assumption  for  our  modeling  purposes.  Making  use  of  the  shift  theorem  we 
find  our  model  magnitude  spectra  are  multiplied  by  a  factor  of  (2-2cos(ot)^''^^\  This 
has  the  effect  of  reducing  the  low  frequencies,  leading  to  a  displacement  spectrum 
peaked  near  the  comer  frequency  as  shown  in  Figure  17. 


Comparison  with  Shear  Crack 

The  above  analysis  can  Just  as  easily  be  done  for  a  circular  shear  crack.  Let  us 
consider  a  circular  shear  crack  of  radius  a'  in  the  X]-X2  plane  subject  to  a  uniform 
applied  shear  stress  of  magnitude  a  which  slips  in  the  xi  direction.  The  relative  dis¬ 
placement  across  the  faces  of  the  crack  is  given  by  (Eshelby,  1957): 

24  02  r-, - ^ 

Au^  = - ^Va'’-r^  (3.21) 
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again  assuming  Poisson’s  ratio  is  1/4.  Again  the  maximum  displacement  which  we 
denote  as  e',  occurs  at  the  origin.  The  energy  difference  between  an  intact  infinite 
elastic  medium  and  one  containing  a  circular  shear  crack  with  the  displacement 
specified  by  equation  (20)  is  (Eshelby,  1957): 
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(3.22) 


We  use  the  same  spectral  amplitude  model  as  given  in  equation  (2),  with  the  parame¬ 
ters  ,  Qg  ,  (UjT,  and  cOg  denoted  with  primes  to  distinguish  them  from  the  tensile 
case.  Applying  the  similar  low  frequency  constraints  and  assuming  a  fraction,  f  of  the 
elastic  energy  is  radiated  and  the  P  to  S  comer  frequency  ratio  is  given  by  g'  we  can 
solve  for  the  four  parameters; 
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Fig.  16  Spectra  ol  an  opening  and  closing  tensile  crack.  Calculated  for  a 
1  km  radius  crack  opening  1  cm  in  0.7  seconds  then  closing. 
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where  where  arc  the  P  and  S  point  source  radiation  patterns  for  a  shear 

crack  (cf.  Aki  and  Richards,  1981),  We  note  that  if  f  =  0.5  we  obtain  values  for 
close  to  the  Brune  (1970.  1971)  value  of  2.34 (3 /a'  For  example  f  =  0.5  and  g'  =  1.4 
yields  cOjf  =2.31  p/ a'. 

We  have  assumed  that  the  general  spectral  shape  for  pure  shear  and  pure  tensile 
cracks  are  similar.  If  these  models  are  correct  then  predicted  observational  differences 
between  the  two  types  of  events  will  be  mostly  due  to  radiation  pattern  differences.  In 
addition  to  differences  in  the  spatial  character,  the  radiation  patterns  should  also  lead 
to  large  differences  in  the  observed  S  to  P  spectral  amplitude  ratios.  At  very  low  fre¬ 
quencies  (oi<<Oi^)  and  very  high  frequencies  (to>>a)jf)  the  P  to  S  spectral  amplitude 
ratio  has  the  same  form  for  both  the  shear  and  tensile  crack  cases.  The  tensile  crack 
S/P  ratios  are: 


I  Q^((t)-^0)l  _  I (a)-^o°)l  _  ^ 

IqP(O)->0)I  ~  r|<^  p5’  lQP(a)-^c»)|  "  rP^  g2p3- 

where  Rq<i,  and  R^  tire  the  S  and  P  wave  radiation  patterns.  The  shear  crack  ratios 
are  of  similar  form  to  (3.26)  with  R^,  R^  and  g'  replacing  R0<j,  R^  and  g.  For  the 
tensile  case  <1  R|<t)l  >  =  2/3  and  <1  Reti)l  >  =  5/3  so  that  at  low  frequencies  the  S  to  P 
spectral  amplitude  ratio  is  about  2.1.  At  high  frequencies  the  ratio  varies  from  about 
2.1  to  about  0.69  as  g  varies  from  1.00  to  1.73.  The  tensile  crack  S/P  ratio  is  shown  in 
Figure  18a  for  three  specific  values  of  g.  The  average  magnitude  of  the  shear  crack 
radiation  patterns  is  <lRe<t)l>  =  0.60  and  <lR0|j>l>  =  0.44  (Boore  and  Boatwright, 
1984)  so  that  at  low  frequencies  the  S  to  P  spectral  amplitude  ratio  is  about  7.1  and 
the  high  frequency  S/P  ratio  varies  from  7.1  to  2.4  as  g  goes  from  1.00  to  1.73.  The 
shear  crack  S/P  spectral  ratio  is  .shown  in  Figure  18b  for  three  values  of  g'.  These 
kirge  differences  in  S/P  ratios  may  be  an  identifying  characteristic  of  tensile  events. 


Composite  Events 

Some  earthquakes  near  Anza,  California  (Vernon  et  al.,  1991)  and  Guerrero, 
Mexico  (Castro  et  al.,  1991)  have  been  observed  to  have  high  frequency  S/P  ratios  less 
than  two.  As  shown  in  Figure  18b  the  lowest  S/P  ratio  expected  from  a  pure  shear 
source  when  averaged  over  the  focal  sphere  has  a  value  of  about  2.4.  These  investiga¬ 
tors  show  that  attenuation  estimates  for  these  areas  are  insufficient  to  explain  the 
anomalously  low  values  by  themselves  and  favor  at  least  a  partial  source  explanation. 
We  use  the  above  developed  tensile  and  shear  models  to  estimate  how  much  tensile 
motion  relative  to  shear  slip  is  necessary  to  cause  anomalous  S/P  ratios. 

The  presence  of  S  wave  nodes  in  the  shear  crack  radiation  pattern  and  the 
absence  of  P  wave  niides  in  the  tensile  crack  radiation  pattern  means  that  even  small 
amounts  of  tensile  motion  in  a  composite  event  may  be  apparent  if  the  observing  sta¬ 
tion  is  at  a  particular  angle  from  the  fault  plane.  However  scattering  and  a  non-uniform 
rupture  veUxaty,  even  for  a  pure  shear  seismic  source  might  be  expected  to  obscure  a 
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Fig.  17.  .S/P  .spcclral  ratios,  (a)  Tensile  crack  model  S/P  ratio  for  three 
values  of  P/S  corner  frequency  ratio,  g  .  (b)  Shear  crack  S/P  ratio 
for  same  throe  values  of  P/S  comer  frequency  ratio,  g’. 


simple  point  source  radiation  pattern,  particularly  at  high  frequencies.  Liu  and  Helm- 
berger  (1985)  have  noted  that  the  high  frequency  P  wave  radiation  pattern  of  earth¬ 
quakes  is  more  homogeneous  than  that  observed  at  lower  frequencies.  Therefore,  we 
will  compare  the  magnitude  of  the  tensile  and  shear  models  when  averaged  over  the 
focal  sphere.  Taking  the  ratio  of  shear  crack  to  tensile  crack  P  wave  displacement 
spectra  at  high  frequencies  (co>>to^'"  and  )  we  find; 

I  (to— )oo)l  _  k'^  c'  P '*7) 

IqP(co-.oo)!  ~  7^  (0)*’)^ 

where  we  have  used  =  k' P / a'  and  co^f^kp/a.  We  define  a  significant  contribution 
of  tensile  P  energy  as  occurring  when  the  shear  to  tensile  P  wave  spectral  ratio  given 
in  equation  (27)  is  less  then  one.  Using  the  average  magnitudes  of  the  shear  and  ten¬ 
sile  point  source  P  wave  radiation  patterns,  a  relative  slip  criterion  is  obtained  from 
(27): 


£  >  0.265 


(3.28) 


For  example  assuming  both  the  shear  and  tensile  cracks  have  the  same  P  to  S  comer 
frequency  ratio  of  1.4  (g  =  g'  =  1.4)  the  relative  slip  criteria  becomes  e  >  0.44e'.  In  other 
words  if  the  shear  crack  has  a  maximum  slip  of  10  cm  the  tensile  crack  must  have  a 
maximum  opening  of  4.4  cm  in  order  to  have  the  same  P  wave  di.splacement  .spectra  at 
high  frequencies.  Of  course  the  tensile  crack  does  not  have  to  consist  of  one  large 
event  but  may  be  made  up  of  many  smaller  subevents  with  total  opening  and  closing 
displacements  equal  to  4.4  cm.  If  the  smaller  events  add  incoherently  however,  more 
total  displacement  may  be  necessary.  In  the  case  of  opening  and  closing  events  the 
factor  of  (2-2cos(0)t))'^^  multiplies  the  tensile  crack  opening  in  equation  (3.28).  This 
factor  varies  from  0  to  2  but  at  very  high  frequencies  has  an  average  value  of  v'2.  The 
relative  slip  criteria  in  equation  (3.28)  does  not  put  any  constraints  on  the  radius  of  the 
tensile  subevents.  Equations  (3.3)  and  (3.22)  show  the  total  energy  of  the  cracks  is 
proportional  to  the  maximum  displacement  squared  times  the  radius.  Thus,  although 
the  total  opening  and  closing  displacements  may  be  large,  the  radius  and  total  energy 
of  the  subevents  may  be  small  comp;ired  with  the  shear  radius  and  radiated  energy. 

As  specific  examples  we  construct  composite  events  consisting  of  a  2  km  radius 
shear  crack  which  has  a  maximum  of  10  cm  of  slip  (corresponding  to  about  a  30  bar 
stress  drop)  and  several  opening  and  closing  tensile  subevents  which  are  assumed  to  be 
randomly  distributed  over  the  fault  plane  and  thus  over  the  duration  of  the  shear  event. 
As  shown  in  Figure  17,  the  subevents  are  peaked  near  their  comer  frequency  and  are 
expected  to  have  a  sigiiiiicant  contribution  only  at  high  frequencies.  Therefore  we  use 
the  approximation  that  the  composite  spectra  arc  the  square  root  of  the  sum  of  the 
scjuares  of  all  the  events  (e  g.  Joyner  and  Boore,  1986).  Tliis  should  be  a  conservative 
approximation  of  the  contribution  of  the  tensile  subevents. 
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F'iii.  IX.  (a)  Compivsiic  spccira  compared  with  pure  shear  crack.  Composite 
spectra  is  compo.sed  ot  a  2  km  radius  shear  crack  with  a  maximum  of  10  cm 
of  slip  plus  10  inchoherenlly  adding  opening  and  closing  tensile  cracks 
(200  m  radius,  maximum  opening  2  cm),  (b)  .S/P  ratio  of  .spectra  in  (a). 
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The  first  composite  event  consists  of  1  shear  crack  and  10  opening  and  closing 
subevents.  The  tensile  subevents  have  a  radius  of  200  m  and  a  maximum  opening  of  2 
cm  with  an  opening  closing  time,  t,  of  0.15  s  (t  -2  a/0.8  p).  Assuming  the  tensile  sub 
events  do  not  overlap,  they  cover  10  per  cent  of  the  area  of  the  shear  crack.  The  total 
energy  radiated  by  the  tensile  subevents  is  also  about  10  percent  of  that  radiated  by  the 
shear  crack.  The  composite  model  spectra  are  shown  in  Figure  19a  compared  with  the 
shear  crack  alone.  The  effect  of  the  tensile  subevents  is  visible  at  frequencies  above  2 
Hz.  Figure  19b  shows  the  S  to  P  ratio  for  the  composite  event  and  the  shear  crack 
alone.  The  tensile  subevents  decrease  the  S/P  ratio  at  frequencies  above  2  Hz  although 
the  average  effect  on  the  S/P  ratio  is  similar  to  that  of  a  shear  crack  alone  with  a  P  to 
S  comer  frequency  ratio  of  1.73  (see  Figure  18b). 

A  second  composite  event  was  constructed  using  smaller  radius  subevents.  The 
second  composite  model  consists  of  the  same  2  km  radius  shear  crack  with  10  cm 
maximum  slip  plus  1(X)  opening  and  closing  tensile  subevents  which  add  incoherently. 
The  tensile  subevents  have  a  radius  of  2()m  with  a  maximum  opening  of  2  cm  and  a  x 
of  0.015  s.  The  tensile  sub  events  cover  just  one  per  cent  of  the  area  of  the  shear 
crack.  Again  the  total  radiated  tensile  energy  is  about  ten  per  cent  that  of  the  shear 
crack.  The  composite  spectra  are  shown  in  Figure  2()a  as  well  as  that  of  the  shear 
crack  alone  for  compari.son.  The  effect  of  these  smaller  radius  tensile  subevents 
becomes  apparent  for  frequencies  above  20  Hz.  Figure  20b  shows  the  S  to  P  ratio  for 
the  composite  event  and  the  shear  crack  alone.  The  tensile  subevents  substantially 
decrease  the  S/P  ratio  above  20  Hz  and  in  the  80  -  100  Hz  range  the  ratio  is  less  than 
two.  The  large  peak  in  the  S/P  ratio  at  about  70  Hz  is  due  to  the  first  null  in  the 
(2  -  2coso>t)'^^  temi  in  the  opening  and  closing  tensile  subevents. 

A  final  composite  model  u.sing  10(K)  tensile  subevents  with  radii  of  20  m  and 
maximum  openings  of  2  cm  was  constructed.  These  tensile  subevents  cover  just  10 
percent  of  the  area  of  the  shear  crack  but  radiate  approximately  equal  energy.  The  final 
composite  model  spectra  is  shown  in  Figure  21a.  The  effect  of  the  tensile  subevents 
becomes  apparent  above  10  Hz.  Figure  21b  shows  the  S/P  ratio  for  the  composite 
event  and  the  shear  crack  alone.  For  this  model  the  S/P  ratio  drops  below  tw'o  above 
about  30  Hz  and  approaches  the  expected  ratio  of  the  tensile  cracks  alone,  a  value  of 
about  1.2,  at  UK)  Hz. 


Conclusions 

A  model  for  the  displacement  amplitude  spectra  of  a  circular  crack  failing  in  ten¬ 
sion  was  presented.  Although  the  spectra  of  the  tensile  crack  model  is  similar  to  that 
expected  for  circular  shear  cracks,  the  radiation  pattern  differences  lead  to  predictions 
of  much  lower  average  S/P  ratios  for  tensile  as  opposed  to  shear  cracks.  These  low  S/P 
ratios  may  be  an  identifying  characteristic  of  tensile  or  tensile  equivalent  seismic 
sources. 

We  construct  sc'me  composite  events  of  shetu"  and  tensile  cracks  that  may  be 
representative  of  slip  on  a  non-pkinar  fault.  We  show  that  a  tensile  crack  must  have 
about  44  per  cent  of  the  maximum  slip  on  a  shear  crack  to  have  equal  average 
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Fig.  19,  (a)  Composite  speclra  compared  with  pure  shear  crack  spectra. 
Composite  is  composed  of  1  shear  crack  (radius  2  km,  maximum  slip 
10  cm)  and  100  incoherently  adding  opening  and  clo.sing  tensile  cracks 
(20  m  radius,  2  cm  maximum  opening),  (b)  S/P  ratios  of  above  .spectra. 
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Fig.  20.  (a)  Cdmpd.site  .spectra  compared  with  pure  shear  crack.  Composite 
spectra  is  composed  ol  one  shear  crack  (2  km  radius,  maximum  slip  10  cm ) 
and  1000  incoherently  adding  opening  and  closing  tensile  cracks  (20  m 

radius,  2  cm  maximum  opening),  (h)  S/P  ratios  of  above  spectra. 


amplitude  P  spectra.  If  the  tensile  subevents  have  radii  within  an  order  of  magnitude  of 
the  shear  event  the  composite  spectra  may  have  the  appearance  of  a  pure  shear  crack 
with  a  high  P/S  corner  frequency  ratio.  If  the  subevents  have  smaller  radii,  the  compo¬ 
site  spectra  show  much  lower  S/P  ratios  at  high  frequencies  than  expected  for  shear 
cracks  alone.  These  anomalously  low  S/P  ratios  may  occur  when  the  tensile  subevents 
radiate  less  than  10  per  cent  of  the  energy  of  the  shear  crack  and  cover  just  a  fraction 
of  the  fault  area. 


Appendix 

The  spectral  shape  assumed  in  equation  (3.1)  is  one  of  a  class  of  spectral  shapes 
of  the  form: 


I  Q(ci))l 


(1  -I-  ((0/0)^)2)V'2 


(3.A1) 


which  can  be  described  by  a  low  frequency  level,  Q„,  a  corner  frequency,  and  a 
high  frequency  falloff  slope,  y.  P  or  S  waves  which  have  displacement  amplitude 
spectra  given  by  equation  (3.A1)  contain  an  amount  of  energy  given  by: 


oo 

W  =  2  p  R^c  j  I  Q((o)!  ^do)  =  (3.A2) 


where  p  is  the  medium  density.  R  the  distance  from  the  source,  c  the  P  or  S  wave 
velocity  and  X  ^  constant  given  by: 


..  r(y-1.5) 

^  r(v/) 


(3.A3) 


where  T  denotes  the  "gamma  function".  Note  that  equation  (3. A3)  is  defined  only  for 
If  y  =  2.0  then  x  =  '^^  we  get  the  same  results  as  the  main  body  of  the 
paper.  If,  for  example.  \\i  =  2.5  (x  =  4/3V7i)  or  vj/=3.()  (x  =  ’^/4)  then  our  corner  fre¬ 
quency  equations  in  the  main  body  of  the  paper  would  change  but  the  low  frequency 
values  are  unaffected.  Making  use  of  equations  (3.A2)  and  (3. A3)  we  can  rewrite 
equations  (3.17)  and  (3.2.‘i)  for  the  case  of  model  spectra  of  the  form  (3.A1).  Equation 
(3.17)  becomes: 


to; 
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and  equation  (3.25)  bccoiiie.s; 


^  1 

2835  f  ^ 
8  n'3  g'^  +  324 


A 

f  " 


(3.A5) 


It  is  apparent  from  the  above  equations  that  a  high  frequency  falloff  slope  of  three 

(Vj/  =  3)  results  in  corner  frequencies  that  are  4''^=  1.59  times  larger  than  when  the 

falloff  slope  is  assumed  to  be  two. 
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SUMMARY 

The  south  Caspian  Basin  is  a  relatively  aseismic  block  within  the  Alpine-Himalayan 
Belt,  but  is  surrounded  by  zones  of  high  seismicity.  We  used  the  focal  mechanisms  of  16 
earthquakes  whose  source  parameters  w'e  determined  from  inversion  of  body  waves,  and  the 
mechanisms  of  15  other  earthquakes  to  determine  the  style  of  faulting  in  the  seismic  belts 
that  surround  the  south  Caspian  Basin.  Earthquakes  beneath  the  Talesh  Mountains  of  NW 
Iran  and  immediately  off-shore  in  the  SW  Caspian  Sea  have  shallow  thrust  mechanisms, 
showing  that  the  continental  crust  of  NW  Iran  is  overthrusting  the  ” oceanic-like”  crust  of 
the  south  Caspian  Basin.  Earthquakes  south  of  the  Caspian  Sea  in  northern  Iran  show  a 
mixture  of  focal  mechanisms.  Both  high  angle  reverse  faulting  and  left-lateral  strike  slip 
faulting  are  observed  in  the  high  Alborz  Mountains.  Farthers  south,  on  the  edge  of  the 
central  Iran  plateau,  oblique  left-lateral  reverse  mechanisms  are  observed.  It  seems  that  the 
NE  direction  of  shortening  between  central  Iran  and  the  south  Caspian  Basin  is  partitioned 
into  pure  left-lateral  strike-slip  and  thrusting  in  the  WNW  trending  high  Alborz,  but  is 
accomnodated  by  oblique  faulting  in  lower  elevations.  Earthquakes  in  the  Kopet  Dag 
.Mountains  east  of  the  Caspian  Sea  al.^o  ^how  a  mixture  of  high  angle  reverse  and  strike- 
slip  faulting  mechanisms  and  thus  may  be  another  example  of  the  partitioning  of  oblique 
slip  into  strike-slip  and  thrust  motion.  Normal  faulting  mechanisms  at  centroid  depths  of 
35  to  50  km  dominate  in  the  belt  of  seismicity  which  extends  across  the  central  Caspian 
Sea.  The  significance  of  the  normal  faulting  earthquakes  is  enigmatic.  It  is  improbable 
that  these  events  represent  the  motion  between  the  southern  Caspian  Basin  and  Eurasia 
for  they  imply  a  sense  of  motion  which  is  incompatible  with  the  observed  topography  and 
folding  in  the  sediments.  Two  shallow  earthquakes  at  about  12  km  depth  in  this  belt, 
onf'  a  small  evfmt  and  tfie  other  a  larg<'  second  sub-event  of  a  mul,,iple  earthquake,  have 
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thrusting  mechanisms  suggesting  that  shortening  occurs  as  the  continental  crust  of  the 
northern  Caspian  is  thrust  over  the  ”oceanic-like”  crust  of  the  southern  Caspian  Basin. 
Shortening  is  also  suggested  by  the  orientations  of  folds  in  the  sedimentary  cover  south  of 
the  central  Caspian  Sea  seismic  belt.  We  suggest  that  this  shortening  does  indeed  represent 
a  NNE  motion  of  the  Ccispian  Sea  relative  to  Eurasia,  but  that  the  motion  is  slow  and  has 
not  produced  many  earthquakes.  The  deep,  normal  faulting  may  be  related  to  bending 
or  down-dip  extension  of  the  incipient  subducted  slab.  If  the  motion  of  the  Caspian  Sea 
relative  to  Eurasia  is  indeed  slow,  then  the  motion  in  the  Alborz  between  central  Iran  and 
the  south  Caspian  Basin  will  be  almost  the  same  as  that  between  Iran  and  Eurasia,  as  has 
previously  been  assumed. 

The  combined  effect  of  the  overthrusting  of  the  south  Caspian  Basin  by  the  Talesh- 
Alborz  mountains  in  the  south,  and  by  the  continental  crust  of  the  northern  Caspian  Sea  in 
the  north  will  lead  to  the  eventual  destruction  of  the  south  Caspian  Bcisin  and  the  possible 
formation  of  an  intermediate  depth,  dipping  seismic  zone  within  the  continental  interior, 
similar  to  that  presently  observed  in  the  Hindu  Kush. 

INTRODUCTION 

This  paper  is  concerned  with  the  active  faulting  in  and  around  the  south  Caspian  Basin, 
which  is  a  relatively  aseismic  block  within  theAlpine-Himalayan  Belt  but  is  surrounded  by 
highly  active  seismic  zones. 

There  were  two  motivations  for  this  study.  The  first  w.  5  to  discover  the  nature  of 
the  deformation  in  the  belt  of  seismicity  extending  across  the  central  Caspian  Sea  from 
the  Caucasus  in  the  west  to  the  Kopet  Dag  in  the  east.  Before  1986,  no  earthquakes 
of  mj  <  5. Shad  occurred  in  this  belt  since  the  installation  of  the  World  Wide  Standard 
Seismic  Network  (WWSSN)  in  the  early  1960’s  and  thus  no  reliable  fault  plane  solutions 
were  available.  As  a  result,  discussions  of  the  motion  of  the  south  Caspian  Basin  with 
respect  to  Eurasia  have  largely  been  speculative.  Nowroozi  (1972)  believed  that  the  south 
Caspian  is  moving  northwestward  relative  to  Eurasia.  Jackson  &  McKenzie  (1984)  pointed 
out  that  the  thrust  faulting  earthquakes  that  occur  throughout  the  Caucasus,  Alborz,  and 
Kopet  Dag  Mountains  all  have  similar  slip  vectors  directed  NE,  suggesting  that  the  south 
Caspian  is  moving  northeastward  relative  to  Eurasia,  although  less  rapidly  than  Iran.  We 
have  studied  the  mechanisms  of  several  earthquakes  in  the  Trans-Caspian  seismic  belt  to 
better  define  the  motion  of  the  south  Ccispian  Basin. 

A  second  motivation  for  the  study  was  to  better  understand  the  tectonic  setting  of 
the  6  June  1990  earthquake  in  the  Talesh  Mountains  of  the  southwest  Caspian.  This  was 
by  far  the  largest  earthquake  7.2)  to  occur  in  the  region  this  century,  and  it  had  a 
focal  mechanism  quite  unlike  that  of  neighboring  events.  .\  detailed  field  and  scismological 
report  on  this  eartlujuake  is  published  idsewhere  (Berberian  et  al,  1991). 

The  principale  new  data  we  present  new  observations  of  earthquake  focal  mechanisms 
and  depths  determined  from  P-wave  first  motion  and  P-  and  SH-waveform  data.  These 
are  much  better  constrained  than  the  anisms  and  in  earlier  studies  (eg.  Jackson  and 
McKenzie,  1984),  which  were  based  on  P  wave  first  motions  alone,  and  extend  the  period 
of  obsetvation  to  1990. 


67 


GEOLOGICAL  AND  GEOPHYSICAL  SETTING 

The  seismicity  of  the  Middle  Ecist  is  dispersed  over  bands  100  -  300  km  wide,  which 
surround  regions  such  as  the  south  Caspian  Basin  and  central  Ira  within  which  the  seismic¬ 
ity  is  relatively  low  (Fig.  21  and  Jackson  &!:  Mckenzien  1984).  The  western  and  southern 
margins  of  the  south  Caspian  Basin  are  bounded  by  an  active  arcuate  fold  and  thrust  belt 
in  the  Talesh  and  Alborz  Mountains  (Fig.  22a).  To  the  east  of  the  south  Ccispian  Sea  lies 
the  west  Turkmenian  lowlands,  which  are  structurally  part  of  the  south  Caspian  Basin. 
Farther  to  the  east  lie  the  seismically  active  Kopet  Dag  Mountains.  The  deeper  water 
(wlOOOm)  of  the  south  Caspian  Basin  is  separated  from  the  shallower  water  («200m)  of 
the  central  and  north  Caspian  Sea  by  the  Apshceron  -  Balkhan  sill,  a  shallow,  seismi¬ 
cally  active  zone  extending  from  the  Caucasus  Mountains  in  the  west  to  the  Kopet  Dag 
Mountains  in  the  east  (Fig.  21,  22b). 

Deep  seismic  sounding  (DSS)  data  from  the  south  Caspian  Basin  reveal  a  crust  consist¬ 
ing  of  a  thick  (15-20  km)  low  velocity  (Vp  3.5  to  4.0  fcms“*)  sedimentary  section  resting  on  a 
high  velocity  {Vp  >  6.7  kms~'^)  layer  (Galperin  et  al,  1962;  Malovitsky,  1968;  Neprochnov, 
1968;  Rezanov  &  Chamo,  1969;  Malovitsky  et  al,  1970;  Fcdynsky  et  al,  1972)  (Fig.  22a). 
DSS  data  to  the  east  of  the  south  Caspian  Basin  indicate  the  existence  of  a  Vp  6.0  kms~^ 
layer,  thought  to  represent  continental  basement,  which  pinches  out  towards  the  southwest 
and  disappears  beneath  the  Turkmenian  lowlands.  The  existence  of  a  6.0  kms~^  layer  is 
also  apparent  in  DSS  data  from  north  of  the  Apshccron-Balkhan  sill  and  in  the  Talesh  and 
Alborz  Mountains  in  northern  Iran  (Berberian,  1983),  to  the  south  of  the  Ccispian  Sea. 
The  high  velocity  crustal  layer  of  the  .south  Caspian  Basin  has  been  assumed  by  many 
to  be  ”oceanic-like  crust”  and  this  view  is  supported  by  observations  of  Lg  blockage  and 
efficient  propagation  for  paths  crossing  the  south  Caspian  Basin  (Kadinsky  -  Cade,  et 
al,  1981).  Seismic  profiling  shows  that  the  sedimentary  cover  of  the  south  Caspian  Basin 
has  undergone  considerable  deformation  (Zonenshain  Le  Pichon,  1986). 

In  the  western  part  of  the  basin  the  fold  axes  trend  NNW-SSE  to  N-S  (Fig.  22a).  Fold 
axes  in  the  north  and  south  parts  of  the  basin  trend  NW-SE  roughly  parallel  with  the 
Apshceron-Balkhan  Sill  and  the  north  Iranian  coast,  respectively  (Fig.  22a). 

The  tectonic  history  of  the  south  Caspian  Basin  has  been  reviewed  by  Berberian  & 
King  (1981),  Berberian  (1983),  and  Zonenshain  &  Le  Pichon  (1986).  The  age  of  the 
” oceanic”  basement  is  uncertain.  Paleoreconstruction  of  the  Iranian  region  (Berberian 
and  King,  1981)  indicates  that  if  the  high  velocity  layer  is  oceanic  it  could  be  a  relic  of 
an  ocean  that  closed  in  either  the  'I’riassic  or  the  Cretaceous,  or  else  a  marginal  sea  that 
developed  behind  an  island  arc.  The  region  has  been  dominated  by  compression  during  the 
Pliocene- Pleistocene.  Some  Soviet  earth  scientists  (e.g.  Muratov,  1972;  Yanshin  et  al,  1980; 
Shlezinger.  1981)  propose  that  the  south  Caspian  Basin  is  young  and  was  formed  by  the 
eclogitizat  ion  and  subsidence  of  granitic  continental  crust  during  the  mid-  to  late  Cenozoic. 
Others,  (e.g.  Arnursski  et  al,  1968;  Sorokhtin,  1979;  Vardapetyan,  1981)  consider  the  south 
C’aspian  Ihisin  to  be  pro-Hercynian  in  age  and  a  remnant  of  the  Early  Mesozoic  Tethys 
sea  floor.  Evidence  cited  in  support  of  the  latter  hypothesis  includes  the  great  thickness 
of  sediment  and  low  observed  heat  flow.  Shikalibeily  &  Grigoriants  (1980)  suggest  that 
the  South  t'aspian  Basin  is  Jurassic  and  is  overlain  by  Cretaceous  volcanic  rocks.  Adamia 
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et  al  (1977)  and  Letouzey  et  al  (1977)  consider  the  south  Ccispian  Basin  as  a  remnant 
back-arc  basinresulting  from  spreading  behind  a  Late  Cretaceous  -  Paleogene  volcanic  arc. 
Despite  the  uncertainty  in  the  age,  most  authors  agree  that  the  high  velocity  ” oceanic-like” 
basement  extends  eastward  beneath  the  western  Turkmenian  lowlands  and  to  the  western 
coastal  region  of  the  Azerbaijan-Ctispian  coast  (Amursski  et  al,  1968).  A  more  detailed 
discussion  of  the  tectonic  history  of  the  region  surrounding  the  south  Caspian  Basin  can 
be  found  in  Berberian  (1983)  and  Zonenshain  and  Lc  Pichon  (1986). 

DETERMINATION  OF  EARTHQUAKE  FOCAL  PARAMETERS 

We  have  used  P-wave  first  motion  polarities  and  P-  and  SH-waveforms  to  constrain 
the  source  parameters  of  moderate  size  earthquakes  occurring  in  the  region  surrounding 
the  south  Caspian  Basin.  The  earthquakes  whose  mechanisms  we  have  studied  plus  those 
whose  mechanisms  we  have  taken  from  Jackson  &  McKenzie  (1984),  the  Harvard  centroid 
moment  tensor  (CMT)  catalogue,  or  Gao  Wallace  (1991),  arc  listed  in  Table  1. 

The  P-wave  first  motions  shown  in  Figure  23  are  from  WWSSN  seismograms  read  by 
the  authors.  In  almost  all  cases  these  were  read  from  the  long-period  vertical  component. 
In  a  few  cases  where  the  long-period  vertical  onset  was  too  small  or  obscured  by  noise,  we 
read  the  first  motion  from  the  short-period  vertical  component  providing  the  onset  was 
clear,  impulsive,  and  recognizable  as  the  instrument  impulse  response.  In  all  cases,  we 
used  a  P-wave  velocity  at  the  source  of  6.8  fcms"*  to  plot  the  station  positions  on  the  focal 
sphere. 

The  P-  and  SH-waveform  data  are  from  the  WWSSN  and  Global  Digital  Seismic  Net¬ 
work  (GDSN)  stations  in  the  distance  range  30'’  -  95'’.  Our  analysis  followed  the  approach 
described  by  McCaffrey  <k  Nabelek  (1987)  and  subsequently  used  in  a  number  of  other 
studies  (e.g.,  Molnar  Sc  Lyon-Caen  1989;  Taymaz  et  al.,  1990).  Source  parameters  were 
determined  using  McCaffrey  Sc  Abers’  (1988)  version  of  Nabelek’s  (1984)  waveform  in¬ 
version  procedure.  This  routine  minimizes  in  a  least  squares  sense,  the  misfit  between 
the  observed  shape  and  amplitude  of  the  long-period  P-  and  Sll-waveforms  and  synthetic 
waveforms.  The  synthetic  seismograms  are  computed  for  a  ooint  source  embedded  in  a 
simplified  Earth  structure  by  combining  the  direct  arrival  (cither  P  or  S)  with  the  near 
source  reflections  (pP  and  sP,  or  sS)  and  near  source  multiples.  The  earthquakes  were 
assumed  to  be  either  embedded  in  a  velocity  structure  which  was  a  simplified  version  of 
a  published  crustal  velocity  model  for  the  same  region,  or  embedded  in  a  crustal  velocity 
structure  of  a  tectonically  analogous  region  if  no  published  velocity  model  exist  in  the 
region  of  the  earthquake.  Seismographs  are  assumed  to  be  located  on  a  homogeneous 
half  space.  When  available,  we  have  taken  the  Harvard  CMT  solution  reported  in  the 
Preliminary  Determination  of  Epicenters  (PDE)  as  the  starting  point  in  the  inversion; 
whereunavailable,  we  have  used  the  first  motion  fault  plane  solution  of  Jackson  Sc  McKen¬ 
zie  (1984).  The  time  durations  of  the  seismogram  inversion  windows  were  chosen  to  be 
long  enough  to  include  the  near  source  reflected  phases  pP,  sP,  or  sS.  We  examined  the 
waveforms  where  we  thought  the  PeP  or  ScS  arrival  might  be  within  the  inversion  window. 
If  we  thought  these  phases  might  be  of  significant  amplitude,  we  truncated  our  inversion 
window. 
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Synthetic  waveform  amplitudes  are  corrected  for  geometrical  spreading,  and  anelastic 
attenuation  is  accounted  for  by  applying  a  Futterman  (1962)  Q  operator  with  t*p  =  1  sec 
and  r,  =  4  sec.  Uncertainties  in  V  mainly  affect  estimates  of  source  duration  and  seismic 
moment,  but  have  smaller  effects  on  source  orientation  or  centroid  depth  (Fredrich  et  al, 
1988).  To  avoid  complications  introduced  by  the  upper  mantle  triplication  or  core  phase 
interference  we  restricted  the  body  waveform  inversion  (with  one  exception  discussed  in 
the  Appendix)  to  P-waves  in  the  30“  to  90“  distance  range  and  S-waves  in  the  35“  to 
84“  distance  range.  In  the  inversion  process,  SH  seismograms  were  given  only  50%  of  the 
weight  of  the  P-wave  seismograms,  and  all  seismograms  were  azimuthally  weighted;  that 
is,  seismograms  from  stations  clustered  together  in  azimuth  were  given  lower  weight  than 
seismograms  from  isolated  stations.  We  assumed  that  all  slip  occurs  at  the  same  point 
in  space,  the  centroid  location,  but  is  distributed  in  time,  the  source  time  function.  This 
is  usually  a  good  approximation  for  events  less  than  about  M,  6.5.  However,  in  one  case 
complexities  observed  in  the  seismograms  suggest  the  existence  of  a  second  source.  The 
source  time  function  is  described  by  a  series  of  overlapping  isosceles  triangles  (Nabelek, 
1984).  We  initially  selected  the  number  and  duration  of  the  isosceles  triangles  using  the 
event  magnitude  as  a  guide.  We  then  eliminated  late  occurring,  small  amplitude  elements 
of  the  source  time  function  if  we  felt  these  were  the  result  of  noise.  The  inversion  routine 
minimizes  the  misfit  between  the  observed  and  synthetic  waveforms  by  varying  the  strike, 
dip,  slip,  centroid  depth,  seismic  moment,  and  source  time  function. 

The  ”  best-fitting”  set  of  source  parameters  is  found  by  the  inversion  routine.  However, 
the  statistical  description  of  the  inversion  misfit  underestimates  the  true  uncertainties. 
Similar  to  McCaffrey  &  Nabelek  (1987),  once  we  had  obtained  the  minimum  misfit  solution, 
we  performed  a  series  of  tests  to  assess  realistic  uncertainties  and  examine  tradeoffs  between 
the  various  factors.  Our  test  procedure  w2L'i  to  fix  the  source  parameter  being  examined 
at  a  series  of  values  that  bracket  the  ”  best- fitting”  value,  then  reinvert  the  waveform  data 
to  examine  what  effect  variations  in  the  fixed  parameters  had  on  the  free  parameters.  In 
summary,  we  found  that  reasonable  variations  in  crustal  velocity  structure  do  not  result 
in  significant  changes  in  strike,  dip  and  rake,  but  can  affect  centroid  depth,  moment,  and 
the  source  lime  function.  Our  tests  showed  that  the  strike,  dip,  rake,  and  centroid  depth 
are  relatively  independent  of  each  other:  that  is,  when  one  parameter  wcis  held  fixed  at  a 
value  within  a  few  degrees  or  kilometers  of  its  minimum  misfit  value  and  the  waveforms 
reinverted,  the  resulting  values  of  the  free  parameters  were  close  to  their  minimum  misfit 
values. 

AN  EXAMPLE  OF  WAVEFORM  ANALYSIS:  THE  EARTHQUAKE  OF  4  NOVEMBER 
1978 

The  4  .November  1978  (1978.11.04)  Siahbil  earthquake  belongs  to  a  group  of  events 
occurring  in  northwest  Iran  which  we  have  studied,  and  whose  mechanisms  are  consistent 
with  overthrusting  of  the  south  Caspian  Basin  by  the  Azcrbaijan-Talesh  block  of  northwest 
Iran.  This  event  has  been  examined  in  several  previous  studies  and  provides  a  good  example 
of  the  waveform  analysis  conducted  in  this  study.  P-wavc  first  motion  polarities  (Fig.  23) 
constrain  a  347“  striking,  high  angle  (79“  ENE)  plane.  This  does  not  conform  to  any 
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mapped  faults  in  the  area  (Berberian,  1983)  and  is  taken  as  the  auxiliary  plane.  The 
shallow  dipping  fault  plane  is  poorly  constrained  by  the  first  motion  data.  Berberian 
(1983)  studied  P-  waveforms  of  this  event  by  forward  modelling  and  concluded  that  the 
earthquake  occurred  on  a  168“  striking  fault  dipping  9“  WSW  at  a  depth  of  20  ±  4  km. 
One  plane  of  the  best  fitting  double  couple  to  the  Harvard  CMT  solution  strikes  177°,  dips 
9°  WSW,  has  a  rake  of  87°,  and  has  a  centroid  depth  of  15  km.  Figure  24  shows  the  match 
between  the  observed  seismograms  and  synthetic  P  and  SH  seismograms  computed  fe”  ♦he 
minimum  misfit  solution.  We  used  the  Harvard  CMT  solution  as  the  starting  point  in 
the  inversion  procedure.  The  velocity  structure  used  in  obtaining  these  results  consisted 
of  a  12  km  thick  layer  {Vp  5.60  Arms' *,  V,  3.24  Arms  ',  p  2600  fcgm"^)  overlying  a  half¬ 
space  [Vp  6.80  kms'  ^,  V,  3.93  fcms“',  p  2900  kgm~^).  The  inversion  routine  requires  that 
the  source  be  located  in  the  half-space.  The  source  parameters  from  our  inversion  (given 
in  the  title  of  Figure  24  and  in  Table  2)  are  not  significantly  different  from  previously 
determined  solutions  for  this  event.  The  inversion  procedure  is  nonlinear  and  the  result 
may  be  dependent  on  the  starting  parameters.  To  verify  that  we  had  not  found  a  local 
minimum,  we  significantly  displaced  the  starting  solution  from  the  Harvard  solution  and 
reinverted  the  waveform  data  to  verify  that  westill  obtained  a  similar  solution. 

flaving  determined  a  minimum  misfit  solution,  we  performed  a  series  of  tests  to  estimate 
the  effect  on  ttie  source  parameters  produced  by  variations  in  the  source  velocity  structure 
and  to  place  bounds  on  the  source  parameters. 

Source  Rej  ion  Velocity  Structure;  Since  the  crustal  velocity  structure  in  the  south 
Caspian  regio;  is  poorly  known,  our  approach  was  to  use  a  simplified  version  of  a  published 
crustal  veloci  /  models  in  the  area  of  the  earthquake  if  they  exist,  or  the  crustal  velocity 
models  of  an  Analogous  tectonic  region  if  no  crustal  velocity  model  in  the  vicinity  of  the 
earthquake  h»'S  previously  been  determined.  We  performed  a  number  of  tests  for  the 
1978.11.04  earthquake  to  see  how  reasonable  variations  in  source  velocity  structure  affect 
estimates  of  ^he  source  parameters  and  just  how  different  the  source  velocity  structure 
could  be  befo  t  it  significantly  affected  the  inversion  results.  Figure  25  shows  the  results 
of  our  tests,  j'he  velocity  models  tested  are  given  in  the  figure  caption. 

The  McCi  ffrey-Abers  inversion  routine  allows  for  only  a  simplified  source  structure. 
The  first  row  in  Figure  25  shows  the  minimum  misfit  solution  from  Figure  24.  In  our 
first  test  (Rov  B),  we  altered  the  velocity  structure  in  such  a  way  so  as  to  give  the  same 
delay  time  b'  ‘ween  the  direct  phase  and  the  surface  reflection  (the  controlling  factor  in 
determining  t  le  depth),  cis  the  delay  determined  for  the  crustal  velocity  structure  we 
thought  most  "ealistic  (used  in  computing  Row  A).  We  then  reinverted  the  whole  data  set. 
The  synthetic  waveforms  for  the  mini-num  misfit  solution  for  this  velocity  structure  are 
indistinguishable  from  those  in  Fiow  A.  The  focal  parameters  strike,  dip,  rake,  and  centroid 
depth  are  not  significantly  different  from  the  results  in  Row  A  and  only  the  seismic  moment 
is  different.  We  next  made  changes  in  the  velocity  structure  which  still  would  represent 
reasonable  crustal  velocity  structures  (Row  C  -  D)  without  attempting  to  preserve  the 
delay  time  above  the  source,  'fhis  again  produced  negligible  changes  in  the  synthetic 
waveforms  from  those  shown  in  How  A.  Row  E  and  F  sfiow  minimum  misfit  synthetics  for 
more  radical  (and  probably  unrealistic)  changes  in  the  crustal  velocity  structures.  These  do 
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show  more  significant  changes  but  even  in  these  crises,  the  source  orientations  have  hardly 
changed  from  those  in  Row  A.  We  conclude  from  these  tests  that  reasonable  variations 
in  crustal  velocity  structure  affect  strike,  dip,  and  rake  less  than  other  factors.  Varying 
the  source  velocity  structure  did  have  a  significant  effect  on  centroid  depth  and  seismic 
moment.  These  arise  from  variations  in  delay  times  between  the  direct  and  surface  - 
reflected  phzises,  and  from  variations  in  the  reflection  coefficients. 

Assessment  of  Source  Parameter  Uncertainties:  Thorough  discussions  of  the  assignment 
of  uncertainties  are  given  in  McCaffrey  &  Nabelek  (1987)  and  Molnar  Lyon-Caen  (1989). 
For  each  event  we  studied,  we  have  made  extensive  tests  to  determine  uncertainty  bounds 
on  thecentroid  depth,  and  the  strike,  dip,  and  rake  of  the  fault  plane.  These  uncertainties 
are  listed  in  Table  2. 

Figure  26  summarizes  only  the  tests  we  have  made  to  place  bounds  on  the  strike,  dip, 
rake  and  centroid  depth  of  the  1978.11.04  event.  The  first  row  compares  the  observed 
seismograms  and  those  computed  for  the  minimum  misfit  solution.  The  following  four 
rows  compare  the  observed  seismograms  with  those  computed  when  one  source  parameter 
is  held  fixed  at  the  smallest  displacement  from  the  minimum  misfit  value  which  resulted 
in  a  significant  degradation  of  the  synthetic  fit.  In  making  the  full  range  of  tests,  we  held 
one  parameter  (e.g.  centroid  depth)  fixed  at  a  value  different  from  the  minimum  misfit 
value,  and  inverted  for  the  remaining  source  parameters.  We  then  successively  increased 
the  difference  between  the  fixed  parameter  we  were  testing  and  its  minimum  misfit  value 
until  a  significant  difference  was  apparent  between  the  observed  seismograms  and  the 
synthetics  computed  for  the  particular  combination  of  source  parameters.  Based  on  a 
range  of  tests  of  these  parameters  which  are  summarized  in  Figure  26,  we  estimate  that 
for  the  1978.11.04  earthquake  the  bounds  on  the  centroid  depth  are  ±  5  km,  strike  ±15°, 
dip  1:4'",  and  rake  ±20°.  For  low-angle  faulting  events  such  as  the  1978.11.04  earthquake 
the  rake  is  not  lightly  constrained.  The  bounds  on  focal  parameters  shown  in  Table  2  are 
of  similar  magnitude  to  those  noted  in  comparable  studies  in  other  regions  (e.g.  Molnar 
A:  Lyon-Caen,  1989;  Tayrnaz  et  al.,  1990). 

THE  EARTHQUAKE  OF  6  MARCH  1986 

We  now  discuss  the  earthquake  of  1986.03.06  in  some  detail  for  two  reasons:  it  was 
the  first  moderately  large  earthquake  to  have  occurred  in  the  central  Caspian  seismic  zone 
since  the  installation  of  the  WWSS.N,  and  it  was  the  only  event  we  studied  whose  seismo¬ 
grams  suggested  to  us  the  possibility  of  a  complex  multiple  rupture  involving  sub-events 
of  different  mechanisms.  The  S-wave  spectrum  from  seismograms  recorded  at  Askhabad 
(Golensky  et  al.,  1989)  and  the  coda  spectrum  (Rautian  -  personal  communication)  are 
complex,  showing  a  constant  low  frequency  level,  a  low  frequency  corner  at  0.08  Hz  above 
which  there  is  a  high  frequency  decay  of  to  a  second  high  frequency  corner  at  1.0  Hz 
above  which  the  spectrum  decays  as  ■^.  Smith  et  al  (1991)  attribute  such  a  spectral 
shape  to  source  complexity  with  the  low  frequency  corru'r  corresponding  to  the  overall 
source  dimension  and  the  high  fre(|iien<  y  corner  relatt'd  in  some  fashion  to  the  size  of  the 
Hubevents  d'he  P-wave  first  motion  polarities  (Fig.  23)  tightly  constrain  a  near  vertical 
nodal  plai.e.  striking  300  .  Our  single  source  minimum  misfit  inversion  solution  is  shown 


in  Figure  27a.  This  solution  is  very  similar  to  the  best  fitting  double  couple  solution  from 
thellarvard  CMT  catalogue.  The  agreement  between  the  observed  and  synthetic  seismo¬ 
grams  is  good  at  most  stations.  However,  waveforms  of  stations  lying  on  the  focal  sphere 
near  the  P-wave  nodal  plane  (e.g.  SHL,  KOD,  PTO,  STU)  show  a  weak  initial  P-wave 
arrival  in  keeping  with  the  nodal  position,  but  also  show  a  large  second  pulse,  shaped  like 
an  inverted  "W”,  that  follows  the  first  motion  by  about  20  seconds.  The  single  source 
shown  in  Figure  27a  cannot  reproduce  the  late  pulse  observed  at  these  stations  since  the 
direct  P  and  pP  are  both  nodal,  and  hence  of  small  amplitude,  and  sP  leaves  the  focal 
sphere  close  to  the  node  in  the  S  radiation  pattern. 

With  our  poor  knowledge  of  the  undoubtedly  complex  crustal  velocity  structure  in  the 
vicinity  of  the  Apshceron-Balkhan  sill,  it  is  impossible  to  adequately  cLssess  the  contribution 
to  seismogram  complexity  resulting  from  three-dimensional  velocity  structure.  However, 
what  appears  to  be  missing  from  the  P-waveforms  of  the  near  nodal  seismograms  is  an 
inverted  ”W”  P-wave  arrival  reminiscent  of  shallow  depth,  reverse  faulting  earthquakes 
(e.g.,  Jackson  &  Fitch,  1981).  This  suggests  to  us  that  the  6  March  1986  earthquake 
may  have  been  a  multiple  event.  Source  complexity  is  also  suggested  by  the  shape  of  the 
source  spectrum  at  regional  stations  in  the  U.S.S.R.  (Golensky  et  al.,  1989;  Ratuian  - 
personal  communication).  We  therefore  reinverted  the  P-  and  SH-waveforms,  allowing  for 
the  existence  of  a  second,  delayed,  reverse-faulting  event.  The  starting  source  model  for  the 
second  waveform  inversion  consisted  of  an  initial  event  with  the  source  mechanism  shown 
in  Figure  27a,  followed  10  seconds  later  by  a  second  reverse  faulting  event.  The  starting 
source  orientation  for  the  second  reverse  faulting  event  was  guided  by  our  interpretation 
of  the  observed  waveforms.  The  reversed  mechanism  must  be  constrained  to  produce  the 
observed  P-waveforms  at  nodal  stations  such  as  SHL,  KOD,  PTO  and  STU,  but  not  disturb 
the  observed  P-waveforms  at  African  stations  such  as  SLR  and  WIN,  nor  disturb  the  nodal 
character  of  the  SH-wavoforms  at  these  stations.  A  45"  reverse  faulting  event  with  a  strike 
near  300",  similar  to  the  strike  of  the  initial  event,  conforms  to  these  requirements.  We 
reinverted  the  waveform  data,  first  constraining  the  mechanisms  of  the  first  event  and 
solving  for  the  fault  plane  parameters,  centroid  depth,  source  time  function,  and  origin 
time  offset  of  the  second  event.  To  produce  the  final  solution  shown  in  Figure  27b,  we 
inverted  the  waveform  data  with  all  parameters  of  both  events  free.  Synthetic  waveforms 
for  this  double  event  solution  fit  observed  waveforms  at  all  azimuths  equally  well.  The 
nodal  character  of  SH-waves  observed  at  African  stations  is  preserved  and  provides  the 
main  constraints  on  the  strike  of  the  second  event. 

To  assure  ourselves  that  the  initial  conditions  we  assumed  in  constructing  the  second 
source  were  not  overly  restrictive  and  biasing  our  result,  we  reinverted  the  waveform  data 
with  both  events  having  a  starting  solution  similar  to  the  normal  faulting  mechanism  shown 
in  Figure  27a,  but  the  second  source  delayed  by  10  seconds.  The  mechanism  for  the  first 
event  changed  little,  and  the  solution  of  the  second  event  changed  from  having  a  normal 
to  reverse  mechanism.  We  reinverted  the  data  with  other  startingdepths  for  the  second 
event  and  again  the  final  result  was  very  near  the  sedution  shown  in  Figure  27b. 

Figure  28  compares  results  for  the  single  and  (!<■  ;ble  source  solution  and  individual 
waveforms  for  the  elements  of  the  double  source.  The  upper  row  compares  the  observed 
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seismograms  with  those  computed  for  the  single  source,  minimum  misfit  solution.  Stations 
SHL,  KOD,  PTO,  and  STU  lie  on  the  focal  sphere  near  the  P-wave  node  for  the  single 
source;  station  WIN  lies  near  an  SH  node.  The  single  source  synthetics  poorly  match  the 
observed  seismograms  at  the  four  P-wave  nodal  stations,  as  discussed  above.  The  P-wave 
seismogram  recorded  at  SLR  is  reasonably  well  fit  by  the  synthetic,  with  the  exception 
of  the  large  back  swing.  The  SH  nodal  character  observed  in  the  WIN  seismogram  is 
preserved  and  the  SH  waveform  observed  at  KOD  is  recisonably  well  fit  with  the  exception 
of  the  large  back  swing  which  is  larger  in  the  synthetic  than  in  the  observed  seismogram. 

The  second  row  of  seismograms  compares  observed  waveforms  at  the  same  stations  with 
synthetic  seismograms  computed  for  the  double  source  solution.  The  fit  at  WIN  is  similar 
to  that  for  the  single  source,  and  fits  for  all  other  stations  are  significantly  improved.  Rows 
3  and  4  show  the  contributions  to  the  seismograms  of  the  individual  components  of  the 
double  source.  Both  events  have  similar  seismic  moments.  We  performed  a  series  of  tests 
as  described  above,  to  estimate  bounds  on  the  centroid  depth,  strike,  dip,  and  rake  of  the 
second  event  on  the  synthetic  and  compare  these  with  the  observed  seismogram.  From  this 
comparison  we  estimate  that  the  centroid  depth,  strike,  dip,  and  rake  of  the  second  event 
are  resolved  to  :;.4  km,  ±15°,  ±5°,  and  ±8°,  respectively.  Note  that  the  seismic  moment 
of  this  second  source  is  almost  as  great  as  that  of  the  first  source.  The  combined  seismic 
moment  of  the  two  subevents  is  5.1X10^^  Nm  compared  to  6.4X10'^  Nm  for  the  Harvard 
CMT  solution. 

It  is  difficult  to  assess  whether  the  three  dimensional  velocity  structure  in  the  source 
region  of  the  1986.03.06  earthquake  could  produce  the  seismogram  complexity  that  we 
attribute  to  a  later  sub-event.  Two  other  moderate-sized  earthquakes  occurred  in  the  cen¬ 
tral  Caspian  Sea  on  16  and  17  September  1989  (1989.09.16,  1989.09.17)  and  our  waveform 
inversion  results  for  these  two  events  are  shown  in  Figures  A 13  and  A 14.  P-waveforms  of 
1986.03.06.  1989.09.16,  and  1989.09.17  at  common  stations  are  compared  in  Figure  29.  As 
a  result  of  the  decline  in  performance  of  the  WWSSN,  data  for  the  1989  events  are  less 
comph^te  than  for  1986.03.06.  The  minimum  misfit  solution  for  both  of  the  1989  events 
have  a  normal  faulting  mechanism  similar  to  the  single  source  minimum  misfit  solution 
for  1986.03.06.  However,  WWSSN  seismograms  recorded  to  the  northwest  (e.g.,  PTO, 
VAL.  NCR.  WES),  while  not  showing  the  same  initial  nodal  character  as  seismograms 
recorded  at  similar  azimuths  for  the  6  March  1986  event,  do  show  some  evidence  for  a 
similar  later  pulse  in  iheP-wave  train  (Figs.  A13,  A14).  Although  the  GDSN  instrument 
response  tends  to  blur  details  of  the  source  time  function  and  WWSSN  data  are  scarce 
for  the  1089  events,  we  inverted  seismograms  for  the  1989.09.16  earthquake  allowing  the 
possibility  of  a  second  source.  This  resulted  in  a  mechanism  where  the  first  source  was 
similar  to  that  in  Figure  .\13,  and  a  second  source  whose  seismic  moment  was  less  than 
10%  that  of  the  initial  source.  Three  double  earthquakes  with  similar  time  delays  and 
mechanisms  for  the  second  events  are  somewhat  suspect,  and  this  may  imply  that  the 
secondary  .trrivals  result  from  crustal  velocity  structure  complexities  in  the  source  region. 
How(!ver  in  a  plane  layered  structure  there  is  no  direct  P  or  SH  arrival  for  the  European 
and  Indian  stations  showing  the  second  pulse.  It  would  therefore  be  necessary  for  the  late 
pulse  to  he  the  result  of  bent  raypaths.  that  is  rays  refracted  after  leaving  the  focal  sphere 
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at  different  take  off  angles.  On  the  other  hand,  there  is  another,  smaller  earthquake  in 
this  central  Caspian  belt  (1987.12.19)  with  a  Harvard  CMT  solution  having  almost  the 
same  reverse  faulting  mechanism  as  we  postulate  for  the  second  sub-event  of  1986.03.06. 
The  earthquake  of  1987.12.19  was  too  small  for  us  to  be  able  to  confirm  the  Harvard  CMT 
solution  using  long  period  P  and  SH  waveforms.  Although  we  prefer  the  double  source 
explanation  for  the  P-wave  complexities  discussed  above  for  1986.03.06,  we  are  not  able 
to  rule  out  the  possibility  that  complexities  in  the  observed  waveforms  are  the  result  of 
source  region  velocity  structure. 

TECTONICS  OF  THE  SOUTH  CASPIAN  BASIN 

The  minimum  misfit  solutions  for  the  other  fourteen  events  we  studied  are  presented  in 
Appendix  A.  Their  orientations  are  consistent  to  within  a  few  degrees  with  the  first  motion 
polarity  data.  The  few  inconsistent  first  motion  polarities  are  discussed  in  Appendix  A. 
Figure  30a  shows  the  fault  plane  solutions  for  these  events  (Table  2).  Figure  30b  shows  fault 
plane  solutions  for  other  events  which  we  have  not  analyzed:  either  first  motion  solutions 
of  Jackson  &  McKenzie  (1984)  which  occurred  prior  to  establishment  of  the  WWSSN 
(1957. 07.02)  or  which  produced  insufficient  long-period  waveforms  for  us  to  analyze  using 
the  procedure  outlined  above.  Harvard  CMT  solutions  that  were  too  small  for  us  to  study 
their  P  and  SH  waveforms,  or  the  20  June  1990  event  (Gao  &  Wallace,  1991).  With  some 
exceptions,  the  patterns  of  focal  mechanisms  shown  in  Figures  30a  and  30b  are  generally 
similar,  giving  us  some  confidence  in  the  solutions  we  were  unable  to  confirm  by  waveform 
inversion.  In  our  discussion  of  the  tectonics,  we  will  rely  primarily  on  the  sixteen  events 
we  have  studied  (Fig.  30a)  and  the  mechanism  of  the  20  June  1990  event  (Gao  Wallace, 
1991)  (large  symbol.  Fig.  30b),  but  we  will  make  cautious  use  of  the  firstmotion  and  CMT 
solutions  of  the  smaller  events. 

Slip  Vectors:  Figure  31  shows  the  horizontal  projection  of  the  slip  vectors  for  the  events  in 
Figure  30.  Choosing  the  slip  vector  required  us  to  identify  one  of  the  two  nodal  planes  in 
each  fault  plane  solution  as  the  fault  plane.  Both  the  Buyin  Zara  earthquake  (1962.09.01) 
(Appendix  A)  and  the  fludbar-Tarom  earthquake  (1990.06.20)  (Berberian  et  al,  1991) 
were  accompanied  by  surface  faulting,  and  the  choice  of  the  fault  plane  is  straightforward 
(Ambraseys,  1963;  Berberian  et  al,  1991).  No  surface  faulting  wtis  recognized  for  the 
Karnaveh  earthquake  (1970.07.30).  However,  relative  relocations  of  the  aftershocks  of 
this  event  defined  an  elongated  pattern  trending  north-northeast  (Jackson  &:  Fitch,  1979), 
approximately  parallel  to  the  strike  of  one  of  the  nodal  planes  (Figure  A2)  and  also  parallel 
to  a  nodal  plane  of  a  large  aftershock  with  a  similar  mechanism  to  that  of  the  main  shock. 
The  aftershock  pattern  suggests  that  the  Karnaveh  eartlujuake  occurred  on  a  northeast¬ 
trending  fault  with  predominantly  left-lateral  strike-slip  motion. 

Other  earthquakes  such  as  the  1957.07.02,  1969.01.03,  and  1983.07.22  events  are  almost 
pure  thrust  and  hence  the  choice  of  the  slip  vector  orientation  is  not  very  sensitive  to  the 
choice  of  the  fault  plane.  For  the  1969.01.03  and  1983.07.22  events  we  chose  the  planes 
dipping  northeast  to  be  consistent  with  the  local  topography;  however,  in  both  cases 
choosing  the  other  plane  would  change  the  orientation  of  the  slip  vector  by  less  than 
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10°.  The  other  events  studied  all  have  one  nodal  plane  that,  if  chosen  as  the  fault  plane, 
implies  a  slip  vector  similar  to  that  of  an  adjacent  earthquake  whose  slip  vector  ambiguity 
is  resolved.  This  determined  our  choice  of  slip  vectors  for  those  events. 

In  the  case  of  strike-slip  events,  this  reasoning  is  dangerous.  The  1990.06.20  event, 
which  produced  surface  faulting,  had  a  slip  vector  almost  perpendicular  to  that  of  a 
nearby  reverse  faulting  earthquake  (1983.07.22)  (See  Berberian  et  al,  1991;  Gao  &  Wallace, 
1991).  Three  other  strike-slip  events  occurred  near  thrust  faulting  earthquakes  (1983.03.26, 
1984.02.22,  and  1990.01.20).  For  all  of  these  we  mark  both  possible  slip  vectors  on  Figure 
31.  The  convention  we  have  used  in  Figure  31  is  to  show  the  motion  of  the  S  or  E  side 
relative  to  the  N  or  W  side. 

The  inset  in  Figure  31  shows  the  orientation  of  the  tension  axes  of  the  five  normal 
faulting  events  in  the  central  Caspian.  All  of  these  events  have  a  shallow  dipping,  NNE 
striking  tension  eixis. 

Talesh  Mountains:  It  is  clear  from  the  earthquake  mechanisms  we  have  studied  that  the 
Talesh  and  NW  Alborz  mountains  overthrust  thesouth  Caspian  Basin  as  suggested  by 
Berberian  (1983)  and  Jackson  &  McKenzie  (1984).  The  slip  vectors  on  the  shore  of  the 
Talesh  Mountains  show  the  motion  of  the  Caspian  Sea  westward  beneath  the  Azerbaijan- 
Talesh  block,  which  moves  east  relative  to  both  Iran  and  the  Arabian-Eurasia  collision 
zone  in  eastern  Turkey.  This  motion  of  the  Azerbaijan-Talesh  block  may  be  accommodated 
inland  by  conjugate  strike-slip  faulting.  The  focal  mechanism  of  the  1976.02.03  earthquake 
(Fig.  30b)  is  consistent  with  left  lateral  faulting  on  a  NNE  plane,  parallel  to  a  prominent 
NNE  lineament  followed  by  the  Araxes  River.  Berberian  &  Arshadi  (1976),  Tchalenko 
(1977),  and  Ambraseys  &  Melville  (1982)  have  studied  faulting  associated  with  two  large 
northwest  Iranian  earthquakes  which  occurred  in  1721  and  1786  (Fig.  22b).  All  agree  that 
the  southwest  side  of  the  fault  w«ls  down-dropped,  but  disagree  as  to  whether  there  was 
a  significant  right-lateral  strike-slip  component.  Stream  diversions  are  visible  on  aerial 
photographs  (Tchalenko,  1977)  suggesting  right  lateral  offset  but  these  have  not  been 
verified  on  the  ground.  The  combination  of  left  lateral  faulting  following  the  Araxes 
lineament  along  the  Iran-USSR  border,  probable  right  lateral  faulting  extending  southeast 
from  the  North  Anatolian  fault  system,  and  low  angle  thrusting  along  the  Talesh  Mountains 
support  Jackson  &  McKenzie’s  (1984)  suggestion  of  the  easterly  expulsion  of  northwest 
Iran  away  from  the  collision  zone  in  the  Caucasus  and  the  overthrusting  of  the  south 
Caspian  Basin.  This  motion  is  probably  responsible  for  the  N-S  trending  folds  seen  in  the 
sedimentary  cover  offshore,  east  of  the  Talesh  Mountains  (Fig.  22a),  and  is  presumably 
favored  by  the  "oceanic-like”  crust  of  the  south  Caspian  Basin. 

Alborz  Mountains:  The  Rudbar-Tarom  earthquake  (1990.06.20)  was  by  far  the  largest 
earthquake  to  occur  in  northern  Iran  this  century.  Berberian  et  al  (1991)  report  at  lerist 
80  km  of  surface  faulting  along  three  discontinuous  fault  segments.  These  are  arranged 
in  a  right-stepping  en-echelon  pattern,  and  are  separated  by  gaps  in  the  observed  surface 
rupture.  The  three  segments  range  in  strike  from  95°  to  120°  and  have  oblique,  left-lateral 
and  reverse  motion  on  faults  that  are  sub-vertical  or  have  steep  dips  to  the  S  or  SSW.  The 


maximum  surface  displacements  were  60  cm  horizontal  (left-lateral)  and  95  cm  vertical  with 
the  south  side  up.  The  moment  tensor  solution  (Gao  &  Wallace,  1991)  for  the  main  shock 
suggests  a  complex  rupture  consisting  of  at  least  three  sub-events,  and  had  a  significant 
non-double  couple  component.  However,  no  change  in  fault  geometry  or  slip  direction 
during  the  main  rupture  could  be  resolved.  The  best  fitting  double  couple  (Table  3)  has  a 
strike  292°,  dip  88°,  and  rake  —9°  which  are  consistent  with  the  mapped  faults.  Berberian 
et  al  (1991)  also  reported  mechanisms  for  two  aftershocks.  The  largest  (1990.06.21)  had 
a  focal  mechanisminvolving  thrusting  (strike  26°,  dip  69°,  rake  87°)  with  a  centroid  depth 
of  17  km  and  a  moment  of  5.7X10*^  Nm  (M^v5.8)  similar  to  those  in  the  adjacent  Talesh 
Mountains.  A  second  large  aftershock  (1990.06.24)  had  a  mechanism  similar  to  that  of  the 
main  shock  (strike  235°,  dip  70°,  rake  -164°),  a  significantly  shallower  centroid  depth  (8 
km),  and  a  moment  of  2.2X10’^  Nm  (Mvv5.5). 

The  1990.06.20  earthquake  occurred  close  to  an  earthquake  that  involved  high-angle 
reverse  faulting  (1983.07.22).  The  slip  vector  of  the  1990.06.20  main  shock  is  in  the  wrong 
sense  (left-lateral)  to  represent  the  relative  motion  between  the  Talesh-Azerbaijan  block 
and  Iran.  We  think  it  probable  that  the  deformation  in  this  part  of  the  NW  Alborz 
Mountains  involves  oblique,  left-lateral  shortening  that  is  partitioned  into  pure  left-lateral 
strike-slip  (i.e.  1990.06.20)  and  pure  thrusting  (i.e.  1983.07.22)  in  almost  the  same  region. 
Such  an  interpretation  is  consistent  with  the  almost  perpendicular  slip  vectors  for  the 
1990.06.20  and  1983.07.22  earthquakes,  and  with  those  slip  vectors  being  approximately 
parallel  and  perpendicular,  respectively,  to  the  regional  topography  and  fold  axes  (see 
also  Berberian  et  al,  1991).  Jackson  A:  McKenzie  (1984)  suggested  that  the  more  easterly 
oriented  slip  vectors  observed  further  south  (e.g.  the  slip  vector  of  42°  in  1962.09.01) 
represent  the  direction  of  regional  convergence  between  Iran  and  South  Caspian.  As  a 
consequence,  there  should  be  a  component  of  left-lateral  slip  in  the  Alborz  Mountains.  In 
the  high  Alborz  Mountains,  it  appears  that  this  oblique  regional  convergence  is  partitioned 
into  strike-slip  (1990.06.20)  and  reverse  faulting  (1983.07.22)  wherecis  it  remains  as  oblique 
slip  in  the  lower  elevations  to  the  south  in  north-central  Iran.  The  component  of  left-lateral 
slip  in  the  Alborz  Mountains  is  small,  which  may  account  for  the  apparent  infrequency  of 
strike-slip  faulting  earthquakes  in  the  range. 

The  confirmed  existence  of  adjacent  strike-slip  and  thrust  mechanisms  with  nearly 
perpendicular  slip  vectors  in  the  NW  Alborz  Mountains  makes  us  unsure  of  the  actual 
slip  vector  for  the  strike  slip  events  near  thrust  events  elsewhere  in  the  Alborz  Mountains, 
which  is  why  we  include  both  possible  slip  vectors  for  1983.03.26  and  1990.01.20.  If  the 
EW  nodal  planes  were  the  fault  planes  in  these  earthquakes,  their  slip  vectors  are  similar 
and  in  the  same  sense  as  for  the  1990  Rudbar-Tarom  earthquake.  In  this  case,  these  events 
too  are  evidence  for  partitioning  of  strike-slip  and  thrust  motion  in  the  Alborz.  South  of 
the  high  Alborz  Mountains,  in  the  Central  Iranian  Plateau,  the  motion  is  apparently  not 
partitioned  and  occurs  as  faults  with  oblique  slip  (e.g.  1962.09.01,  1980.12.19,  1980.12.22). 

Most  slip  vectors  in  the  Alborz  are  directed  NE  (Fig.  31)  consistent  with  convergence 
between  central  Iran  and  the  south  Caspian  Basin,  .^part  from  the  strike  slip  events, 
which  we  havediscussed  above,  the  only  exception  is  a  small  (Afvv5.l)  normal  faulting 
event  (1988.08.23)  for  which  only  a  Harvard  CMT  solution  is  available  (Fig.  30b).  We 
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are  uncertain  how  to  assess  this  event  and  cannot  evaluate  the  reliability  of  the  solution. 
There  axe  no  other  normal  fault  solutions  like  it  in  the  Aiborz  (The  1970.07.30  Karnaveh 
earthquake,  although  it  had  a  normal  component,  had  a  slip  vector  consistent  with  the 
NE  motion  between  Iran  and  the  south  Caspian  Basin).  Below  a  certain  threshold  in 
magnitude  it  is  quite  common  for  small  events  to  have  a  greater  variety  of  mechanisms 
than  the  larger  events  (e.g.  Lyon  -  Caen  et  al,  1988).  Taymaz  et  al  (1991)  report  a  well- 
constrained  reverse-faulting  event  of  Mw5.5  in  a  region  of  NW  Turkey  that  is  dominated 
by  normal  faulting.  With  only  one  anomalous  event  (that  is  probably  not  well  constrained) 
in  the  Aiborz,  we  cannot  cissess  its  significance  further.  The  component  of  normal  faulting 
seen  in  1970.7.30  suggests  that  the  low  elevations  and  thick  sediments  in  the  Turkmenian 
Lowlands,  (Fig.  22a)  may  be  partly  related  to  extension-induced  subsidence. 

The  Kopet  Dag  Mountains:  The  Kopet  Dag  has  an  abrupt  northeast  topographic  front 
(Fig.  22a)  that  was  associated  with  a  large  earthquake  in  1948  near  the  city  of  Askhabad. 
Field  and  teleseismic  data  for  this  earthquake  are  ambiguous  (Rustanovich  &:  Shirakova, 
1964;  Tchalenko,  1975;  and  Jackson  &  McKenzie,  1984),  but  several  authors  associate 
the  linearity  of  the  abrupt  topographic  scarp  with  right-lateral  strike-slip  motion,  whereas 
elsewhere  in  the  Kopet  Dag  there  is  certainly  thrust  faulting  with  northecist  directed  slip 
vectors.  One  other  event  (1984.02.22,  Figure  3Ca)  further  NW  is  consistent  with  right  lat¬ 
eral  strike-slip  on  this  trend.  The  Kopet  Dag  may  thus  be  another  example  of  partitioning 
of  strike-slip  and  thrust  motion  that  are  postulated  to  occur  in  the  Aiborz  Mountains. 

Central  Ctispian  Seismic  Belt:  In  the  belt  of  seismicity  which  extends  across  the  central 
Caspian  Sea  following  the  Apshceron  -  Balkan  sill,  normal  faulting  mechanisms  with  an 
ESE  strike  and  focal  depths  of  30-50  km  dominate.  It  is  improbable  that  the  normal 
faulting  events  represent  the  motion  between  the  southern  Caspian  and  Eurasia.  If  the 
steep  nodal  planes  were  the  fault  planes,  then  the  sense  of  motion,  north  side  down,  is  the 
reverse  of  that  seen  in  the  bathymetry,  i.e.,  deep  water  and  thick  sediments  to  the  south;  if 
the  nodal  planes  dipping  gently  south  were  the  fault  planes,  some  expression  of  extension 
might  be  expected  at  the  surface.  Yet  seismic  reflection  surveys  report  folding  of  the  surface 
sediments  (Zonenshain  A’  Le  Pichon,  1986),  suggesting  shortening.  This  is  compatible  with 
the  high  angle  reverse  faulting  mechanism  that  we  postulate  for  a  possible  second  sub-event 
of  the  1986.03.06  earthquake  at  a  shallower  depth  than  the  normal  faulting  earthquakes. 
The  C.MT  solution  for  the  1987.12.19  earthquake  (Fig. 30b)  shows  reverse  faulting  with  a 
NW-SE  strike,  and  supports  this  interpretation.  However,  it  was  a  small  (Mw4.9)  event, 
and  in  the  light  of  the  anomalous  small  event  (1988.08.23)  in  the  Aiborz  Mountains  (Fig. 
30b),  we  do  not  place  much  confidence  in  it. 

We  suspect  that  shortening  occurs  sls  the  continental  crust  of  the  northern  Caspian  is 
thrust  over  the  ’’ocean’c”  crust  ''  the  southern  Caspian  Basin,  and  that  the  normal  faulting 
events  occur  in  tlie  deep  (possible  .mantle)  basement  of  the  material  being  overthrust.  It 
is  possible  that  a  subducting  mantle  slab  with  very  low  seismicity  dips  north  beneath  the 
.\pshreron-Halkhan  sill,  accounting  for  the  north-dipping  T  axes  of  the  normal  faulting 
solutions.  .Vnother  po.ssibilit v  is  that  the  normal  faulting  is  the  result  of  plate  bending 
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exceeding  the  elastic  strain  limit  as  the  north  side  of  the  Caspian  Sea  is  thrust  over  the 
south  Caspian  Basin.  This  would  be  analogous  to  the  normal  faulting  seaward  of  the 
oceanic  trench  and  is  a  possible  cause  of  the  normal  faulting  beneath  the  Ganges  Basin 
south  of  the  Himalayas  thrust  front  (Lyon-Caen  At  Molnar,  1985).  The  low  level  of  seismic 
shortening  across  the  central  Caspian  suggests  that  most  of  the  Iran-Eurasia  shortening  is 
accommodated  in  the  much  more  active  Alborz  thrust  belt.  The  seismicity  of  the  trans- 
Caspian  belt  continues  northwest  into  the  Caucasus  and  southeast  into  the  Kopet  Dag. 
Reverse  faulting  and  thrusting  on  ESE  striking  planes  are  dominant  in  both  of  these  areas. 

DISCUSSION 

The  two  most  important  new  results  that  our  study  demonstrates  are  the  partitioning  of 
thrust  and  strike  slip  motion  in  the  Alborz,  and  the  existence  of  a  belt  of  normal  faulting 
earthquakes  acrcos  the  central  Caspian  Sea.  Both  have  important  implications  for  the 
tectonics  of  the  region. 

The  partitioning  of  oblique  shortening  into  pure  thrust  and  pure  strike-slip  motion  is 
not  particularly  rare,  though  there  is  some  debate  as  to  why  it  occurs  (e.g.  McKenzie 
&  Jacksoii  1983,  Zoback  et  ah,  1987,  Mount  Ar  Suppe  1987,  Molnar  1991).  McKenzie  A: 
Jackson  (1983)  point  out  that,  if  fault-bounded  blocks  move  in  response  to  forces  imposed 
on  their  bases  by  continuous  flow  in  the  mantle  lithosphere  beneath  the  seismogenic  upper 
crust,  then  oblique  slip  on  faults  oriented  obliquely  to  the  strike  of  the  deforming  zone 
cannot  accommodate  large  finite  deformation  because  they  rotate  about  a  vertical  axis  as 
they  move.  One  way  in  which  the  faults  can  take  up  large  finite  strains  is  if  the  motion 
on  them  is  partitioned  into  pure  shortening  and  pure  strike  slip.  A  clue  to  this  process 
may  lie  in  the  observation  that  the  motion  is  partitioned  in  the  high  Alborz,  where  strains 
are  obviously  large,  but  not  in  the  lower  elevations  of  north-central  Iran,  wher>.  strains 
are  presumably  smaller.  The  faults  may  have  become  organized  into  a  stable  geometry  as 
the  strain  increased.  This  may  explain  why  partitioning  is  common  in  other  places  where 
-strainsare  large,  such  as  in  California  (e.g.  Mount  &  Sui)pe,  1987),  the  NW  Zagros  and 
NE  Kopet  Dag  Mountains  of  Iran  (Tchahuiko  &  Braud,  1971;  Tchalenko,  1975;  Jackson  A: 
McKenzie  ,  1984)  and  in  island  arcs  (e.g.  Eitch,  1972),  but  oblique  slip  is  by  no  means  rare 
in  places  wh'^re  total  strain  is  probably  r<datively  small,  such  as  in  the  lower  elevations  of 
central  Iran  (Jackson  A:  McKenzie,  1984).  If  partitioning  really  docs  occur  because  it  is 
the  favored  waj  for  crustal  rocks  to  accommodate  large-scale  and  large  magnitude  oblique 
motion,  then  knowledge  of  the  orientations  of  stress  in  the  upper  crust  adjacent  to  such 
faults  will  provide  little  insight  into  the  dynamics  of  the  deformation,  as  pointed  out  by 
Molnar  (1991). 

The  significance  of  the  normal  faulting  earthquakes  in  the  belt  across  the  central 
(Ja.spian  Sea  is  enigmatic.  Eor  roiisons  expressed  above  we  do  not  believe  that  th^y  rep- 
ri„ent  a  southward  motion  of  the  Caspian  Sea  rclaiive  to  Eurasia.  Shortening  at  depths 
shallower  than  the  normal  faulting  is  implied  by  th(!  orientation  ol  folds  in  the  sediment, 
a  possible  (but  arnluguous)  second  sub-event  of  tin;  19S3.03.0G  earthquake  that  involved 
reverse  faulting,  and  a  CJM'k  solution  for  a  small  (and  probably  not  well-c,mstrained)  earth¬ 
quake  that  involved  thrusting  on  1987.12.19  (Eig  3()b).  We  sustjcct  that  this  sliortening 
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does  indeed  represent  a  NNE  motion  of  the  south  Caspian  basin  relative  to  Eurasia,  but 
that  the  motion  is  slow  and  has  not  produced  many  earthquakes.  In  this  respect  the  con¬ 
vergence  may  be  similar  to  that  in  the  Makran  of  SE  Iran  (Jackson  McKenzie,  1984) 
and  along  the  Washington  -  Oregon  coast  of  the  NW  United  States  (Heaton  &  Kanamori, 
1984),  both  of  which  have  rather  low  levels  of  seismicity,  but  may  be  capable  of  generating 
infrequent  large  earthquakes.  Our  preferred  interpretation  of  the  normal  faulting  events  is 
that  they  are  related  to  either  bending  or  shallow  down-dip  extension  of  a  slab  of  mantle 
lithosphere  attached  to  and  beneath  the  Caspian  Sea,  as  it  is  overthrust  by  Eurasia,  but 
this  interpretation  cannot  be  regarded  as  conclusive.  If  the  motion  of  the  Caspian  Sea 
relative  to  Eurasia  is  indeed  slow,  'Jieii  the  motion  in  the  Alborz  between  central  Iran  and 
the  south  Caspian  Basin  will  be  almost  the  same  as  that  between  Iran  and  Eurasia,  as  was 
assumed  by  Jackson  &  McKenzie  (1984)  in  their  analysis  of  the  tectonics  of  Iran. 

We  conclude  that  the  active  tectonics  implied  by  the  observed  faulting  in  the  Talesh, 
Alborz  and  central  Caspian  Sea  will  lead  to  the  eventual  destruction  of  the  south  Caspian 
Basin  and  perhaps  the  formation  of  an  intermediate  depth,  dipping  seismic  zone  in  the 
continental  interior,  similar  to  th.at  presently  observed  in  the  Hindu  Kush. 
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Table  1  :  Earthquake  Hypocentral  Data 


d 

Date 

m 

y 

Origin 

Time 

Latitude 

(“N) 

Longitude 

(°E) 

Depth(km) 
P  pP 

mb 

Ms 

Source* 

Number  of 
Stations 

2 

7 

1957 

00:42:23 

36.14 

52.70 

10 

6.8 

ISS 

1 

9 

1962 

19:20:40 

35.63 

49.87 

27 

7.2 

ISS 

3 

1 

1969 

03:16:37 

37.10 

57.83 

4 

5.4 

ISC 

163 

30 

7 

1970 

00:52:20 

37.85 

55.94 

22 

5.7 

ISC 

306 

14 

2 

1971 

16:27:32 

36.62 

55.74 

4 

5.3 

ISC 

236 

3 

2 

1976 

16:40:41 

39.93 

48.41 

67 

5.2 

ISC 

226 

25 

5 

1977 

11:01:47 

34.91 

52.06 

39 

41 

5.3 

ISC 

258 

4 

11 

1978 

15:22:19 

37.71 

48.95 

37 

36 

6.0 

6.1 

ISC 

418 

1 

10 

1979 

07:37:01 

40.07 

51.87 

42 

47 

5.0 

4.4 

ISC 

150 

9 

12 

1979 

09:12:04 

35.14 

56.87 

23 

5.2 

5.5 

ISC 

248 

4 

5 

1980 

18:35:19 

38.09 

49.07 

36 

36 

5.3 

6.1 

ISC 

323 

22 

7 

1980 

05:17:08 

37.36 

50.35 

37 

36 

5.3 

5.1 

ISC 

238 

3 

12 

1980 

04:26:15 

37.17 

50.47 

44 

22 

5.1 

4.7 

ISC 

233 

19 

12 

1980 

01:16:56 

34.50 

50.67 

29 

43 

5.5 

5.8 

ISC 

354 

22 

12 

1980 

12:51:21 

34.49 

50.67 

39 

18 

5.4 

5.2 

ISC 

299 

4 

8 

1981 

18:35:43 

38.21 

49.41 

40 

29 

5.4 

5.2 

ISC 

358 

19 

11 

1981 

14:10:37 

40.78 

49.26 

42 

42 

5.1 

4.2 

ISC 

236 

26 

3 

1983 

04:07:19 

36.06 

52.28 

33 

25 

5.4 

4.9 

ISC 

261 

22 

7 

1983 

02:41:01 

36.98 

49.23 

43 

33 

5.6 

5.0 

ISC 

299 

22 

2 

1984 

05:44:37 

39.52 

54.11 

0 

5.1 

5.8 

ISC 

296 

29 

10 

1985 

13:13:40 

36.75 

54.81 

13 

8 

6.0 

6.0 

ISC 

518 

6 

3 

1986 

00:05:37 

40.37 

51.60 

28 

63 

6.1 

6.2 

ISC 

614 

11 

6 

1986 

20:15:41 

40.28 

51.71 

41 

49 

5.2 

4.6 

ISC 

305 

10 

4 

1987 

06:43:20 

37.21 

57.70 

6 

35 

5.0 

ISC 

226 

7 

9 

1987 

11:32:27 

39.37 

51.76 

39 

5.5 

5.5 

ISC 

442 

19 

12 

1987 

08:27:36 

40.72 

52.05 

80 

67 

4.9 

ISC 

198 

23 

8 

1988 

05:30:51 

35.42 

52.28 

35 

25 

5.0 

4.8 

ISC 

249 

16 

9 

1989 

02:05:09 

40.34 

51.53 

55 

6.4 

6.5 

PDE 

563 

17 

9 

1989 

00:53:40 

40.20 

51.75 

51 

6.1 

6.1 

PDE 

533 

20 

1 

1990 

01:27:10 

35.83 

52.95 

25 

5.5 

5.9 

PDE 

348 

20 

6 

1990 

21:00:10 

36.96 

49.41 

19 

6.4 

7.7 

PDE 

352 

*ISS  International  Seismologiral  Summaries  (  -  1963) 
ISC  International  Seismological  Center  (1964  -  1988) 
EIDR  Earthquake  Data  Reports  (1989  -  1990) 

I’DE  Preliminary  Determination  of  F^picenters  (1990-  ) 
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Table  2  :  Source  Mechanism  Data 


d 

Depth 

(km) 

Mo 

xlO*’'  Nm 

Mvv 

Strike 

Dip 

Rake 

Slip* 

Vector 

Phases 

P  SH 

1 

9 

1962 

10±4 

367.8 

7.0 

lOr  1  10*' 

52°  ±  3° 

70°  ±  15° 

42° 

10 

4 

3 

1 

1969 

7±4 

1.9 

5.5 

132*’  1  15*’ 

50°  ±  5° 

95°  ±  8° 

33° 

6 

4 

30 

7 

1970 

11±4 

42.9 

6.3 

293**  1  20° 

56°  ±  5° 

210°  ±  10° 

23° 

23 

16 

14 

2 

1971 

11±5 

4.0 

5.7 

336*’ 

39° 

93° 

62° 

9 

4 

11 

1978 

21±5 

19.0 

6.1 

346°  1  15° 

79°  ±  4° 

95°  ±  20° 

IT 

22 

16 

4 

5 

1980 

15±5 

40.2 

6.3 

181°  1  10° 

84°  .+  10° 

267°  ±  15° 

91° 

21 

17 

19 

12 

1980 

14±8 

14.0 

6.0 

115°  1  20° 

41°  ±  8° 

120°  ±  15° 

25° 

10 

11 

22 

12 

1980 

1515 

2.8 

5.6 

113°  1  20° 

56°  ±  8° 

125°  ±  15° 

23° 

4 

8 

4 

8 

1981 

20  +  8 

2.4 

5.6 

154°  1  15° 

35°  +  15° 

32°  ±  20° 

64° 

15 

22 

7 

1983 

1018 

1.9 

5.6 

120°  1  20° 

35°  ±  8° 

83°  ±  15° 

30° 

10 

7 

22 

2 

1984 

2714 

5.1 

5.7 

106°  1  20° 

60°  ±  10° 

174°  ±  20° 

16° 

7 

21 

29 

10 

1985 

1315 

21.8 

6.2 

246°  1  15° 

66°  ±  10° 

71°  ±  10° 

16° 

31 

14 

6 

3 

1986a 

31 

25.2 

6.2 

299°  ±  15° 

88°  ±  5° 

275°  ±  8° 

29° 

35 

15 

6 

3 

1986b 

1314 

21.8 

6.2 

114°  ±  15° 

63°  ±  5° 

87°  ±  8° 

23° 

35 

15 

9 

1987 

3018 

2.3 

5.5 

305°  ±  10° 

10°  ±  10° 

103°  ±  20° 

22° 

14 

12 

16 

9 

1989 

31110 

68.4 

6.5 

80°  ±  15° 

26°  ±  5° 

225°  ±  15° 

38° 

22 

14 

17 

9 

1989 

35110 

21.7 

6.2 

277°  ±  15° 

50°  ±  5° 

249°  ±  15° 

38° 

21 

20 

*Slip  vector  direction  is  for  S  or  E  side  of  fault  relative  to  N  or  W  side 
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Table  3  :  Source  Mechanisms  From  Other  Studies. 


d 

Date 

m  y  Depth  (km)^ 

Ma 

xlO^^  Nm 

M^ 

Fault  Plane 
Strike  Dip  Rake 

Slip 

Vector 

Source* 

2 

7 

1957 

256° 

44° 

259° 

30° 

M72 

3 

2 

1976 

298° 

72° 

184° 

28° 

JM84 

25 

5 

1977 

112° 

60° 

112° 

22° 

JM84 

1 

10 

1979 

35(F) 

0.4 

5.0 

304° 

20° 

-63° 

214° 

HVD 

9 

12 

1979 

15(F) 

2.8 

5.6 

350° 

44° 

121° 

40° 

HVD 

22 

7 

1980 

30 

2.6 

5.5 

310° 

70° 

00 

00 

40° 

HVD 

3 

12 

1980 

16(F) 

1.2 

5.3 

281° 

57° 

47° 

11° 

HVD 

19 

11 

1981 

33(F) 

0.6 

5.1 

131° 

73° 

23° 

41° 

HVD 

26 

3 

1983 

10(F) 

1.2 

5.3 

104° 

61° 

17° 

14/96° 

HVD 

11 

6 

1986 

50 

0.8 

5.2 

291° 

43° 

-123° 

201° 

HVD 

10 

4 

1987 

15 

0.5 

5.1 

292° 

45° 

30° 

22° 

HVD 

19 

12 

1987 

33(F) 

0.3 

4.9 

322° 

36° 

110° 

28° 

HVD 

23 

8 

1988 

15(F) 

0.6 

5.1 

348° 

32° 

-41° 

204° 

HVD 

20 

1 

1990 

33(F) 

10.0 

5.9 

357° 

66° 

172° 

0 

00 

0 

HVD 

20 

6 

1990 

14 

880.0 

7.2 

292° 

88° 

-4° 

112° 

GW91 

*  (F)  denotes  fixed  depth  for  moment  sensor  solution 
»  M72  first  motion  fault  plane  solution  from  McKenzie  (1972) 

JM84  first  motion  fault  plane  solution  from  Jackson  and  McKenzie(1984) 

HV'D  Harvard  Centroid  Moment  Tensor  Catalogue 

C1VV9I  .Moment  tensor  solution  from  Gao  and  Wallace  (1991) 
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FIGURE  CAPTIONS 


Figure  21,  Distribution  of  epicenters  for  earthquakes  greater  than  magnitude  4  occurring 
in  the  region  surrounding  the  Caspian  Sea  during  the  period  1962-1988.  Epicenters  are 
taken  from  the  NEIS  catalogue. 

Figure  22.  Map  of  the  Caspian  Sea  region  a)  showing  the  1,000,  2,000  and  3,000  m 
topographic  contours,  region  of  ” oceanic-like”  crust,  contours  of  sediment  thickness,  and 
fold  axes  in  sedimentary  cover;  b)  showing  region  of  ” oceanic-like”  crust,  major  fault 
systems,  and  seismicity  (M  >  4)  for  the  period  1962-1988.  Geologic  features  and  location 
of  the  ”oceanic-like”  crust  are  adapted  from  Shikalibeily  &  Grigoriants  (l98l),  Berberiaii 
(1983)  and  Jackson  &  McKenzie  (1984). 

Figure  23.  Lower  hemisphere  equal  area  projections  of  P-wave  first  motion  polarity  data. 
Station  positions  have  been  plotted  with  a  velocity  of  6.8  Arms"*  beneath  the  source,  the 
same  as  used  in  the  waveform  inversions.  Filled  circles  are  compressional  first  motions; 
open  circles  are  dilational;  nodal  onsets  are  marked  with  crosses.  Large  symbols  are 
polarities  read  on  long-period  WWSSN  instruments,  small  symbols  are  polarities  based 
on  clear  short-period  seismograms  whose  onset  displays  the  short-period  impulse  response 
(see  Jackson  <k  McKenzie,  1984).  P-  and  T-  axes  are  denoted  by  the  letters  P  and  T.  Nodal 
planes  are  those  of  the  minimum  misfit  waveform  inversion  solutions.  Above  each  focal 
sphere  is  the  event’s  date  (year,  month,  day).  Comments  on  individual  first  motion  polarity 
plots  accompany  the  discussion  of  each  event  in  the  appendix,  except  for  1978.11.04  and 
1986.03.06  which  are  discussed  in  the  main  body  of  the  text. 

Figure  24.  The  P  and  SH  radiation  patterns  of  the  minimum  misfit  solution  for  the  4 
November  1978  (1978.11.04)  earthquake  are  shown  in  the  center  of  the  upper  and  lower 
portions  of  this  figure,  respectively,  and  the  corresponding  values  of  strike,  dip,  slip  (in 
degrees),  centroid  depth  (in  km),  and  seismic  moment  (in  units  of  10**Nm)  are  listed 
beneath  the  event  header.  The  focal  spheres  are  shown  with  the  P  and  SH  nodal  planes,  in 
lower  hemisphere  projections.  The  P  and  T  axes  are  marked  by  the  solid  and  open  circles, 
respectively.  Surrounding  the  focal  sphere,  the  observed  P  and  SH  waveforms  (solid  lines) 
are  compared  with  synthetic  waveforms  (dashed  lines)  computed  for  the  minimum  misfit 
solution.  For  the  waveform  comparison,  amplitudes  have  been  normalized  to  that  of  an 
instrument  with  a  gain  of  3000  at  an  epicentral  distance  of  40°.  The  time  windows  for  the 
waveform  inversion  are  denoted  by  the  solid,  vertical  lines  on  the  waveform.  The  label  to 
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the  left  of  each  waveform  identifies  the  station  code,  an  upper  case  letter  corresponding 
to  the  location  on  the  focal  sphere,  and  a  lower  case  letter  denoting  the  instrument  type 
(w  =  WWSSN  LP,  d  =  GDSN  LP).  The  source  time  function  is  shown  below  the  P-wave 
focal  sphere.  Below  the  source  time  function  is  the  waveform  time  scale.  The  waveform 
amplitude  scales  are  to  the  left  of  the  focal  spheres. 

Figure  25.  This  shows  the  results  of  the  source  region  velocity  structure  test  for  the  4 
November  1978  event.  To  the  left  of  the  waveforms  are  the  P  and  SH  focal  spheres  and 
source  time  function  corresponding  to  the  minimum  misfit  solutions  for  the  particular  test. 
The  strike,  dip,  rake,  and  centroid  depth  are  listed  above  the  focal  spheres.  The  upper 
row  (A)  of  the  observed  (solid  line)  -  synthetic  (dashed  line)  comparison  is  for  the  velocity 
structure  used  in  determining  the  minimum  misfit  solution  shown  in  Figure  24.  Rows  (B 
-  D)  shows  waveforms  for  minimum  misfit  solutions  obtained  for  velocity  structures  B  - 
D  below;  all  of  which  might  represent  reasonable  crustal  velocity  models.  These  synthetic 
waveforms  are  indistinguishable  from  those  in  Row  A  and  the  source  orientations  have 
hardly  changed.  Rows  E  -  F  compare  minimum  misfit  synthetics  for  more  radical  changes 
in  the  source  velocity  structure,  neither  of  which  is  probably  reasonable  for  the  Caspian 
region.  Row  F  is  for  a  mantle  source,  but  the  centroid  depth  resulting  from  the  inversion 
is  at  the  base  of  the  crust.  We  conclude  from  these  tests  that  the  focal  parameters  strike, 
dip,  and  rake  are  insensitive  to  reasonable  variations  in  the  velocity  structure. 


Velocity 

Depth  To 

density 

Structuie 

Interface  (kin) 

(ferns') 

(ferns'  ') 

kgm-3 

A 

0.0 

5.6 

3.2 

2600 

12.0 

6.8 

3.9 

2900 

B 

0,0 

6.2 

3.6 

2700 

13.2 

6.8 

3.9 

2900 

C 

0,0 

6.2 

3.6 

2700 

18.0 

6.8 

3.9 

2900 

D 

0.0 

6.5 

3.8 

2800 

E 

0.0 

5.6 

3.2 

2500 

F 

0,0 

6.2 

3.6 

2700 

25.0 

8,1 

4.7 

3300 
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Figure  26.  This  figure  summarizes  the  tests  we  have  made  to  cissess  the  uncertainties 
for  the  strike,  dip,  rake,  and  centroid  depth  for  the  4  November  1978  earthquake.  The 
format  for  this  figure  is  the  same  cis  for  Figure  25.  Row  A  compares  selected  waveforms 
for  the  minimum  misfit  solution  shown  in  Figure  24.  Row  B  shows  the  comparison  for 
the  minimum  misfit  solution  obtained  when  the  depth  is  held  fixed  at  13  km  but  all  other 
parameters  are  free.  Rows  C,  D,  and  E  show  the  similar  comparison  except  when  the 
strike,  dip,  and  rake,  respectively,  are  held  fixed  at  the  values  shown  after  the  source  time 
function.  The  rake  for  this  event  is  poorly  constrained. 

Figure  27.  The  minimum  misfit  solution  for  the  6  March  1986  earthquake  in  the  central 
Caspian  Sea.  The  figure  format  is  the  same  tis  that  for  Figure  24.  a)  single  source  solution; 
b)  double  source  solution. 

Figure  28.  This  figure  compares  observed  waveforms  and  synthetic  waveforms  computed 
for  various  inversion  solutions  for  the  1986.03.06  earthquake.  The  figure  format  is  the  same 
as  for  Figure  25.  Row  A  compares  waveforms  for  the  single  source  inversion  solution;  row 
B  compares  waveforms  for  the  double  source  inversion  solution.  Note  the  improvement  in 
the  fits  for  the  four  nodal  P-wave  seismograms  (SHL,  KOD,  PTO,  STU).  Rows  C  and  D 
show  the  contribution  of  the  two  individual  sources  of  the  double  source  solution. 

Figure  29.  Comparison  of  waveforms  at  same  common  seismograph  stations  for  the  three 
largest  central  C8Lspian  Sea  earthquakes  (1986.03.06,  1989.09.16,  1989.09.17).  Dotted  lines 
denote  the  synthetic  seismograms  computed  for  these  from  the  minimum  misfit  solutions 
shown  in  Figures  27b,  A 13,  A 14  and  given  in  Table  2.  The  amplitude  scale  (d-GDSN, 
w-WWSSN)  is  to  the  left  of  each  row  and  all  seismograms  have  a  common  time  scale.  All 
ststions  are  WWSSN  except  GDH  which  is  a  GDSN  station. 

Figure  30.  Minimum  misfit  fault  plane  solutions  of  the  earthquakes  analyzed  in  this  paper 
(Table  2)  (a),  and  fault  plane  solutions  for  other  earthquakes  that  we  did  not  study  (Table 
3),  from  Jackson  &:  McKenzie  (1984),  Gao  A:  Wallace  (1991),  and  the  Harvard  CMT  cat¬ 
alogue  (b).  Focal  spheres  are  lower  hemisphere  equal  area  projections,  with  compressional 
quadrants  shaded.  Next  to  each  sphere  is  the  date  of  the  event  (year,  month,  day).  In 
(a)  the  centroid  depth  in  km  is  noted  near  the  epicenter.  The  focal  mechanism  of  the 
1990.06.20  event  is  denoted  by  the  larger  fault  plane  solution  in  10b  and  its  depth  is  noted 
near  the  epicenter. 
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Figure  31.  Horizontal  projection  of  the  slip  vectors  for  earthquakes  whose  focal  mechanisms 
are  shown  in  Figure  30.  Heavy  solid  arrows  are  slip  vectors  from  events  whose  P-  and  SH- 
waveforms  we  analyzed  in  this  study  (Fig.  30a)  and  from  the  1990.06.20  earthquake  (Gao 
(fc  Wallace,  1991).  Light  open  arrows  are  slip  vectors  from  first  motion  solutions  of  Jackson 
&  McKenzie  (1984)  or  from  Harvard  CMT  solutions  (Fig.  30b).  The  sense  of  slip  shown  is 
of  the  south  or  east  side  of  the  fault  relative  to  the  north  or  west  side.  We  have  shown  the 
slip  vector  implied  by  our  minimum  misfit  solution  for  1971.02.14  and  labeled  it  62°  (?)  to 
indicate  that  we  do  not  have  a  great  deal  of  confidence  in  this  solution.  The  inset  in  the 
upper  right-hand  corner  of  the  figure  shows  the  orientation  of  the  T  axes  of  the  normal 
faulting  events  in  the  central  Caspian  Sea 
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APPENDIX 

In  this  appendix  we  show  the  observed  seismograms  for  the  events  not  discussed  in 
detail  in  the  main  text,  and  compare  them  with  synthetic  seismograms  computed  for  the 
minimum  misfit  inversion  solution.  For  each  one  of  these  events  we  have  conducted  a  series 
of  tests  identical  to  those  shown  in  Figure  26  to  place  bounds  on  the  strike,  dip,  slip,  and 
centroid  depth.  Focal  parameters  for  the  minimum  misfit  solution  from  the  inversion  and 
bounds  on  these  parameters  from  the  test  are  given  in  Table  2.  The  format  of  each  figure 
is  the  same  as  in  Figure  24. 

Figure  Al.  1962  September  1  -  The  Buyin  Zara  earthquake  {Ms  7.2,  mj  6.9)  devastated 
the  area  south  of  Qazvin  in  northern  Iran,  killing  around  12,200  people.  Faulting  was 
discontinuous,  but  was  observed  in  a  zone  striking  103°  over  a  distance  of  about  80  km 
(Ambraseys  1962;  Ambraseys  Melville  1982).  Displacements  generally  involved  both 
thrusting  on  planes  dipping  south  and  also  left-lateral  strike  slip.  Fault  segments  typically 
had  a  strike  of  around  100°.  Average  amplitudes  were  about  1.4  m  vertical  displacem.ent 
and  0.6  m  strike-slip. 

Many  of  the  seismograms  of  this  earthquake  were  off-scale,  and,  because  it  occurred  in 
the  early  days  of  the  WWSSN,  not  many  seismograms  were  available  anyway.  Those  we 
were  able  to  obtain  are  shown  in  Figure  Ala.  The  coverage  of  the  focal  sphere  is  not  very 
good,  with  most  stations  in  the  NW  and  only  two  in  the  SE.  The  first  motions  are  better 
distributed  (Fig.  23),  and  require  a  large  reverse  slip  component.  To  obtain  the  minimum 
misfit  solution,  shown  in  Figure  Ala,  we  used  only  three  SH  waveforms:  at  BAG,  PMG, 
and  STU.  We  gave  the  PMG  SH  waveform  only  half  the  weight  it  would  normally  have, 
because  of  its  distance  (101°),  at  which  the  ScS  phase  might  be  expected  to  interfere.  We 
justified  its  retention  because  of  its  similarity  to  the  SH  waveform  at  BAG  (distance  65°, 
at  which  ScS  does  not  interfere  with  S).  We  did  not  use  the  MDS  SH  waveform  in  the 
inversion  because  of  its  distance  (93°),  and  because  we  had  no  independent  check  on  the 
ScS  contamination;  it  is  shown  because  the  polarity  of  its  onset  suggests  an  SH  nodal  plane 
between  its  position  and  that  of  STU  on  the  focal  sphere.  By  comparing  the  amplitude 
of  the  seismograms  at  BAG  with  those  of  other  stations,  we  deduced  that  its  gain  was 
erroneously  given  zis  1500  on  the  WWSSN  film  chip,  and  that  it  should  have  been  3000: 
the  inversion  was  not  sensitive  to  this  change,  which  we  have  corrected. 

The  fit  to  observed  waveforms  in  Figure  Ala  is  reeisonable,  and  demonstrates  a  long 
time  function  of  about  24  seconds.  At  a  rupture  velocity  of  3  kms~^  this  can  account  for 
72  km  of  faulting,  which  is  compatible  with  what  was  observed  (Ambraseys  &;  Melville 
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1982).  The  moment  of  3.68X10’^  Nm  {Mw  7.0)  can  account  for  1.4  m  of  displacement  on  a 
fault  72  km  long,  extending  to  a  depth  of  10  km  with  a  dip  of  52°.  This  amount,  its  sense 
(which  is  predominantly  reverse  with  a  smaller  left-lateral  component),  and  its  strike,  axe 
also  compatible  with  surface  observations. 

We  were  particularly  interested  in  the  slip  vector  of  this  earthquake,  and  so  carried  out 
tests  on  the  resolution  of  the  rake,  which  controls  the  strike  of  the  auxiliary  (north-dipping) 
plane,  and  hence  slip  vector.  These  tests  are  shown  in  Figure  Alb. 

The  minimum  misfit  solution  (rake  =  70°)  is  shown  in  row  2.  In  rows  1,3  and  4  the 
rake  was  fixed  at  some  other  value,  and  the  other  parameters  were  allowed  to  change  in  the 
inversion.  The  rake  cannot  have  a  value  below  60°  without  producing  the  wrong  polarity 
for  SH  at  STU,  thus  limiting  the  maximum  left-lateral  component.  The  minimum  left- 
lateral  component  is  less  well  resolved:  the  SH  fit  at  STU  improves  if  the  rake  exceeds  70°, 
and  a  value  of  100°  produces  a  better  fit  than  the  minimum  misfit  solution  (row  4):  but  it 
requires  a  right-lateral  component,  which  is  incompatible  with  surface  observations.  The 
minimum  misfit  slip  vector  is  42°,  but  we  conclude  it  could  have  a  value  between  30°  and 
50°  (corresponding  to  a  rake  of  70°  —  10/  -f  15°).  Note  that  there  is  some  trade-off  between 
strike  and  rake  in  this  case.  The  dip,  depth,  and  moment  change  little  when  the  strike  and 
rake  are  perturbed.  We  estimate  the  mechanism  of  this  event  to  be:  strike  101°  ±  10°;  dip 
52°  ±  3°;  rake  70°  -  10°/  +  15°;  depth  10  ±  4  km. 

Figure  A2.  1969  January  3  -  P-wave  first  motion  polarities  (Fig.  23)  require  high  angle 
reverse  faulting,  but  cannot  constrain  the  strike  and  rake  well.  Waveform  inversion  and 
testing  constrain  the  focal  depth  to  7  ±  4  km.  The  SH  waveforms  at  SHL,  CHG,  and  NAI 
constrain  the  strike  to  132°  ±  15",  dip  50°  ±  5°,  and  rake  to  95°  ±  8°. 

Figure  A3.  1970  July  30  -  P-wave  first  motion  polarities  for  this  event  were  studied  by 
Jackson  and  McKenzie  (1984).  They  plotted  first  motions  of  this  event  with  a  mantle  source 
velocity  of  Up  8.1  kms~^  because  they  could  not  fit  orthogonal  nodal  planes  otherwise.  We 
suspect  that  some  of  the  polarities  at  European  stations,  which  they  gave  as  dilatational, 
are  in  fact,  nodal.  With  these  changes,  the  nodal  planes  for  the  minimum  misfit  solution 
are  consistent  with  the  first  motion  polarities  plotted  with  the  focal  velocity  of  6.8  kms~^. 
We  estimate  the  centroid  depth  to  be  11  ±  4  km. 

Figvire  ,\4.  1971  February  14  -  P-wave  first  motions  and  P-wave  waveform  inversion  require 
this  event  to  result  from  high-angle  reverse  faulting.  First  motion  data  suggest  a  nodal 


108 


plane  roughly  as  shown.  P-waveforms  constrain  the  centroid  depth  to  11  ±  5  km  but  do 
not  greatly  aid  in  constraining  other  focal  parameters,  particularly  strike  and  slip  vector 
azimuth.  No  suitable  SH  waveforms  were  available  for  this  event. 

Figure  A5.  1980  May  4  -  First  motion  polarities  require  this  event  to  have  one  almost 
vertical  nodal  plane  and  this  is  confirmed  by  the  waveform  inversion.  Inversion  and  testing 
constrain  the  centroid  depth  to  15  ±  5  km,  the  strike  to  28°  ±  10°,  the  dip  to  7°  ±  10°,  and 
the  rake  to  297°  ±  15°. 

Figure  A6.  1980  December  19  -  First  motion  P-wave  polarities  require  high  angle  reverse 
faulting,  but  do  not  constrain  the  strike  and  rake  well.  Waveform  inversion  and  testing 
constrain  the  centroid  depth  to  14  ±  8  km,  the  strike  to  115°  ±  20°,  the  dip  to  41°  ±  8°, 
and  the  rake  to  120°  ±  15°. 

Figure  A7.  1980  December  22  -  First  motion  P-wave  polarities  require  high  angle  reverse 
faulting,  but  do  not  constrain  the  strike  and  rake  well.  Waveform  inversion  and  testing 
constrain  the  centroid  depth  to  15  ±  5  km,  the  strike  to  113°  ±  20°,  the  dip  to  56°  ±  8°,  and 
the  rake  to  125°  ±  15°.  The  first  motion  at  NAl  is  incompatible  with  the  planes  as  shown 
(Fig.  23)  but  lies  very  close  to  a  nodal  plane  and  within  the  acceptable  error  in  its  dip. 

Figure  A8.  1981  August  4  -  First  motion  data  constrain  the  steeply  dipping  plane  (Fig. 
23).  Two  first  motion  polarities  (JER,  1ST)  are  inconsistent  with  the  minimum  misfit 
solution  (strike  154°  ±  15°,  dip  35°  ±  15°,  rake  32°  ±  20°).  Both  are  near  a  nodal  plane  and 
are  at  regional  distances,  so  that  their  precise  position  on  the  focal  sphere  is  uncertain. 
For  comparison  with  our  results,  the  Harvard  CMT  solution  gives  strike  159°,  dip  26°,  rake 
40°,  and  centroid  depth  of  25  km  for  this  event,  which  is  within  the  uncertainty  bounds, 
identical  to  our  solution. 

Figure  A9.  1983  July  22  -  This  event  occurred  one  minute  prior  to  the  1983  Coalinga 
earthquake  (Mw  5.75)  in  California,  whose  seismograms  caused  interference  at  some  sta¬ 
tions.  We  chose  seismograms  from  stations  where  the  Coalinga  P-  or  SH-  phases  do  not 
arrive  within  the  inversion  window.  P-wave  first  motion  polarities  (Fig.  23)  require  high 
angle  reverse  faulting,  but  do  not  constrain  the  strike  and  rake  well.  Waveform  inversion 
constrains  the  strike  to  120°  ±  20°,  the  dip  to  35°  ±  8°,  and  the  rake  to  83°  ±  15°.  The 
nodal  planes  shown  in  Figure  23  are  for  the  minimum  misfit  solution  and  violate  the  first 
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motion  data  for  JER  and  1ST.  Both  of  these  stations  are  at  regional  distances  and  their 
precise  position  on  the  focal  sphere  is  uncertain. 

Figure  AlO.  1984  February  22  -  We  could  not  read  any  reliable  first  motion  polarities  for 
this  event.  Only  a  limited  number  of  P-waveforms  are  available  and  all  but  one  of  these 
are  from  the  GDSN.  The  mechanism  (strike  106“  ±  20“,  dip  60“  ±  10“,  rake  174“  ±  20“)  is 
primarily  constrained  by  the  SH-waveforms.  The  depth  is  27  ±  4km.  The  Harvard  CMT 
solution  is:  strike  109“;  dip  83“ ;  rake  175“,  and  centroid  depth  40  km.  Except  for  the 
depth,  this  is  indistinguishable  from  our  waveform  inversion  solution. 

Figure  All.  1985  October  29  -  First  motion  polarities  (Fig.  23)  require  this  event  to  have 
a  relatively  steep  nodal  plane  dipping  NW.  The  source  orientation  (strike  106“  ±  15“,  dip 
30“  1  10“,  rake  126“  ±  10“)  is  well  constrained  by  the  abundant  P-  and  SH-waveforms.  The 
depth  is  13  ±  5km.  The  Harvard  CMT  solution  for  1985.10.29  is  strike  97“,  dip  31“,  rake 
72“ ,  and  centroid  depth  15  km.  The  strike  and  dip  of  the  CMT  solution  is  indistinguishable 
from  our  waveform  inversion  solution;  however,  the  rake  is  significantly  different.  We  are 
more  confident  in  our  waveform  inversion  results. 

Figure  A 12.  1987  September  7  -  First  motion  polarities  constrain  a  WNW  striking,  steeply 
dipping  nodal  plane.  Waveform  inversion  and  testing  constrain  the  centroid  depth  to  30  ± 
8  km,  the  strike  to  305“  ±10“,  dip  10“ ±10“,  and  rake  103“  ±20“.  The  Harvard  CMT  solution 
for  1987.09.07  is  strike  312“,  dip  14“,  rake  106“,  and  centroid  depth  29  km,  indistinguishable 
from  our  waveform  inversion  solution. 

F’igure  A 13.  1989  September  16  -  First  motion  polarities  imply  normal  faulting  for  this 
event  but  poorly  constrain  the  mechanism.  Waveform  inversion  for  a  single  source  gives 
a  mechanism  involving  either  faulting  on  a  WNW  striking,  steeply-dipping  plane,  or  on 
a  ENE  striking,  shallow-dipping  plane.  Because  of  the  similarity  of  the  waveforms  of 
this  event  and  1986.03.06  (Fig.  29),  we  also  inverted  waveforms  of  this  event  for  a  second, 
delayed  source.  This  cannot  be  a.s  definitive  as  in  the  case  of  1986.03.06  because  the  GDSN 
instrument  response  tends  to  blur  source  details  and  because  of  the  lack  of  WWSSN 
stations  which  have  a  higher  frequency  response.  The  inversion  with  a  second  source 
produced  marginally  better  fit  to  the  waveforms  than  the  single  source  solution,  and  gave 
a  moment  for  the  second  source  that  was  approximately  10%  of  the  moment  of  the  first 
event.  Our  preferred  single  .source  solution  from  the  waveform  inversion  is  strike  80“  ±  15", 
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dip  26°  ±  5°,  rake  225°  ±  15°  and  depth  31  ±  10km.  For  comparison,  the  Harvard  CMT 
solution  is  strike  104°,  dip  36°,  rake  228°,  and  centroid  depth  34  km.  The  strike  and  dip 
of  the  CMT  solution  is  slightly  different  from  our  waveform  inversion  solution,  which  we 
have  more  confidence  in. 

Figure  A 14.  1989  17  September  -  Waveforms  of  this  event  are  similar  those  of  to  the 
1989.09.16  earthquake.  First  motion  polarities  imply  normal  faulting  but  poorly  con¬ 
strain  the  orientation.  Waveform  inversion  for  a  single  source  yields  a  mechanism  that 
corresponds  either  to  faulting  on  a  WNW  striking,  steeply  dipping  plane,  or  on  a  ENE 
striking,  shallow  dipping  plane.  Figure  29  compares  P-waveforms  of  this  event  and  the 
1986.03.06  and  1989.09.16  earthquakes  at  common  stations.  Because  of  the  scarcity  of 
WWSSN  seismograms,  we  have  only  inverted  this  event  assuming  a  single  source.  The 
resulting  orientation  is  strike  128°  ±  15°,  dip  44°  ±  5°,  and  rake  295°  ±  15°  and  depth  35 
±  10km.  The  Harvard  CMT  solution  for  1989.09.17  is  strike  142°,  dip  58°,  rake  287°, 
and  centroid  depth  34  km.  The  dip  of  the  CMT  solution  is  different  from  our  waveform 
inversion  solution,  however,  we  are  more  confident  in  our  result. 
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